| ors eeRRED PERIODICAL 
be 1 Bee UNIVERSITY OF HAWAIF 
ES ae : i 1 conmpeemmces LIBRARY 


Near 
few 


- _— 

1! a 

i : < 1 ey be ae Sa i 4 ‘ ) * hee 

? 4 ae shi = 

“Ge 5 : Soe? Cr i } ; i ores 

A 3 . 4 ». 
pee, H meet & 2 } a ; om } 

wa } 124 net (ey Ve HAST EU P| 

: } ‘ 7 8 t BRGY ig j . i : J e fy : 4 ) 

| , \ 


Volume ECG . “DECEMBER, 1957 Nuniber 4 


Published Quarterly 


io TABLE OF CONTENTS 
: ee Cad ad yy 
Veer sate Fike 
Ussor CONTRIBUTIONS /¢ / 

The Synthesis and Analysis of Digital Systems by Boolean Matrices..............--- Joseph O. Campeanu 231 
Simulation of Transistor Switching Circuits on the IBM 704..........-..5-..000 seen eee R. J. Domenico 242 
An Optimum Character Recognition System Using Decision Functions............... PROS C.K. Chow 247 

An Analysis of Certain Errors in Electronic Differential Analyzers: I—Bandwidth Limitations........... 
EPCS ee ea RENNER hin rede eh on Glee oy She oS Socairarciralt edges aoa Salsas hehe UI Cs DOW Aly Nae oo 
Synthesis of Vector NOM RES Jia oe ae na tae rae sere Corer R. E. Horn and V.G. Fauque 261 
Switching Functions of Three Variables.......... EU AeA Pie tutes yecteee re altel fo eee Sat a ee ee D. W. Davies 265 
Analvcicvonmoequemtialsliachimes twee 2c ce ecient e+. s es sere ae D.D. Aufenkamp and F. E. Hohn 276 

CORRESPONDENCE 
The Logical Combination of Punched Paper Tapes....... 2... +- scene eevee wetness Robert M. Mason 285 
Demonstration of Conditional Stability on an Analog Computer............... A. Nathan and Y. Mahler 287 
On the Use of Redundant Integrators in Analog Computers.......................00-- Norman R. Scott 287 
GTB L dcr s eee aE CC eh Fer SE ig Seor Hie aay onlay 'e s, sek ovens Makin gy oP 0 Gi SOS, eee roel eae aks shone 288 
TRAD = OS Fil eae en a eres So 289 
aT 
EVONTC TOT COUTTS III CCLUCUIN Cane Men ug Sessa ac ee ierng nce hike at asmtrae af Aros, ebay one: Sauce cape Harry D. Huskey 291 
FTE AW Ieee frost oc 9 arth en cys ha eee ee Gy LAGER S SoMe ON ono eee ag viene eS Follows Page 308 
“el gee 


IRE PROFESSIONAL GROUP ON ELECTRONIC COMPUTERS 
The Professional Group on Electronic Computers is an association of IRE members with professional interest 


in the field of Electronic Computers. All IRE members are eligible for membership, and will receive all 
Group publications upon payment of a fee of $2.00 per year, 1957. 


PGEC OFFICERS 
Werner Bucuuorz, Chairman 


W. H. Ware, Vice-Chairman H. W. Norpyke, Secretary-Treasurer 


PGEC ADMINISTRATIVE COMMITTEE 


Term Ending 1958 Term Ending 1959 Term Ending 1960 
J. M. BRooMALL D. C. BoMBERCER R. O. ENpREs 
D. R. Brown W. BucuHo.z W. L. Evans 
W. H. BurKHART J. C. LaPomntE R.F. Jounston 
J. D. Nor H. P. MEssINcER S. R. OLson 
R. THORENSEN W. H. Ware C. W. ROSENTHAL 
L. B. WADEL R. A. Roccenbuck 
: COMMITTEES 
Membership Committee Sectional Activities Committee 
L. C. Norrey, Chairman S. B. Disson, Chairman 
Student Relations Committee Technical Program Committee 
R. W. Mervinie, Chairman W. D. Wincer, Chairman 
Publications Committee Constitution and By-Laws Committee 
N. M. Bracuman, Chairman J. C. LaPointe, Chairman 


rs 
\ 


PGEC EDITORIAL BOARD 


J. P. Nasu, Editor 


Netson M. BLAcHMAN R. E. Mracuer STANLEY RocERS 
M. Rupinorr J. R. WEINER 


IRE Transactions® on Electronic Computers 


Published by the Institute of Radio Engineers, Inc., for the Professional Group on Electronic Computers 
at 1 Fast 79th Street. New York 21, N.Y. Responsibility for the contents rests upon the authors and not 
upon the IRE, the Group, or its members. Price per copy: IRE-PGEC members, $1.20; IRE members, 
$1.80, nonmembers, $3.60. Yearly subscriptions rate: nonmembers, $17.00: colleges and public libraries, 
$12.75. Address requests to The Institute of Radio Engincers, 1 East 79th Street, N.Y. 21, N.Y. 


Notice to Authors: Address all papers and editorial correspondence to J. P. Nash, Missile Systems 
Division, Lockheed Aircraft Corp., 3251 Hanover Street, Palo Alto, Calif. To avoid delay, 3 copies of papers 
and figures should be submitted, together with the originals of the figures which will be returned on request. 
All material will be returned if a paper is not accepted, 


CopvkIGHT © 1958—TuHeE INSTITUTE OF RAptO ENGINEERS, Ine. 


All rights, including translation, are reserved by the IRE, Requests for re- 
publication privileges should be addressed to the Institute of Radio Engineers 


1957 


IES 


A CS 1S 


The Syntl 


IRE TRANSACTIONS GN ELECTRONIC COMPUTERS 


and Analysis of Digital 


23] 


saa by Boolean Vee 
JOSEPH 0. CAMPEAUt 


Summary—In this paper methods are described by which 
Boolean matrices can. be used to synthesize digital systems. The 
matrices offer a means by which the design of such systems can be 
systematized much in the same way as do matrix methods when ap- 
plied to electrical circuit design. They also present a means by which 
the problems of optimum logical design and programming can be 
approached. 


INTRODUCTION 
[: THE DESIGN of electrical circuits the use of 


matrix methods has been found helpful. The prob- 

~-lems-encountered in the design of electrical circuits 

are similar in some respects to those which exist in the 
design of digital systems. 

It is, therefore, not unreasonable to expect that 
matrix methods might prove helpful when applied to 
digital systems.! This has been found true in two ways. 

The first occurs in the design of digital systems where 
the process to be performed is already determined and 
where the only remaining problem is the writing of the 
logic to perform the task. In this type problem the 
matrix methods appear helpful. For example, one prob- 
lem to be worked later in this paper concerns a system 
which is to add 3 or 5 to any number shifted into a bank 
of flip-flops, or multiply this number by 3 or 5, depend- 
ing on the state of two-phase control flip-flops. 

What the matrix methods do when applied to prob- 
lems of the general type of which the above is an exam- 
ple is to make the work involved assume more of a 
“crank-turning” nature than would otherwise be the 
case. 

The second way in which the matrix methods appear 
valuable is in an approach to the design of optimum 
digital systems and also in optimum programming. 

Often, in the design of digital systems, more than one 
method can be determined which will accomplish the 
desired task. One of the methods, however, will require 
less diodes or flip-flops or both than do the other meth- 
ods. It is the task of the person designing a system to 
devise methods to accomplish the desired task which 
will require the least amount of equipment. He tries to 
design an “optimum” system. 

Much in the same way, once a general purpose digital 
computer has been built there remains the job of pro- 


* Manuscript received by the PGEC, August 27, 1956; revised 
hantuseript received, August 28, 1957. W ork oi this paper was carried 
uit while the author was employed at the Hughes Aircraft Co., 


Culver City, Calif. 

+ Litton Industries, Beverly Hills, Calif. 

1 For the original presentation of this material see J. O. Campeau, 
The synthesis and analysis of counters in digital syste ms by 
ae matrices,” Master’s thesis, University of Calif fornia, Los 
Angeles; June, 1955. 


gramming it to perform some specific process. By prop- 
erly arranging the information in the memory, the pro- 
grammer can often perform the program with less orders 
and in less time than he could have with some other 
arrangement. 

By using matrix methods it is possible to reduce the 
problems of optimum design and programming to a sys- 
tem of simultaneous Boolean equations where the un- 
knowns are the elements of the matrix which represents 
the desired optimum design or program. 

In general, it has been found that these equations 
are too voluminous to be solved by the methods at 
hand, so that the results obtained are not of practical 
significance at the present time. 

If, in the future, techniques are worked out which 
allow the solution of these equations, then a technique 
which will allow an analytical approach to optimum 
design and programming will be available. 

In order to conserve space no detailed developments 
will be presented in this paper. A more complete treat- 
ment of the material can be found in the paper refer- 
enced above. 


BooLEAN MATRICES 
Consider the following set of Boolean equations: 
CF Ce ail AT 
4182 + 41x. = V2. (1) 


In the above set of equations the union or logical 
“or” and intersections, or logical “and” operations, are 
referred to rather than sums and products (as will 
generally be the case unless it is stated that arithmetical 
operations are intended). Eq. (1) can be written (ex- 


panded) into minimal polynomials: 
0X12 + 141%. -b 1X 1x -f- 141% == his 
0X1 *2 a. 141% +- 14% -f- Ox 1x2 Vos (2) 


Notice that here all minimal polynomials appear 
multiplied by either a 1 ora 0. 

Finally the coefficients of the minimal polynomials 
could be written along with the variables in the follow- 


ing form: 
, art sl pet | 3) 
= / 3 
OFC ay rapa) ae) 3 : 


The components in the above vectors can only have 
values of one or zero and so the vectors can be called 


“Boolean vectors.” The array of ones and zeros is not 


aos MAE LNANSACTIONS OF I 


It will be 


and is the array of 


a matrix in the ordinary sense of the word. 


called a “Boolean matrix” the co- 
efficients of the ordered minimal polynomials in (2). In 
the the 


equations written as a 


form shown above original set of Boolean 


has been “matrix” equation. 
Now we determine a method which will allow the pass- 
ing from the Boolean matrix (3) to the set of (2). 

In order to do this, a formalization of the method 
which will generate the minimal polynomial form of 
any system of Boolean equations is required. More 
specifically, in order to determine the Boolean matrix, 
the coefficients of the minimal polynomials which make 
up the function are required. 


To begin, a row matrix A with elements a; is defined: 
aS j a 4 ih 3, make 2% 2” (4) 
where the problem under consideration involves 7 vari- 


ables (the above example involved two variables). In 
the case for 7 =2 the matrix A assumes the form 


F Re. 
—~ 
=) 
ae 
bo 
W 

Nene 

“- 
nr 

a 


Then associated with the ordinary row matrix A we 
define another (Boolean) matrix 4 with rows and 2" 
columns with elements aj; such that 


n 


Bees) 2 las, Be Dy acne 21 Ry » -X6) 
i=1 
Here arithmetical addition and multiplication are 
intended and the a;;’s are defined implicitly. This for- 
mula amounts to decimal-to-binary conversion. 
Using (5) and (6) for »=2 the Boolean matrix A 
assumes the form 


OF ie Omet 
A= : (7) 

OOF 1 4 
The 4 matrix defined in (4) and (6) can thenbe 
used to pass from a system of equations such as (1) to a 
Boolean matrix (3). A “multiplication” of a Boolean 
vector x by a Boolean matrix B can be defined in the 


following manner. 
Start with (1) generalized to 1 dimensions: 


Ji( 41,42, °° t=1,---,n. (8) 


The elements of a matrix B are determined from (8) 


by 


as ia) eye 


4: ads worl ao 


bi; = filarj, ao;, ++ + (9) 


5) @nj) 


Gert aoe 


Using (9), for example, the element by derived from 
(1) would be 


bis — Silay, a) = Aa, 0) =) Pet Ore 


It can be seen that (9) provides a means by which 
the array of coefficients in (3) can be obtained from the 
set of (1). In fact, (9) is merely the desired formalization 
of the usual process used to obtain the minimal poly- 
nomials, 


ELECTRONIC COMPUTERS 


December 


Then given the matrix equation 


Bx = y, (10) 


a set of equations consisting of the unions of selected 
minimal polynomials of the form of (2) can be developed 
by using the relationship 


a a5 bi; I (a4 j%4 + Gij%) t= 1,-++,n; (11) 
j=l 
A : ary 
(the = stands for “equal by definition”). 
Eq. (11) will be used as the definition of the “multi- 
plication of a vector x with elements x; (¢=1, ae) 


by a matrix B with elements b,; (¢=1,---, 2, 
f=, 2”).” In (11) the a,,’s are the elements of 
the -1 matrix defined in (4) and (6). The [| term refers 
to intersections over the range of the dummy index k 
while the >> term refers to unions over the dummy in- 
dex 7. 

Some comment is appropriate at this time concerning 
the above equation. The [| terms are each a particular 
minimal polynomial involving the variables x1, x2, +++, 
x,. For example, consider the minimal polynomial gen- 
erated for 7=2 where the problem involves two vari- 
ables x1, and x, (so that m=2). The [] term in (11) for 
j=2 can be expanded to read: 


2 
ibe (ax2%% + Gok.) = (ayer, + @ 12%) (Ga2X2 + GooKo), 
acl 


From (7): 
Pe = i 


a22 = (0), 


Thus the above expression becomes 


Il (ax2%% + Gyot,) = (1-4, + 0-41)(0- a2 -+ 1-2) 
k=l 


ade 


This expression will be recognized as one of the minimal 
polynomials of the variables «x; and x». 

In (11) the minimal polynomials generated are each 
multiplied by the appropriate element of the Boolean 
matrix. It will be recalled that the elements of the 
Boolean matrix are defined to be the coefficients of the 
ordered minimal polynomials making up a set of Boo- 
lean equations. Thus by multiplying the I] terms (the 
minimal polynomials) by the correct coefficient and 
then summing the results (indicated by the >> over the 
j’s) the original equations can be obtained. 

In considering (10) and (11), the question might 
be asked, why use the J] terms (which involve the 
elements of the 4 matrix defined earlier) in order to get 
the required minimal polynomials? PP urthermore, why 
define a special type of multiplication in order to yer und 
the Boolean matrix of a system of equations? For exam- 
ple, given a set of variables x, - » Xn, a vector could 
be defined whose components would be the ordered 
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minimal polynomials of the variables of the variables 
(this vector would have 2" components). Then this 
vector could be combined with the Boolean matrix of 
the system using ordinary matrix multiplication and 
the desired system of equations obtained without using 
the special operation defined above. 

The answer to the above question is that in the sec- 
tions to follow, the Boolean matrix will be used to de- 
scribe the operation of any digital system in which all 
operations are controlled by a clock (that is, any syn- 
chronous digital system). In using the Boolean matrix 
in this way it becomes necessary to consider repeated 
applications of the matrix, since each application of the 
matrix advances the digital system by one clock time. 

For example, the following operation might be de- 
sired: B(Bx). Here the B matrix is applied to the vector 
x. The result Qwhether the special multiplication is used 
or not) is a system of equations. There will be exactly 2 
of these equations, one for each of the variables in the 
problem. The expression Bx is thus another vector with 
n components (just as was the original vector x). 

It can be seen that in order for the second application 
of the B matrix to have any meaning at all a multipli- 
cation which deals with the elements of a vector x and 
generates the minimal polynomials as part of its opera- 
tion is required. 

Thus even though the first Pea teeth of the B 

matrix can be accomplished using the usual definition 
of the multiplication of a vector and a matrix, this is not 
adequate when repeated applications of a Boolean 
matrix are to occur. 

The multiplication defined (11) transforms one 
Boolean vector with 7 components into another vector 
which also has components so that repeated applica- 
tions of any Boolean matrix can be accomplished. 

As a side comment, it is also possible to define the 
multiplication of a vector x by a Boolean matrix B 
using modulo two addition and intersection: 


A ae, n 
Vi + bi LY Cs © as) aN eS 7, 
1 


j= 


The above definition provides a slightly neater for- 
mula. Note that here the summation is taken to mean 
modulo two addition rather than union. It can be shown 
quite easily that for a given Boolean matrix the right 
side of the above equation and of (11) are identical. 

Also it is interesting to note that the Boolean matrix 
for a given system will be identical to the “truth table” 
which can be developed for the same system. 

The reason for this property will be discussed at 


greater length below. 


Properties of Boolean Matrices 


In the previous section the Boolean matrix was in- 
troduced along with the operation which combines 
any Boolean matrix B which is 2X2” in size with a Boo- 
lean vector of 2 elements to produce a set of » Boolean 
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equations. In this section the Boolean identity matrix, 
matrix addition, matrix multiplication, and the Boolean 
determinant and the Boolean inverse matrix will be 
presented. 

All of these things were developed in an attempt to 
parallel as much as possible the type of operations en- 
countered in ordinary matrix operations. 

The Inner Product: In order to facilitate the proofs to 
follow it is convenient to introduce an abbreviation. 
Given two Boolean vectors a and B with components 
a; and 6;, define 


eat (iB; + a.8;). 


a= 


ae = (12) 


This will be called the “inner product” of the two Boo- 
lean vectors 

From the definition of the identity matrix given in 
(4) and (6) it is obvious that if a; and a; are Boolean 
vectors whose components are respectively the jth and 
ith columns of the identity matrix then: 


(aj, a;) =1 
= (0) 


for7 = 
(13) 


Using the inner product (11) can be written in the 
form 


for 7 ¥ 1. 


on 
“* eS bj(x, aj) (14) 
pal 
where y and x are Boolean vectors whose components 
are respectively y; and x; (¢=1, -,n), and 6; and a; 
are vectors whose components are respectively the jth 
columns of the B and identity matrices. 
Matrix Addition: Define the union (addition) of two 
Boolean vectors x and y with elements x; and y; as a 
new vector with components 


Se = Xe t+ yi AOA. seep, (15) 
where 
z=x+y. 
Matrix addition is given as 
i=, eae 
Siig tet sie ie aan ae (16) 


where 0,;, ci; and d,; are the elements of the matrices 
Brande Cy and rand 


Bre 


From the above definition of matrix addition it can 
be shown that 


(BA GC). = Bs Cs, 


but that B(x+y)4#Bx-+By in general where B and C 
are any two matrices and x and y are any two vectors 
with the same number of components. 


Matrix Multiplication: Feqs. (11) or (14) define what 
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is meant by the multiplication of a matrix C by a vector 
x to produce another vector y: 


Cx = y. 


Suppose now that it is desired to multiply the vector 
y by another matrix B: 


By = B[Cz]. 


This operation can be performed on the vector Cx = y 
using (14) again. 

What is now required is a method of combining only 
elements of B and C in such a way as to produce a new 
matrix D with the property that 


B(Cx) = Dz, 


When this operation (whatever it is) has been developed 
it can be called “Boolean matrix multiplication” since 
it performs exactly the same function with Boolean 
matrices as does ordinary matrix multiplication with 
matrices. The above equation would then further be 
written: 


BlCx| = Da = [BC lx. 


From (4) and (6) it can be seen that no two columns 
of the A matrix are identical element for element. Con- 
sider any Boolean vector x. For each of the 2" possible 
x’s in m space : 


(x, aj) as ie > 


for one and only one j if the a,;’s are the column vectors 

of the A matrix. This fact will now be used in the deri- 

vation of the formula for Boolean matrix multiplication. 
Use (14) to.express the matrix equation 


B[Cx] = z 
in vector form: 
9” 
of ys RGR, aj) 
j=1 
pie Be 
ox +S al > r(x, a,), oi 

et sE=a' r=1 


What will now be done is to show that all operations 
to be performed on the above equation are valid for 
any vector x in ” space, and hence are valid in general. 

It was mentioned above that only one column of the 
A matrix will match element for element with any par- 
ticular vector x. Call this column a». In the above equa- 
tion (x, a,) is nonzero only for r=. Hence all the terms 
for which r=p are equal to zero and can be dropped 
out without affecting the validity of the equation: 


ae =, b;(Cp, aj). 


j=1 


{ if ser 
0 if rx pS’ 


Now since 


(a,, *) = 


December 


the right side of the above equation can be multiplied by 
(a», *) and the terms 


on 


| tiles) | a) 
r=1 j=1 
rep 


added without affecting its value. 


ot 


y=[ Sheed fora +E] E bleed ona 


r=1 j=1 
TAP 


9” 


=D] Leste, on (on, 


rs 


Now define a matrix D with column vectors d, such 
that 


9” 
d; = >) bc, a) ret des oe (17) 
j=1 


Then the above equation becomes 
g” 
Eee es d,(a,, x) 
r=1 


and this is equivalent to the matrix equation 
Dx = ¥, 
where 
B(Cx) = (BC)x = Dx = y, 


Eq. (17) is the formula for Boolean matrix multipli- 
cation of any two matrices B and C to produce the prod- 
uct D. In terms of the elements of the matrices: 


2 
dg a 


Jak 


OFF be (CerOnz + Cures) At 

k=1 ly = ey 

In the above formula the a;;’s are the elements of the 
A matrix defined in (4) and (6). 

The Inverse Matrix and Other Properties: It can be 
shown that the A matrix defined in (4) and (6) is the 
identity matrix. Thus 


Av.= "4, 


for any vector x. Furthermore, it is possible to define 
characteristic vectors of any Boolean matrix as being 
those vectors which for a given matrix have the property 
that 

Be = e 


where é is a characteristic vector of the matrix B. 

The operations of multiplication of a vector or a 
matrix by a scalar are completely analogous to the usual 
operations and the formulas are the same. 

It can be shown that the inverse of a given matrix 
B with column vectors b; has an inverse B~ (if it exists 
at all) with column vectors b,-! given by the following 
formula: 


on 
aor 
r 


= eal rar aj(bj, a) r= 1,++s, Qn 
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where a; and a, are column vectors of the identity matrix 
defined in (4) and (6). 

The scalar expression “1/|B|” which multiplics the 
vector sum is developed so that it resembles the cor- 
responding expression with ordinary matrices. 

First define “division” of two Boolean variables a/b 
as being the problem: given a=bc and knowing a and 8, 
what 1s the value of c? Then define the Boolean determi- 
nant of a Boolean matrix B with column vectors bj; as 


g” on 


| B| = if yy (bj, ai) 


i=1 j=1 


where the a,’s are the column vectors of the identity 
matrix. 

It can be shown that if a Boolean matrix B has an 
inverse that |B] =1; otherwise |B] =0. Thus the ex- 
pression “1/|B|” in the above formula for the inverse 
will have value 1/1=1 if an inverse exists and value 
1/0 =undefined if no inverse exists. 


The Row Boolean Matrix 


Given any Boolean matrix it is possible to define a 
row matrix. The row matrix is obtained by regarding 
each column of the Boolean matrix as a binary number. 
and by writing down the decimal! equivalent to this 
number. A matrix B with elements 0;; becomes a row 
matrix B with elements 6; where 

n ; = 
b; = be 21D a) F 1, hea 2”, (18) 
i=l 

An example will illustrate this. Consider the Boolean 
matrix presented in (3). This matrix and its row equiva- 
lent matrix have been shown below: 


Qantas tant 
ja 
Opidiee iano 
Pen eee sae e | (19) 


In addition to this, given a Boolean vector x it is 
possible to reduce it to a single number by the same 
technique. For example, the vector 


1 
x= | | becomes x = (1) 
0 
while the vector 
1 
x= | becomes x = (3). 


In the above equations arithmetical multiplication 
and addition are to be used. 


Trt Application oF BOOLEAN MATRICES TO THE 
DEsIGN OF DIGITAL SYSTEMS 
The discussion in this section will be limited to digital 
systems which operate synchronously; that is, systems 
in which all operations are initiated by a clock pulse 
Which appears at a fixed rate. Furthermore, all the 
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storage in the system will be assumed to be accomplished 
by identical delay elements. This does not sacrifice any 
generality since it can be shown that any synchronous 
digital system can be reduced to an equivalent system 
which will have only identical delay elements and logic. 
The total number of binary storage elements in the two 
systems will be the same, and the states assumed by 
the storage elements in the two systems will be identical 
at all times. 


The Representation Boolean 


Matrices 


of Digital Systems by 


The synchronous digital system composed of delay 
elements and logic can be described by a system of 
Boolean equations in the following way. The delay ele- 
ments in the system can be called Q1, Qo, - + +, Qn, and 
the variable +1, for example, can represent the state 
assumed by the delay element Q; during this interval. 
In general x; (¢=1, -- -, m) can be the present state of 
Q;. 

When the clock pulse appears, each of the 1 delay 
elements in the system will sample its input logic, and 
when the next clock pulse appears each delay element 
will have assumed the state, either “1” or “0” which its 
input logic had at the present clock pulse. The state the 
delay element Qi will have assumed when the next 
clock pulse appears can be called y1, and in general 
yi (t=1, - ++, m) can represent the next state of Q;. 

In this manner it is possible to represent the state of 
each delay element in the system at the next clock pulse 
in terms of the states of all the elements in the system 
at the present clock pulse. 

Consider the following digital system shown in Fig. 1. 


Fig. 1—Digital system. 


The boxes labeled Q; and Q, are delay elements both 
with identical delays of 1 bit time. These delay ele- 
ments sample the input logic on the left when the clock 
pulse arrives. When the next clock pulse arrives the 
output state of the delay element will be either a “1” or 
a “0” depending on what the input state was. 

In each box Q;, the output on the upper right is x; 
(¢=1, 2) while the one on the lower right is #;. The ele- 
ment with the “++” performs the logical “and” or inter- 
section operation. 

The system of equations which give the state of the 
delay elements Q; and Q: at the next clock pulse (y1 and 
yo) in terms of their present states (x; and x») is then 


X,-- ve = Vi 


NpXo -- Byv. = Vo. 
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It will be noted that this system of equations is identical 
to (1). It will be recalled that the Boolean matrix for 


this sytem of equations was 
Oh ak) koa Rs 
B= 
ks WS We 8 
and that the equivalent row Boolean matrix is given 
from (18) as 
Ly, 


This row matrix then corresponds to the digital sys- 
It is in- 


Bes th 3 8 


tem whose block diagram appeared in Fig. 
structive to trace through the states the system passes 
into when it is placed in various configurations. In order 
to do this write down the row identity matrix A. Directly 
beneath this write the row matrix for the digital system 
5 

2 3), 


Ni doe gan del 


Then to find what the next state of the system will 
be locate the present state in the matrix A. Directly 
beneath this in B will be found the next state of the 
system. 

For example, directly beneath “3” in the A matrix is 
“1” in the B matrix. From (18), the “3” stands for 
%,=1 and x,2=1. The “1” stands for 4,=1 and y= 
This result can be checked by referring to the ieee 
diagrani of the digital system. _ 

Table I summarizes the operation of the system. 


TABLE I 
Present State Next State 
0 0 
1 3 
2 3 
3 1 


It is possible to draw a picture which describes the 


behavior of the system: 
Za 7 Sy) 
Sy Gi} 
“< 


Thus, if the system is placed in the state “0” (in 
which Qi=0 and Qo=0) it-will remain in that state. 
If the system is placed in the state 2 (in which OQ, =0 and 
Q»=1) it will move to state 3 (Q:=1 and Q,=1), and 
so forth. 

The above example illustrates the use of the row 
Boolean matrix to determine the state diagram of any 
The reason that the row identity 
matrix can be used with the row matrix form of the 
Boolean matrix of the system to determine the state 
diagram of the system is as follows. What is actually de- 
sired in determining the state diagram is a truth table 
of the system. The truth table for a system consists of 


given digital system. 
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all possible states the system can assume listed in a 
column and the corresponding transformed states listed 
opposite in a second column. From the definition of the 
row Boolean matrix given in (4) itis obvious that all pos- 
sible states of the system are present in this special 
matrix. 

Consider (14). This equation gives the matrix equa- 
tion Bx=y in terms of the vectors x and y and the 
column vectors of the B matrix and the identity matrix. 
It can be seen that for a given x the expression (x, a,) is 
nonzero for one and only one value of j, say 7=p. Thus 
essentially (14) for a particular x reduces to 

y= b, 


here (x, a;) is nonzero only for 7=p. When the row 
lanes. matrix is written above the column B matrix 
it can be seen that this is equivalent to mechanizing 
(14). Thus for a given x we search through the identity 
matrix until a column is determined which matches x 
exactly (equivalent to determining the number p de- 
scribed above). The pth column of the matrix B then 
becomes the new state of the system, y. 


The Synthesis of Digital Systems 


In this section methods which allow the use of Boolean 
matrices in the synthesis of digital systems with pre- 
scribed operations will be introduced. The Boolean 
matrices offer a method of synthesis which deals only 
with the coefficients of the minimal polynomials making 
up the logic rather than with the minimal polynomials 
themselves. This eliminates the need for writing the 
minimal polynomials until the final logic has been deter- 
mined. ; 

The first problem to be considered will be that of de- 
signing a system which will replace any two numbers 
shifted into it with their product. In the simplest case, 
which is not trivial, the system desired would be one 
which accepts any two binary numbers between zero 
and three. The storage elements comprising the system 
can be called Qi, Q2, Qs, and Qs. Then the least signifi- 
cant bit of the first number will be stored in Q,, the 
most significant bit of the first number in Qs, the least 
significant bit of the second number in Q3, and the most 
significant bit of the second number in Q,. 

In this case, since the storage elements are by the 
nature of the problem grouped two by two rather than 
all together, it is more expedient to write the identity 
row matrix A in the form shown below rather than in 
the usual manner. The least significant bit of the prod- 
uct is to appear in Q; and the most significant bit to ap- 
pear in Qs. Directly beneath each row in A the number 
it is to transform into is written down. The second row 
of numbers can be shown to be the B matrix for the 
system. Then the A and B matrices become 


ete ie 
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Notice that directly below the two elements which 
are present in the modified A matrix is their product. 
For example, in the last column the elements in the A 
matrix are 3 and 3. Directly beneath them in the B 
matrix is 9, the product 3X3. Then using (18) the B 
matrix corresponding to the row matrix B can be deter- 
mined. This amounts to writing beneath each element 
in B acolumn in B in binary form. This matrix B is 
given as 


O20eOsOROrt Ula Us0 US OFT eT 
OF On 020202021. T0210 1021 6) 0) 
00-00-0400! 0202.00. 1el0 021-0 
OF 020206020) 002020 0600.00 1 


Notice, for example, that directly beneath the 9 in the 
B row matrix appears the column[ ! ]in the matrix B. 


0 
0 
1 


Eq. (11) can.then be used to obtain the system of 
equations which specify the inputs to the delay elements 
in the system in terms of the present states of these 
delay elements: 


XX QV 3X4 +. NNN 3X4 -+- Hy Vov3HX4 -+ XiVoV3N4 = V1, 
HpVovghs + xypvough, + eee + ayvotsa4 

+ aptevgry + Fyxoxgry = yo, 

PVN 3A, +- NVR 3Vy + RyVoX 3X, = Ng5 


XiXoN3X4 = Ya. 


These equations can be simplified to read: 


Mins = Vay 

Ko 1e4 Nox gh, - Kovg%s + Xow yL3X4 = Yo, 
Leavers + 3h Kort, = Ys, 

NiXov3Xq = Ya. 


These equations give the next state of each of the 
delay elements in the system in terms of the present 
states of these elements. The present state of the 
th delay element Q; is x; and the next state of this delay 
edd) 

he logic for an equivalent system using what are 
called “JK. flip-flops” can be developed using the 
lormula: 


cloment is y;, @=1,--: 


J jx; -|- Kix; yy = L Ses) Ph TUS (20) 
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The logic using the JK flip-flops is then 


Qi: JQi = 0 
KQi = &3 

Qe: JQ2 = 41%4 : 
KQo = &y€4 + dads + vives 

Q3: JQs = won, 
KQ3 = x1 + #2 + % 

Qu: JQ, = 0 


KO = Sua Sef 23. 

The JK flip-flop is essentially a conventional flip-flop 
with inputs to each of the tube grids. 

As a second example consider the following problem: 
it is desired to design a system with only four delay 
elements Q;, Q2, Q3, and Qs. Any number X =0, 1, 2, or 3 
can be shifted in binary form into the elements Q: and 
Qs. At the same time the delay elements Q3 and Qs are 
also set into one of four states. 

Then, depending ou the states Q3; and Qs were just 
set to, only one of the following things is to occur: 


1) 1f-@;=0, O:=0, replace X by 32% 
2) 1£0,=4, i010; replace.x by 5x4 
3) TiO; =O, Ov=17 replacees bys-os 
4) If Q;=1, Q:=1, replace X by 5+%. 


This problem may be summarized in a block diagram: 


Contxol—_> 


ea 


The steps in the problem solution are as follows: 


1) Write down the modified identity row matrix A. 
2) Develop the row matrix B by writing beneath 
each column of the modified identity matrix the 
desired result. 

Develop the matrix B from B by using (18).This 
amounts to writing down beneath each of the ele- 


3) 


ments of the B matrix its binary equivalent. 


4) From (11) write the logical equations which de- 


—— 


scribe the desired system. Alternately, since the 
elements of the B matrix are the coefficients of 
the ordered minimal polynominals making up the 
desired functions the equations can be written 
directly. 
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These steps have been carried out below: 
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and the a,,;’s are. the elements of the identity matrix. 


fors # 1 


ie! 
Pi ae ry Rin 


m=1 


= 0 for: = 1. (23) 


Now consider the matrix A. Each column of A will 
represent some combination of the inputs to the system, 
the @;’s. Some of the combinations represented by the 
columns of A will never occur since they may contain 


Xe eee 
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X ys eOC8e se RELCT ee Gi ae PEELE. ae X1X2X3X4 
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X3Na(X1 + Xe) + XsXa(X1 + Xo) + XiXsX_ 4+ XiX2NXs = Vs, 


It is possible to formalize the synthesis procedure 
which has been carried out above. Say the inputs to the 
system are the set of numbers &, - ++, @, and that the 
outputs are a second set of numbers by, - ++, 54. Say 


A 
A oo b 
ALT et ae | z — ‘ 
i is DIGITAL re aes 
pe rity = A 
A SYSTEM 
“Saat esate Pires 
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that each of the 4;’s takes on a total of g; values the 4,’s. 
For each @; determine a number k; such that 


2h) = (2.3) max, <* 2*i, 


(21) 


Even though a particular input @; may not take on all 
values 0 to 2‘ all values in this range will be considered 
in the discussion to follow. 

Furthermore, it will be assumed that whatever the 
G;’s and $;’s represent, some code will have been used al- 
ready to make them all positive integers. This does not 
sacrifice generality since it is common practice to em- 
ploy such codes in all digital systems in the sense that 
all quantities are always represented by collections of 
z’s and o’s internally. 

Since each of the inputs will be allowed to take on a 
total of 24* (¢=1, ---, p) values, and since there are a 
total of p inputs, the 4,’s, there will be a total of 


gl Des possible combinations of the @,’s, although 
some of the combinations will not be allowed (in the 
sense that cach @; may not take on all values from 0 
to 2**) in the actual system. 

From an identity matrix with >>?_, k; rows as defined 
in (4) and (6) develop a second matrix A with elements 
a;; Where 


1 = 


i. 5 Do ke 
t= 


(22) 


Ht 


P 
Dene See 
i=1 


Xie -+ X1X9X3 = We 


columns the “unallowed” columns in the discussion to 
follow. 

Now form a second matrix B directly beneath A as 
follows. Each column of A represents some combination 
of the input variables, the @,;’s. Beneath each column of 
A form a column of B where the element b,; will be 
the appropriate value of the particular output variable 
é,. The order in which the 6,’s are arranged in a particu- 
lar column is unimportant, but the same ordering must 
be used for all columns of B. 

Directly beneath unallowed columns of A any ele- 
ments at all may be placed in the matrix B. This is to 
say, these elements are arbitrary. 


ny 

The matrix B will have 2 2 columns and q 
rows (since there were to be g output quantities). In 
each row of B there will be some largest b,;. For each - 
row determine an m; such that 


is ae < (Vata << Le (24) 


Now using m; determined above, form a matrix B 
from the matrix B. The elements of B, b,;, are deter- 
mined from the elements of B as follows: 


qitmi 


bz = a) UN ih es 


n= Ti 


(25) 
where 
i 
= > for: 4 1 
r= 


== ()) LOY e—* 


and where 
b,; = 1 or O only. 


Then B is the Boolean matrix of the desired logic. 
Those columns of B which correspond to “unallowed” 
columns of the matrix A may have any value. These 
columns can be used to simplify the logic which is to 
be developed. 
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Using (11) form the logic for the desired digital sys- 
tem. As was mentioned above, the unallowed columns 
of B, and correspondingly, the minimal polynomials 
developed from them, 
1 or a 0, depending on which produces the simplest 


can be multiplied by either a 


logic. 

Note that the above discussion also shows how any 
synchronous digital system can be represented by a 
Boolean matrix. 

Given a synchronous digital system (in which all 
operations are initiated by a single pulse which comes 
at some nominal repetition rate), the Boolean matrix for 
the system can be determined as follows: 


1) In the same manner as was described above, use 
(21)—(23) to determine a matrix A for the system. 
Beneath each element in the matrix A place the 
corresponding element in the B matrix. The ele- 
ment in B is determined by subjecting the system 
to the particular set of inputs represented by each 
element of A. 

Using (25) develop the Boolean matrix for the 
system. 


2) 


3) 


It should be mentioned that in doing the above opera- 
tions all storage elements in the digital system (includ- 
ing cells on a drum, for example) must be censidered as 
part of the Boolean vector which gives the state of the 
system at any time. 

Cells on the surface of a drim can be ccncilered to be 
part of shifting registers made up of unit delay ele- 
ments. 


THE PROBLEMS OF OpTIMuUM DESIGN AND 
PROGRAMMING 


It was mentioned in the introduction that Boolean 
matrices offer an approach to the problem of optimum 
design and programming. In that section the meaning of 
optimum design and programming was outlined. Essen- 
tially, what is meant is that to perform a given task it 
is always desired to minimize the time taken by the 
equipment, the number of storage elements (flip-flops), 
and the number of logical elements (diodes, for exam- 
ple) required. In the case of the optimum program it is 
desired to so organize the information in the program 
that the problem can be worked in the minimum amount 
of time. 

The material to be presented in this section is an 
outgrowth of material presented earlier. A complete 
treatment of this subject will not be attempted at this 
time, but will be deferred until a later date. This dis- 
cussion will serve as an introduction to the concepts 


involved. 


Optimum Design 


In the previous section on the design of digital systems 
using Boolean matrices it was shown that any syn- 
chronous digital system could be represented by a single 
Boolean matrix. 


of Digital Systems by Boolean Matrices ao 

Suppose that a design has been worked out for a 
given digital system which will perform a given task. 
After converting all parts of the system to delay ele- 
ments and logic by using formulas of which (20) is 
typical, the Boolean matrix for the system can ie writ- 
ten down. The number of variables involved in the 
system of equations the matrix represents will be equal 
to the number of delay elements in the system. From 
this it can be seen that the Boolean matrix which de- 
scribes the digital system is in general extremely large. 

The matrix representation for the given digital sys- 
tem will then be 


Bu = y (26) 


where each component of the vector x is the present 
state of one of the delay elements in the system. The 
corresponding component of the transformed vector is 
the next state the delay ape will assume. 

The Boolean matrix B will be an array of 1’s and 0’s 
with 2 rows and 2” columns when there are a total of 
n delay elements in the system. 

In the operation of the digital system each clock 
pulse which advances things one step then corresponds 
to one multiplication of the vector x by the Boolean 
matrix B. After two clock pulses have occurred the 
states of all the delay elements in the system are given 


by z where 
B(Bs) = = (27) 
After p clock pulses the states of all the delay elements 
are given by some vector z where 
Bie (28) 


and where the vector « represents the initial states of all 
the delay elements. 

In other words the “power” to which the matrix B 
is raised is determined by the number of clock pulses 
which have elapsed. Suppose that for the system under 
consideration here a total of pulses must elapse before 
the desired process is completed. For example, for addi- 
tion, one word time, in general, must elapse before the 
addition is completed and this might consist of 21 clock 
pulses. 

The transformation which the given digital system 
has performed after # clock pulses is then given by the 
matrix B”. Suppose that it is desired to determine 
whether or not there exists a second digital system 
which will perform the transformation B?” taking a dif- 
ferent period, 
ber of delay elements x as did the original digital system 
with matrix B, and require a minimum of 
diodes. 

Call the matrix of this as yet undetermined digital 
system C. Then in order that the system with matrix C 
perform the desired transformation B? after q¢ clock 
pulses it is necessary that 


say g clock pulses, require the same num- 


number 


(Bei Geer (29) 
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for all values of the initial condition x. This then means 
that element. for element the following matrix equality 


must hold: 
Br = C2, (30) 


If the system under consideration has » delay ele- 
ments then the matrix (30) amounts to a total of 12” 
Boolean equations involving the known elements of 
B and the unknown elements of C. 

The requirement that a minimum number of diodes 
be used in C can be taken into the problem by a second 
set of equations which specify the allowable matrices 
(only those which describe systems with a minimum 
number of diodes). 

The solution to (30) which expresses the unknown ele- 
ments of C in terms of the known elements of B and the 
above equations will then specify the matrix C and will 
represent the desired design with the minimum number 
of logical elements (diodes). 

In considering the problem of optimum design at 
least two parameters are present: the number of de- 
cision and storage elements available. Generally it is 
desired to maintain a certain ratio between the number 
of storage and decision elements available. For example, 
it may be desired to hold the ratio at 25 diodes per flip- 
flop. Eeeee Bert 
An important question which can now be considered 
is: to do a given job what is the minimum number of 
memory cells which are required? In order to answer 
this we must first differentiate between operations in 
the computer which sample the inputs and operations 
which process the inputs. 

Thus a digital computer might accept a set of inputs 
and then spend many clock pulses processing the in- 
puts to obtain an output. The multiplications required 
in the problem are performed sequentially and not at 
the same time. 

All memory elements which are used to hold and se- 
quence the inputs can be eliminated if all processing 
of the inputs occurs at the same time. Such an operation 
would be similar to an analog computer in which all 
operations occur simultaneously. The only memory ab- 
solutely essential is that which is necessary to store in- 
formation from one input sample to the next. Thus, if 
any of the quantities developed depends on the past 
history of any other quantity present in the computer, 
some memory is required. 

The amount of memory required is also a function of 
the amount of information present. Thus, if a particular 
function were always exactly predictable to within the 
last two or three bits, then only these bits would have 
to be st he rest of the quantity could be approxt- 
mated as some function of time. 

In order to determine the least amount of memory 
which will handle a given problem the following pro- 
cedure can be used. Say a given digital system has 
matrix B and that it takes p clock pulses to perform the 
desired transformation B?. The matrix will in general 
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be » by 2” in size-when the total memory capacity of 
the digital system is m bits. What we desire to determine 
is whether a second matrix C can be constructed which 
will satisfy (30). The matrix C will then represent a digi- 
tal system with # storage elements which does the same 
thing as B in q steps. In addition to this also impose the 
restriction on the matrix C that only a certain class of 
systems will be allowed. These allowable systems will 
still have » storage elements, but m of the 2 storage 
elements will 1) not be part of the output storage ele- 
ments (those elements in which the “answer” appears) 
or part of the input elements (where the original data 
is placed) and 2) not be involved in any of the other 
n-m “active” storage elements in the system. 

The problem is that C must have as many total stor- 
age elements as B in order to form the matrix equality 
shown in (30). In order to determine whether a system 
with only m storage elements can also accomplish the 
transformation B’, what is done is to specify that m 
of the storage elements in C will have nothing to do 
with the rest of the system. 

Such storage elements, since they have no part in the 
system, can be removed from the system without af- 
fecting its operation. It is not difficult to specify the 
equations which will only be satisfied by the elements of 
matrices which satisfy the above requirements. For 
example, for »=2, all matrices C for systems in which 
one of the two storage elements has nothing to do with 
the second storage element will have the following form: 


The fact that 


C23 = C21, 
Co4 = C22, 


can be shown to be equivalent to the requirement that 
X, does not appear in the logic for X. This can be shown 
from (11). 

By the use of the above technique it is possible to set 
up a set of simultaneous Boolean equations. The un- 
knowns in the equation will be the elements of the 
matrix C which will satisfy (30) so that it performs the 
same function as does B and also has m “unused” ele- 
ments which can be dropped out. 

Various values of 7 can be tried in this process, start- 
ing from m=n, until‘an m is determined for which a 
solution can be found to the above mentioned system 
of equations. A system with 2-m storage elements will 


_ then be the one which can accomplish the transforma- 


tion B? with the minimum number of storage elements. 
Generally, of course, storage is required in any system 
in which the values any variable will take on as a func- 
tion of time are not completely predictable and where 
it is desired to generate a function which depends on the 
previous values of the variable. When the values are 
completely predictable, “fixed” memory which cannot 
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be altered can be used, as opposed to the use of a 
changeable memory. This type memory can be con- 
tained in the logic for the system. 


Optimum Programming 


The problem of optimum programming arises after a 
general purpose digital computer has been designed and 
built. The programmer is given some specific problem 
which the digital computer is to perform. He then has 
the job of organizing the information in the memory 
and of causing the operations to be sequenced so that 
the required problem can be solved in a minimum time. 
In the sense that the digital computer is already de- 
signed, the programmer’s bounds are fixed. He must 
work within the order code provided with the computer. 
Suppose the programmer has already written a program 
which will accomplish the desired job. This program 
“specifies the initial state of the digital computer. 

Start with all the delay elements in the computer in 
the “0” state and write the program into the computer’s 
memory. This “writing” of the program can be repre- 
sented by the application of a matrix C to a vector all 
of whose components are zero (call this vector x). Then 
the vector which will represent the computer with the 
program written will be y where 


Cx = y. (31) 


The matrix C might be called the “program” matrix 
since it writes the program into a computer. 

Asin the case of the optimum design, reduce the given 
digital computer to an equivalent system which will 
contain only delay elements. Each storage cell on the 
drum, if the computer uses a drum, will then be a single 
delay element. Then construct a single vector « whose 
components are the delay elements comprising the com- 
puter. The logic in the computer acts at each clock pulse 
to transform the vector x into another vector y whose 
components are the respective states of each delay ele- 
ment in the computer at the next clock pulse. The matrix 
which represents this logic can be called B. 

At each clock pulse the logic in the computer trans- 
forms according to the matrix B introduced above. 
After one clock pulse has gone by, for example, the state 
of the computer can be represented by z where 


“BC =e: (32) 
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In other words, the matrix C describes the program 
written into a computer all of whose storage elements 
are initially in the “0” state, and the matrix B repre- 
sents the action the computer has taken on the program 
after one clock pulse. 

Suppose the program takes p clock pulses to perform 
its task. After p clock pulses the state of all the memory 
elements in the computer can be given by some vector z 
where 


BeCe =e. (33) 


It is desired in the design of the optimum program 
to determine whether or not there exists a second 
program which will accomplish the same task but which 
will required less time. Call the program matrix for 
this unknown program D, and call the desired program 
time g. Then if the computer starts with all its memory 
elements in the “0” state, the result of writing the 
program as in (31) will be a vector y where 


Dias: (34) 


After q clock pulses have elapsed, the state of all the 
memory elements in the computer can be represented 
by a vector z, similar to (33) where 


BiDx = 2. (35) 


Now in order that both programs produce identical 
results, it is necessary that the results zg in (33) and (35) 
both be the same. In order that this be true it is neces- 
sary that element for element the two matrices given 
below be equal: 


BoD = Bec, (36) 


As in the case of the optimum design of a digital sys- 
tem, the above matrix equation amounts to a system of 
simultaneous Boolean equations involving the known 
elements of B, C, and the known numbers p and g. The 
unknowns in this system of Boolean equations will be 
the elements of the as yet undetermined program matrix 
BS. 

The solution to the matrix equation, the program 
matrix D, will then represent the desired optimum 
program. If no solution to these equations can be found, 
then no program exists which will perform the desired 
task in qg clock pulses. 
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Summary—When the configuration of a circuit and the equiva- 
lent representations of the transistors are known, a computer pro- 
gram can be written to yield the performance of the circuit and the 
mean values of the circuit parameters. Nonlinearity of the tran- 
sistors is accounted for by piece-wise linearization of an equivalent 
circuit. Rules of interconnection have been devised to combine this 
procedure efficiently with the manipulation of the matrix equations 
that define the linear external circuitry. The genera! method can be 
extended to combinations of basic circuits. 


INTRODUCTION 


N transistor switching circuits, certain combina- 
tions of device characteristics are to be 
most detrimental to proper circuit performance. 
These combinations are relatively hard to find in new 
transistors, and usually occur after the device has been 


known 


in operation for thousands of hours. A method of simu- 
lation has been developed to design transistor switching 
circuits with an over-all reduction in design time. This 
method permits the designer to obseive the operation 
of his circuit with a transistor that has end-of-life char- 
acteristics allowing him to obtain more accurate per- 
formance data. The transient behavior of a transistor 
and its associated circuitry may be analyzed for any 
type of input excitation and under any load conditions. 


THE TRANSISTOR EQUIVALENT CIRCUIT 


The transistor may be classified in the general class 
of four-terminal networks as either a linear active, or 
nonlinear active network, depending on the amplitude 
of excitation applied. In general, input signals of greater 
than five to ten millivolts will put the transistor into the 
nonlinear active class. As the amplitude of the excitation 
is increased to values normally encountered in switching 
circuits, nonlinear circuit-analysis techniques must be 
used to accurately describe its operation. Of the many 
methods advanced to analyze nonlinear systems, it 
appears most feasible, for digital-computer applications, 
to linearize the transistor piece-wise, if a good small-sig- 
nal representation of the transistor can be made avail- 
able. 

Many good small-signal equivalent circuits have ap- 
peared in the literature. Of these, a pi equivalent pre- 
sented in “Transistor Electronics” appears to be most 
adaptable.! 


# Manuscript received by the PGC, June 1, 1957; revised 
manuscript received, September 30, 1957. This paper was presented 
at the Second National Simulation Conf., Houston, Vex., April 12, 
1957. 

+ IBM Corp., Poughkeepsie, N. Y. 

1Lo, Endres, Zawels, Waldauer, and Cheng, “Transistor Elec- 
tronics,” Prentice-Hall, Inc., Englewood Cliffs, N. J., p. 264; 1955. 


The equivalence has been modified from that shown 
in the text, to conform with the requirements needed 
for a switching circuit. These modifications will be 
discussed later. The final form used, which has proven 
very successful, is shown in Fig. 1. The current 14, 


db Yb'b 


Ca BASE b o 


’ | ; 
Ig=oth- wg 1 19 


66 
EMITTER 


Fig. 1—Equivalent circuit for a p-n-p transistor. 


which may be termed a minority carrier current, is 
related to the emitter-base voltage through the diode 
equation 


in = to(edV es’ — J) (1) 
where 


g il 
| KF .0.026yolts 


at room temperature, and 79 is the minority carrier cur- 
rent which flows when the diode is reverse biased. The 
generator which is connected between emitter and base 
is modified by (1—a@o), since the configuration is derived 
for common emitter. However, the circuit is equally 
valid for both common base and common collector. 
The current generator was included to account for the 
emitter-base voltage drop of 0.2 to 0.3 volts when the 
transistor is conducting. Accurate values for this voltage 
are necessary when predicting the response of a network 
driven from a voltage source. 

If the diode equation is differentiated with respect to 
Vs, we obtain a conductance 


d 
G = ——y in = Nice Ve” = Min + ie). (2) 
aV ep 


This conductance is identical to the conductance dis- 
cussed in “Transistor Electronics”? and defines the 
slope of the diode curve at any point. 

The emitter capacitance C, can be considered the 
parallel combination of two condensers. These are the 
barrier capacitance Cyz and the diffusion or storage 
capacitance C,. 


“Toid., p. 241. 
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(Op = Crr + Cr: (3) 


The barrier capacitance Crg*4 results from a contact 
potential formed by the junction of two opposite im- 
purity-type semiconductors, P and N type for the 
emitter and base, respectively. This potential determines 
a depletion or transition region which acts as a condenser 
whose plates are the boundaries of the region. An ex- 
pression which relates the capacitance to the applied 
voltage is 


(4) 


where V2 is the applied emitter-base voltage; Vo the 
contact potential; 2 a constant dependent on the type 
of junction; and K a function of emitter area, dielectric 
constant, conductivity, and width of the depletion 
region. For a drift transistor, typical values® of these 
Benstants are Vo=0;5 volts, -7'=0.4, and K=15) uui- 
(volts)”. 

The storage capacitance acquires a charge equivalent 
to that stored in the base region when the transistor is 
conducting. An expression which relates this capacitance- 
to 1, is 


G A(tn + 40) ; 
. (5) 
2nfs 


where f, is independent of 7.8. 
The common-base cutoff frequency face 1s related to 
two frequencies f, and fr, using the expression® 


Se a 
oes fe a 


where 1/f, and 1/fr are the time constants of the storage 
capacitance, the barrier capacitance, and the emitter 


(6) 


conductance G. 


fans () 
site 2nC; 

2 (8) 

i ie 2rCr ep 

Iq. (6) may be rewritten as 

face ver f (9) 

{eres 

fe 


*L. J. Giacoletto and John O’Connell, “A variable capacitance 
mrmanium junction diode for uhf,” RCA Rev., vol. 17, pp. 221-238; 
h, 1956, 
‘KR. D. Middlebrook, “An Introduction to Junction Tranststor 
rv,” Jolin Wiley & Sons, Inc., New York, N. Y., pp. 58-64, 160-— 
ean Rye [ 
— *D. DeWitt and A. Berger, IBM Product Development Lab., 
Mouchkeepsie, N. Y. (private correspondence); August, 1957. 
“A. L. Kestenbaum and N. H. Ditrick, “Design, construction, 
a high frequency performance of drift transistors,” RCA Rev., 
Vol. 18, pp. 16-19; March, 1957. 
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Since f; is independent of 2, it can be shown, using (9), 
that face approaches f, as 7, is increased. I’or values of 
i, greater than 5 to 6 ma, the following relationship 
holds: 


hace fe. (10) 


This value of f, is used in (5) to determine the storage 
capacitance. 

While f; is independent of 7;, it is dependent on the 
collector voltage. This dependence is introduced by the 
decreasing of the effective base width with increasing 
collector voltage.® 

Using a derivation similar to that for Crg, it can be 
shown that the collector capacitance is related to the 
collector voltage V4 by 


is 
Ch gaVeersa)" 


(11) 


c 


where K, Vo, and 7 are as defined for Crz. 

The time constant 1/wa’ adjusts the phase character- 
istic of the collector generator which would otherwise 
act as a minimum phase network. Measurement of the 
common base short-circuit current gain of the transistor 
shows that the phase angle is actually greater than 45° 
when the magnitude is down 3 db from the low fre- 
quency value. To a good approximation, 

n'G 


C, (12) 


(Qyah = 


where 7’ is a constant for a transistor with a value be- 
tween 2 and 3. 

The conductances g, and g2 are of secondary impor- 
tance in the transient solution. Both values are very 
small and experience has indicated that they can be 
approximated by constant values measured at a cur- 
rent and voltage close to the operating point. 

The conductances gy, and g., are bulk conductivities 
in the base and collector, respectively. The base bulk 
conductance is very important to the transient response 
while g.-, appears to have little or no effect when the 
transistor is kept out of saturation. The value of gy-5 is 
obtained from a measurement of common-emitter input 
impedance with the collector short circuit. The measure- 
ment is made at 30 mc. Plots have been made as a func- 
tion of coliector current and voltages and the conduct- 
ance has been found to be relatively constant. 

The two remaining nonlinear parameters of the tran- 
sistor are the current gain and the cutoff frequency. Both 
of these can be measured asa function of emitter current 
and plotted to determine a functional relationship for 
them. 

The system of differential equations which describe 
the equivalent circuit may be written in nodal form in 
terms of the difference of voltages. 


7 Lo, ef al., op. cit., pp. 251-253. 
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The value of 7, 7, and 7, may be obtained by first inte- 
grating the differential equations for the differences of 
the internal node voltages. The resulting difference 
voltages and (16) through (20) are used to obtain %, 7,, 
and 2. 


Veen Vig Vm Ve) (16) 
Ve = (Ve. — Vor) + Ve (17) 
ip = (Vir — Via)goo (18) 
hem Ve =" Ve) ake (19) 
de = 15 + 4. (20) 


Two requirements must be satisfied before attempt- 
ing a successful piece-wise linearization of a nonlinear 
system: the coefficients of the equation as a function of 
the operating point must be known; and sufficiently 
small intervals of integration must be maintained so 
that the linear assumption is not violated. The func- 
tional relationships for the coefficients can be deter- 
mined from direct measurement and/or from the physi- 
cal considerations discussed previously. A numerical 
method of integration is ideal for limiting the inte- 
gration interval if the maximum value is directly con- 
trolled as input information. 


Tuner NEtworK ANALYSIS PROGRAM 


A computer program which will meet this require- 
ment is the SHARE program PK-NIDA (a numerically 
integrating differential analyzer). The mathematical 
basis for solution of the differential equations is Adam's 
method of numerical integration. Several logical fea- 
tures of the routine control the accuracy and speed 
of solution. The maximum integration interval and 
the allowable error are under control of the user. In- 
tegration may proceed at an interval which ranges be- 
tween the assigned maximum interval and 271! of this 
value. The allowed error specifies the accuracy required 
between successive predicted and corrected values of 


For example, if an allowed 


the variables. error of 2710 
is used, accuracy to ten binary places is required. If pre- 
dicted and corrected values disagree within these ten 
binary places, the program halves the integration inter- 
val, reverts back to the last good values, and proceeds 


with this smaller interval. If predicted and corrected 
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values agree to 16 or more binary places, the program 
will double the integration interval and the integration 
continues at this larger interval. This halving or dou- 
bling of the interval can occur at any time during the in- 
tegration process, except when first starting or when se- 
vere discontinuities are encountered. For these cases a 
set of starting equations are employed which use an 
interval 2~'!6 of the maximum. The values of the vari- 
ables and their derivatives are stored in tables which 
contain the last eight computed values. These tables 
are continually updated as new values are determined, 
with the latest values placed at the top of the tables. 

In the investigation of methods of evaluating actual 
circuit performance, it would appear on first sight that 
the method for analyzing the nonlinear transistor is 
incompatible with the analysis of the linear external 
circuitry. Manipulation of matrix equations is ideal for 
linear systems, while a method such as the one just de- 
scribed is best for a nonlinear system. However, an effi- 
cient combination of these methods has been achieved 
by a program which separates the transistors from the 
remainder of the circuit and produces a solution based 
on the most efficient analysis for each part. 

The linear circuit elements are converted to matrix 
form in accordance with a set of simple rules which may 
be conveniently programmed.® Fig. 2 shows a simple RC 
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Fig. 2—A simple RC network and generator. 


and generator used to demonstrate this 
method. The differential equations for the system are 
shown. If the value of the generator is known, the net- 
work may be solved in matrix form by 


network 


CXEC GX 0 
where 


C is the capacitance matrix shown; 

G is a combination of a conductance matrix and a 
unit matrix; 

X isarate vector; 

X isacurrent and voltage vector. 


*J. C. Morgaun, Share Routine PE-CAML, IBM Product Dev. 
Lab., Poughkeepsie, N. Y. 
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If the C matrix is inverted and a matrix multiplication 
performed, the final expression X=AX is obtained, 
where .1 is the result of premultiplication of the —G 
matrix by the C7! matrix. For a machine to perform 
these operations, it must be given a precise description 
of the circuit as input information. This description 
must include the following: 


1) The number of independent nodes (so that the 
size of the C matrix can be determined). 

2) The number of inductive currents. 

3) All branch elements and the direction of positive 
current flow through them. 


It is worth noting that the C matrix cannot be singular 
since it must be inverted. The singularity indicates that 
there is a dependent node associated with this matrix, 
and this dependency can be eliminated only if some 
stray capacitance is inserted at the node. 

The forming of this matrix equation is performed en- 
tirely by the machine, based on the input data. Infor- 
mation is punched into cards, each card containing the 
type of element, its magnitude, the direction of positive 
current flow, and the nodes to which it is connected. 
With these data, the machine forms the C matrix, —G 
matrix, and the 4 matrix. Since the final 4 matrix con- 
tains only elements of the linear system, it is unneces- 
sary to recompute it unless some component value is 
changed. As an example, we can consider how the ma- 
chine forms the C matrix for_this network. Since there 
are two independent nodes, this information is used to 
determine a 2 by 2 matrix. The condensers are located in 
the matrix by indicating the nodes to which they are 
connected, and enumerating these nodes in the order 
that determines positive current flow. For this network, 
the information is: C; between nodes 1 and 0, and C2 be- 
tween nodes 2 and 0. These condensers are then inserted 
with a negative sign in locations 1, 1 and 2, 2 of the 
matrix. If a condenser had been located between nodes 
1 and 2, with positive current in this direction, it would 
have been inserted negatively in locations 1, 1 and 2, 2 
and positively in locations 1, 2 and 2, 1. 

If a transistor switching circuit like the one shown in 
Vig. 3 is considered, the circuit may be divided into two 
parts. One part contains the external linear circuitry, 
and the other the nonlinear transistor. Since it is known 
that matrix equations afford the most efficient method 
for handling linear systems, a look at the general matrix 
equation will reveal a method of linking the transistor 
.to this solution to obtain a solution for the entire system. 

Fig. 4 shows this equation along with the modified 
Circuit which will be considered. The matrix equation is 
X=: .1X, where the rate vector X contains all the de- 
iivatives of the independent node voltages and currents. 
the vector X contains the magnitude of these same 
currents and voltages, as well as the fixed bias voltages 
aid the voltages or currents supplied by generators. 
The A matrix contains the values of all the passive cle- 
Mcnts of the linear circuit. Integration of the rate vector 
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circuit. 
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X= AX 


Fig. 4—Modified configuration of the switching circuit. 


over a small interval of time results in new values for the 
X vector. If the integration interval is controlled so that 
the X vector, containing the transistor node voltages, 
does not change enough to violate the linear treatment 
of the transistor, it is possible to justify the insertion of 
a current generator at each node to which the transistor 
is connected. These current generators are shown in the 
circuit diagram. The X vector will now contain values 
for the node voltages and currents and the external 
voltage generators, plus the transistor current gener- 
ators. 

Fig. 5 shows the actual matrix which would be con- 
structed by the machine for this circuit. It can be seen 
that the equations for the inductive currents have been 
added to the A matrix after the multiplication of the 
—G matrix by the C-! matrix. 

The interconnection between the solutions of the 
linear and nonlinear parts of the circuit is effected 
through the X vector. Initially, the transistor current 
generators are held constant and the rate vector is 
evaluated. This rate vector is then integrated over a 
controlled interval of time to obtain new values for the 
X vector. The voltages in the X vector which are asso- 
ciated with the nodes to which the transistors are con- 
nected are transferred to the routine which solves the 
transistor equations. Using these voltages, the differen- 
tial equations are then evaluated and integrated over 
the same time interval to produce new values for the 
current generators. These new values are inserted back 
into the X vector, and a transfer is made to the matrix 
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Fig. 5—-Constructed matrix equation. 


routine where a new rate vector is computed. This iter- 
ative process is continued until no further change is ob- 
tained at the output for some input excitation. 

Since minimum storage requirements and program- 
ming ease are of primary importance in this simulation 
method, these two points have been considered in de- 
veloping the program. The first has been accomplished 
by accommodating all transistors with one NIDA 
routine. The second has been achieved by locating 
automatically all variables pertinent to a transistor, 
given only the number of the transistor and the location 
of its node voltages and current generators in the X 
vector. : 

To initialize the system, the above information must 
be given for each transistor. In addition, the user must 
supply the following: data needed to create the A 
matrix, the maximum integrating interval, the allowed 
error, the number of transistors in the circuit, and 
parameter values for the transistors. 

With the matrix equation written and the transistor 
subroutine initialized, the program transfers to the 
linear circuit routine and computes the values of the rate 
vector from the matrix equation. It then proceeds to 
predict and correct new values for the X vector over a 
fixed integration interval equal to the maximum value 
inserted into the program. These new values are then 
stored back into the X vector and the program transfers 
to the transistor. 

This routine handles the transistors one at a time. 
The values ofthe base, emitter, and collector voltages 
are obtained from the X vector and are used first to 
modify any parameters which are a function of them. 
Other paramet 
ter Current. 

Since only one NIDA routine is used for all tran- 
sistors, it isimportant to retain the integration interval 
that was used on the transistor in the previous pass. This 
is necessary since the eight values stored in the table 


ers are modified as a function of the emit- 


for each variable are spaced in time according to the 
interval. The importance of this becomes apparent if 
we consider the operation of two successive transistors 
in the routine. If the first of the two transistors were 
operating in a region where no discontinuities were pres- 
ent, NIDA would be using an integration interval near 
or equal to the maximum value. If the second transistor 
had in the meantime encountered some difficulty, it 
would be using much smaller integration intervals. As- 
suming that the routine had just completed the first 
transistor and passed on to the second, NIDA would 
still have a large interval stored. It would interpret the 
tables as being spaced in time according to this large 
interval instead of the small interval for which they 
were computed. Using the wrong interval would give 
an erroneous value for the next point. With the proper 
interval restored into NIDA, by address modification, 
the differential equations are evaluated and integrated, 
and the tables updated. The program will continue with 
this transistor until it has progressed exactly one maxi- 
mum integration interval. It is important to note here 
that this interval must be maintained equal for all tran- 
sistors in order to maintain a time relationship between 
the linear and nonlinear part. When the computation 
of the last transistor is completed, the values of all 
emitter, collector, and base currents are stored back into 
the X vector, and the program returns to the matrix 
equations. 

Figs. 6 and 7 show the computed response of the cir- 
cuit in Fig. 3. A voltage source with a fall time as indi- 
cated was used to drive the circuit in both cases. Iden- 
tical transistors with the indicated characteristics were 
used to compare the effect of a change in capacitance on 
the collectors of both transistors. As would be expected 
the addition of the capacitance reduced the overshoot 
and increased the rise time, but did not affect the turn- 
on delay caused by having the transistor heavily cut off. 

The program thus far has attempted to show the feas- 
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. 6—Unloaded response of the circuit in Fig. 3. 


ibility of such an approach in determining the response 
of transistor switching circuits. It should not be inferred 
that the method is a substitute for the preliminary de- 
sign of circuits, or a means of gaining insight into the 
operation of a circuit. It should be considered instead as 
a substitute for the final bench setup used to obtain 
data. The program needs further work before its use 
can be widely advocated. At present the program 
requires 15 to 20 minutes to obtain a complete response, 
and its capacity is a circuit containing up te ten tran- 
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Fig. 7—Loaded response of the circuit in Fig. 3. 


sistors. The next phase of this program will concentrate 
on methods of reducing the computation time and 
making the routine more flexible. 
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An Optimum Character Recognition System 
Using Decision Functions’ 
C. K. CHOWt 


Summary—The character recognition problem, usually resulting 
from characters being corrupted by printing deterioration and/or 
inherent noise of the devices, is considered from the viewpoint of 
Statistical decision theory. The optimization consists of minimizing 
the expected risk for a weight function which is preassigned to meas- 
ure the consequences of system decisions. As an alternative, mini- 
mization of the error rate for a given rejection rate is used as the cri- 
terion, The optimum recognition is thus obtained. 

The optimum system consists of a conditional-probability densi- 
Sities computer; character channels, one for each character; a re- 
jection channel; and a comparison network. Its precise structure and 
ind ultimate performance depend essentially upon the signals and 

‘oe structure. 

Explicit examples for an additive Gaussian noise and a ‘‘cosine’’ 
“9156 are presented. Finally, an error-free recognition system and a 
7 >sible criterion to measure the character style and deterioration 
Se presented. 


; Manuscript received by the PGEC, June 3, 1957. 
Burroughs, Corp., Paoli, Pa. 


INTRODUCTION 


HARACTER recognition has been receiving con- 

siderable attention as the result of the phenomenal 

growth of office automation and the need for 
translating human language into machine language.!? 
Broadly speaking, the character printed in conventional 
form and size on the document (checks, etc.) is first con- 
verted to electrical signals, and sufficient information is 
then extracted from the latter. The purpose of the 
recognition system is based on the observed data and 
on a priort knowledge of the signal and noise structure 


1k. R. Eldredge, F. J. Kamphoefner, and P. H. Wendt, “Auto- 
matic input for business data processing system,” Proc. Eastern Joint 
Computer Conf., pp. 69-73; December 11, 1956. 

2 Ie. C. Greanias and Y. M. Hill, “Considerations in the design of 
character recognition devices,” 1957 IRE Nationa CONVENTION 
Rercorp, pt 4, vol. 5 pp. 119-126. 
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to identify which of the possible characters is present, 
or to reject if the data are ambiguous. 

The over-all performance of the recognition system 
depends not only upon itself, but also upon the number 
of characters to be recognized, the character style, and 
noise statistics. In this paper the character style and 
noise statistics are assumed given and adequate, and 
the purpose of the paper is to obtain an optimum recog- 
nition system. For convenience, the recognition prob- 
lem is considered one of statistical inference, so that use- 
ful results in decision theory can be applied.*> To ac- 
complish this, the notion of risk is employed and proper 
weights are assigned to various types of error, rejection, 
and correct recognition to measure the consequences of 
decisions. This results in an optimum system which 
minimizes the expected (average) risk function and in- 
cludes a possible alternate system with a minimum 
error rate. The results reveal the explicit structure of an 
optimum system which is determined by the a priori 
noise statistics, the signal structure, and the preassigned 
weights. 


SYSTEM APPROACH TO THE PROBLEM 
One practical application of a character recognition 
system for business documents is to read arabic numer- 
als and selected symbols printed in magnetic ink. A 
method? for achieving this is shown in Wig. 1. The char- 
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Fig. 1—A recognition system. 


acter is first passed through the field of a permanent 
magnet where it is magnetized in a given direction be- 
fore being scanned by the read head. From the read 
head, the printed character is converted into an elec- 
trical signal corresponding to the differentiation of the 
plane area of the character. The function of the recogni- 
tion system is to examine the amplitude-time signal 


SA. Wald, “Statistical Decision Functions,” John Wiley & Sous, 
Inc., New York, N. Y., 1950. 

4D. Van Meter and D. Middleton, “Modern statistical ap- 
proaches to reception in communication theory,” IRE TRANS., 
vol. PGIT-4, pp. 119-145; September, 1954. 

§D. Middleton and D. Van Meter, “On optimum multiple-al- 
ternative detection of signals in noise,” IRE TRaAns., vol. IT31, pp. 
1-9; September, 1955, 
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obtained by the read head and to decide which of the 
possible characters is being recognized. 

It is convenient, at times, to deal with the sampled 
data rather than the continuous time waveforms. By the 
sampling theory, if the number of samples is sufficiently 
large, little information carried by the continuous signal 
is lost. As shown in Fig. 1, the signal from the read 
head is first amplified and then fed into a tapped delay 
line. This serves as a means for sampling and acts asa 
temporary storage device to convert the series informa- 
tion into parallel information. Although not essential, 
sampled data are used in the following discussion. 

Let the vector v= (1%, 2», - vs) (subscript s being the 
number of samples) denote voltages on the taps of the 
delay line at the instant of sampling. (See Appendix I 
for the meaning of the Symbols.) The vector a;= 
(du, G2, + + +, Ais) denotes the true sampled signal asso- 
ciated with the 7th character where i=1, 2,-::,6,€ 
being the number of possible characters to be recog- 
nized. The vector v constitutes the input to the recog- 
nition system. It is assumed that the characters are dis- 
tinct, 7.e., all a;’s are different. 

Ina simple form, the recognition systeni may consist 
of ¢ separate channels, one for each character. Each 
channel obtains a weighted sum of »,’s, with properly 
chosen weights, 0;;. The output of the ith channel is 


s 
X ;(v) CoS ae bi ;0;. 
j=1 (1) 
This operation may be realized by a summing amplifier 
and possibly with some inverters to provide negative 
weights, if required. One possible set of weights is: 


(2) 


The recognition system is known as a correlation net- 
work when the weights are defined by (2). 

If the printing is perfect, and the reading devices are 
noiseless, the observed data v will be identical to one of 
the a;’s and therefore, it can easily be shown that the 
right channel of the correlation network has the largest 
(algebraic) output. Consequently, the recognition sys- 
tem identifies the character with absolute accuracy by 
the channel having the highest output. However, in 
practice, there are always, to some degree, deteriora- 
tions in printing and inherent noise in the devices. 
Therefore, the observed data v generally will not be 
identical to any of the a,’s. In view of this, ambiguities 
arise which may result in possible misrecognition. To 
safeguard against the occurrence of error, the recogni- 
tion system should have provisions for examining the 
degree of ambiguity and making rejects when required, 
This function can be achieved in various Ways; €.¢., 
whenever the next highest output of the correlation 
network exceeds some preassigned fraction of the high- 
est output, the system will reject, otherwise the system 
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identifies the character by the channel having the high- 
est output. 

The system described above merely represents one of 
many possible recognition systems and is not necessar- 
ily optimum. A basic problem in the design of recogni- 
tion systems is to evaluate the system performance in 
the presence of priuting deterioration and inherent 
noise and to obtain an optimum system. Optimum per- 
formance depends primarily upon the character style 
and permissible deterioration. Greanias and Hill in a 
recent paper’ describe the effects of character style and 
printing deterioration on the character recognition 
problem from the viewpoint of matching the character 
with an ideal character and further propose definitions 
for character quality and style factors. In this paper, 
the discussion is confined to the problem of obtaining 
an optimum recognition system for a given set of ade- 
quately styled characters and known statistics of char- 
acter deterioration. The recognition problem is consid- 
ered to be that of testing multiple hypotheses in the 
statistical inference. Consequently, the design and 
evaluation of a recognition system is comparable to a 
statistical test. Results of decision theory can be ap- 
plied.” >. 

In order to judge the relative :nerit of recognition 
systems, some criterion of evaluating system perform- 
ance must be established. The error rate of the system 
for a given rejection rate is used as the performance 
criterion for cases where no distinction is made among 
-‘misrecognitions. Cases may arise where different mis- 
recognitions have different consequences; e.g., the regis- 
tering of a four as an asterisk may not be as serious an 
error as registering it as a nine. The criterion of mini- 
mum error rate is then no longer appropriate. Instead, 
the criterion of minimum risk? is employed. Proper 
weights are assigned to measure the consequences of 
errors, rejections, and correct recognitions. These 
Weights indicate the loss incurred by the system for 
every possible decision. The loss, which should be re- 
garded as negative utility, may actually represent loss 
in dollars or unit of utility in measuring the conse- 
quence of the system decision. The over-all performance 
of the system is judged by its expected (or average) risk. 

In the following discussion, an optimum system 
which minimizes the expected risk is derived, and a sys- 
tem having minimum error rate is obtained. Examples 
are presented for illustration purposes. An error-free 
system and a possible criterion for judging character 
stvle and deterioration are also presented. 


THe EXPECTED Risk 


The vector a; (aa, @x, *, Gis) n the s-dimensional 
space denotes the true sampled signal associated with 
the ith character (¢=1, 2, - > 
respectively, the number of possible characters to be 
recognized and the number of samples. Let p= (fi, po, 
"++, pb.) be the a priori distribution of characters (pi 
is the @ priori probability that the ith character occurs). 


-, c), where ¢ and s are 
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Then, evidently, 


c 
Dodi lee, pee: (3) 
21 
The received data are denoted by a s-components vector 
v=(v1, v2, +++, Us). It is the signal corrupted by factors 
such as the deterioration of printing and inherent noise 
of the devices. A priori noise statistics and the manner 
in which various signals and noise are combined deter- 
mine precisely the conditional probability density 
F(v}a;) of the observed data v when a; is the incoming 
signal. 

The space of decisions available to the recognition 
system consists of ¢+1 possible decisions do, di, da, 

-, d.. The quantity d;(70) is the decision that the 

ith character is present while do is the decision for re- 
ject. A basic problem in statistical decision theory is the 
selection of a proper decision rule 6. The rule is expressed 
as a vector function of the data v, namely, 6(v) = (5(do| v), 
5(d,| 2), d(d2|v) + + - 6(d.]v)) with ¢+1 components, and 
satisfies the restriction that: 


MS. 5(d; | v=o for all 2, (4) 
i=0 
and 
6(d; | Oy teal for all 7 and »v. (5) 


The quantity 6(d;| v) is the probability that, for a given 
observed data v, the decision d; will be made. 

In order to judge the relative merits of the decision 
rules it is necessary to assign the weight function 
W(ai, d;). This is a function of a; and d;, which is the 
loss incurred by the system if the decision d; is made 
when a; is the true signal. This measure of consequence 
for various d; under various a; is a datum of the problem 
and is given in advance. Let the weight function be: 

W(a:, dj) = wi t meee a (6) 
iy Oj a Patino te O 
where w;i(140) is the weight of correct recognition. of 
the 7th character; wi;(t4+j7 #0) is the weight of misread- 
ing the zth character as the jth one, and wio(7#0) is the 
weight of rejecting the 7th characters. Therefore, it is 
required that 


(i #7 #0). (7) 


Wij > Wg) = aes 


Usually, w;; is much larger than wy since the most seri- 
ous consideration in design of a character recognition 
system is the occurrence of undetected errors. 

In general, w,;’s may all differ, so that various mis- 
recognitions, rejections, and correct recognitions can be 
properly weighted. The expected risk for any decision 
rule 6 is 


R(p, 8) = Do Dd) | 8(d;| 0) piw.jF(v] asda, (8) 


i=1 j=0Y V 


with integration over the entire observation space V. 


THe Minimum Risk System 


The problem is then to choose a decision rule to 
minimize the average risk. By using (4), and since 


IvF(v| a i)dv =1 for all 7, (8) may be written as:] 


R(p, 8) = Rol) + RilP, 4), (9) 
where 

Ro = >> pwio, (10) 

ical 
Ri = | do d(d;| 2)Z,(v)dz, (11) 
V j=0 
and 
iv (wi; — wo) piF (0 | a3 7 = 1, 2 é 

Z;(v) = t=1 , (12) 


0 at OL tree Ds 


The symbol Rp will express the expected risk when re- 
jection is made for all recognition and R, is that part of 
R which may be adjusted through the choice of 6. Evi- 
dently 


Ri(P, 8) ak min [Z,(v) ]dv,, (13) 


ven ee w 


and the equality sign holds if, and only if, the decision 
rule is chosen as: ; 


5*(d,| v) = 4, 


8*(d;| v) = 0 (14) 


whenever 


min [Z;(v)] = Z,(2). (15) 
z] 

This is the optimum decision rule 6* (the Bayes strat- 
egy) which minimizes the expected risk and is non- 
randomized since its components are either zero or one. 
Therefore, R, for this decision rule is always nonposi- 
tive, and its expected risk (the Bayes risk) is no greater 
than Ro. The expected risk for the optimum decision 
rule, 6*, is 


R(p, 8*) = DD prwin +f min [Z;(v) ]do. (16) 
=] ‘ j 


J 


Eqs. (14) and (15) reveal that the optimum system 
for character recognition consists of a computer which 
evaluates F(e|a,)’s; (¢=1, 2,---,c) for an observed 


data v; computes the various Z,(v) (7=1, 2,---, c); 
examines and compares these Z;(v) (7=0,1,2,---, G): 


selects the smallest (algebraically) one, say Zi (v); and 
finally makes the decision d, [having the same subscript 
as Z,(v)]. Of course, this method of setting up the com- 
puting procedures is not unique; e.g., any ordering-pre- 
serving transformation may be used. In any event, the 
system must be equivalent to the above. 
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OBSERVED . 
DATA 


STATION NO.I 


F(vlag) 


STATION NO.2 


COMPARISON STATION (MINIMUM DETECTION) 


| cca el eats 


do | dx dc 


Fig. 2—Functional diagram of an optimum system. 


The functional diagram of the optimum system is 
shown in Fig. 2. Station No. 1 consists of c similar com- 
ponent networks. Each network receives the observed 
data v= (u,v, - - + , vs) and computes the corresponding 
conditional probability density F(v|a;) as its output. 
This operation depends only upon the a priori knowl- 
edge of signal and noise structure and on the observed 
data v; it does not depend upon the weight function 
W(ai, dj) or on a priori probability distribution of sig- 
nals, p. 

The outputs of station No. 1 are fed to station No. ys. 
which consists of ¢ character channels N2;(j70), and 
one rejection channel, No. They perform the linear 
operation of weighting each input and then the sum- 
ming of all weighted inputs. The weights are Diss 
and p,wio's. The output of the rejection channel is 


Xo(v) = Ds wiopik (0 | a:), (17) 


while the outputs of character channels are 
X,(v) = 2) wip FO| a) (G=1,2++-6). (18) 
t=} 


The comparison station receives X’s from station 
No. 2, examines all its inputs, and makes decision by 
selecting the algebraically smallest of the c-++1 X's. Tf 
the rejection channel has the smallest one, the system 
rejects. If one of the character channels has the smallest 
output (say X,(v), (k¥0)) then the system recognizes 
the signal as the Ath character. 

Since X;(v) ({=0,1, +++, ¢) is equal to Zj(v) +X (v),. 
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this system makes decisions in accordance with 6* as 
defined by (14) and (15), and thus minimizing the ex- 
pected risk. 


PROBABILITIES OF ERROR AND REJECTION 


The expected risk provides a means for evaluating 
the performance of a recognition system. At times, it 
may be desirable to compute the probabilities of error, 
rejection, and correct recognition as an auxiliary set of 
“merit figures. They are obtained for any decision rule 6 
as follows: 

The probability of correct recognition is: 


P.8) = | 27 d(di| 2) pF(o| a:)do; (19) 
Ve onl 
the probability of rejection, or rejection rate, is: 
P,(6) = f acao| 0) >» pF (o| a;)dv; (20) 
Vv val 


and the probability of misrecognition, or error rate, is: 
PUG el ep po (21) 


Eqs. (19)-(21) result directly from the fact that 
f5(do|v) F@v|a,)dv and fr5(d,|v) F(vlaidv are, respec- 
tively, the conditional probabilities “f rejection of the 


ith character and correct recognition of the 7th character. 


ry 


i 


CRITERION OF MINIMUM ERROR RATE 


Cases may arise in which the criterion of judging the 
system performance is the magnitude of its error rate 
for a given rejection rate. In using this criterion, the 
optimum recognition system is the one which, for a 
given rejection rate, a, has a minimum error rate. The 
optimum decision rule is obtained as: (See Appendix II 
for proof.) 


o**(d, |v) = 1 (k ¥ 0) (22) 
whenever 
PF (v| a) = PF (0 | a;) for ally ~ k, and 
piF(v| ax) > BD) pF (o| ad), (23) 
feat 
and 
5**(do| v) = os (24) 


Whenever 
BY) pF(v| a) > piF(o| a;) forallj(j = 1,2 +--+ 0), (25) 
i=] 


Where B(0 <6 <1) is a nonnegative constant determined 
by the condition that P,(6**) =a; namely, 


4k F(v| a;)dv = a. 
fi (dol ») Dp (| ai) (26) 


rhe constant B increases with increasing a, and P.(6**) 


and P.(6**) are monotonic decreasing functions of a. 
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etre 4; de do 


Fig. 3—Functional diagram of the system having the minimum 
error rate. 


The change in the constant 6 provides a control over 
the error-reject ratio. 

The optimum rule provides the basis for the func- 
tional diagram of the system of minimum error rate as 
shown in Fig. 3. The first station is identical to that of 
the minimum risk system (see Fig. 2) which computes 
the conditional probability densities F(o| ai)’s. 

The second station for this system is somewhat simi- 
lar to that shown in the functional diagram for the 
minimum risk system (see Fig. 2). The c-character 
channels perform the weighting operation and have 
biF(v| as) as outputs. The rejection channel, No, per- 
forms the operation of weighting and summing and has 
paver biF(v| a;) as its output. All of these c+1 outputs 
are nonnegative. The comparison station then examines 
these outputs and selects the largest. If the output of 
the rejection channel is the largest, the system rejects; 
otherwise the system will identify the character by the 
channel having the largest output. 

It can be shown that the system depicted in Fig. 2 re- 
duces to the system shown in Fig. 3 when 8 is replaced 
by (Wa—Wr)/(Wm—wW-), and the following weight func- 
tion is used. 


or* 


WwW, 10% 1-=7, > 0 
W(a:,d)) ={w, fori#0,7 = 0. 

Wm, fOr Paq 52, 

EXAMPLES 


1) As an illustration, consider a condition where the 
signals and noise are additive, and the noise has inde- 
pendent normal distribution. To be explicit, the prob- 
ability density function of the noise of the jth compo- 
nent of the 7th character is taken as: 


292 IRE TRANSACTIONS ON 
I ( (v — a,;)*) 

Ld tj/ , 

oa exp 4 — (28) 
\/ 270 ij 20%; 

where v;—a,; is the noise and o;;? is the given variance. 


The conditional probability density Flay, under 
the assumption that noise is statistically independent, is 


L danaes Orisa aay 
ie siete 


20537 


se a 
km) i Il Ci; 


j=1 


(29) 


Therefore, the last expression dictates the precise struc- 
ture of station No. 1 for the optimum system. Each com- 
ponent circuit performs the operations of taking differ- 
ences, squaring, and taking ex- 
ponential. This station is common for a minimum risk 
or minimum error rate system. The structures of the sec- 
ond station and comparison station are indicated in 
Figs. 2 and 3. 

2) In this example, the signal and noise structure are 
such that the conditional probability density for a 
given length, | v| , of vis directly proportional to the co- 
sine of the angle @ between vectors v and a; for |6| <1/2 
and is zero elsewhere, and that the distributions of | o| 
for given a; are identical for all 7. [It is denoted as 
f(f2| ).] In other words, Fla :) can be written as 


weighting, summing, 


F(o| a;) = for a;:v > 0 


elsewhere, 


(30) 


where p is a constant independent of 7 and is determined 
by the fact that fv Fv a;)dv=1, and a;-v denotes the 
scalar product of vectors a; and v. 

An inspection of the optimum decision rule (6* or 
6**) reveals that the system remains optimum if the 
first station is to compute T[F(v| a;)] instead of F(v| ai), 
where 7 is defined as 


: v 
T(v| a;) = ——— F(v] a; 
(2 | a;) of(le) (2 | a 
ay £ 
—-y bn; f ct > 0 
ba Geshe IRS BET Siac loess 
0 for a;-v < 0, 


where ,;’s are constants [see (2) ]. This operation can be 
easily realized. Each component of the first station is 
simply a correlation network followed by a half-wave 
rectifying circuit. The circuit passes the positive output 
of the correlation network unaltered and converts its 
negative output into zero. 

The e that the recognition 
system described in the second section of this paper is 


above results also indicat 


not optimum for the particular signal and noise struc- 
ture as given in examples 1 or 2. 


§ This particular signal and noise structure was anes eS the 
author’s colleague, I. M. Sheaffer, Jr., Burroughs Corp., Paoli, Pa. 


ELECTRONIC COMPUTERS December 


ERROR-FREE SYSTEM 


For convenience, let V; denote the set (or region) of 
all possible observation v when the 7th character is pres- 
ent, and let V; be the largest subset of V; so that V;’s 
are nonoverlapping. If noise distributions are so trun- 
cated and the signal vectors a,’s are so placed that all 
Vs are nonempty, then an error-free system for char- 
acter recognition does exist. eas for any observed 
data v belonging to V’,, only F(@ izero while all 
the character can then be identified 
with certainty. On the other hand, if the data v do not 
belong to any one of the V,’s, then more than one of the 
F(o| a;)’s will be nonzero. This results in the data being 
ambiguous for recognition purpose, and an error-free 
system will reject. Symbolically, the error-free decision 
rule is: 


others are zero; 


S(d,|v) = 1 if F(v| a) > 0 


and 


F(vla) =0 foralli#¥ k, (32) 


and 


5(do| v) = 1, otherwise, 


the rejection rate is determined by the probability 
measures of V’;’s, namely /7,F(v|a,)dv. The latter is de- 
termined by the character style and allowable deteriora- 
tion. The character style may be considered ideal and 
the control over the printing perfect, if the resultant 
S¥, i)dv is unity for all z, and all characters with al- 
lowable deterioration can then be recognized with 
neither an error nor reject. In this sense, the probability 
measures of V’;’s may be used to evaluate the combined 
quality of the character style and printing. 


U 


CONCLUSION 


The decision theory has been successfully applied 
to the problem of character recognition. By employ- 
ing the concept of risk, differences in consequences 
for various decisions have been taken into considera- 
tion. A rejection channel has been introduced to exam- 
ine the degree of ambiguity of input signal and make 
rejections when necessary. 

As developed, the structure and performance of an 
optimum system depend upon the signal and noise sta- 
tistics; therefore, a priori knowledge of these statistics is 
required. Usually, a realistic estimate of noise statistics 
is not easy to obtain. However, it is sincerely felt that 
the requirement for high grade performance in character 
recognition warrants the expenditures in this direction. 

Quite often an optimum system may prove to be too 
expensive for mechanization. Nevertheless, the results 
presented in this paper are considered useful in that they 
provide insight into the recognition problem and furnish 
an ideal system, which actual recognition circuitry may 
be patterned after. 

Although it is recognized as being beyond the scope 
of this paper, it is worth mentioning that one practical 
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ipproach to the over-all problem would be to design 
idequately the character style and to control properly 
she printing process so that a reliable system would not 
ye too far fetched or difficult to ultimately realize. 


APPENDIX I. 
List oF SYMBOLS 


i= (Ga, G2, ++ +, Gis), S-dimensional vector as- 
sociated with the ith character, (¢=1, 2, 
vec): 
ai;=jth sample of the signal of the 7th character. 
| a; zee a;;?)1/?, the length of vector aj. 
¢=nurnber of characters. 
0 = decision that rejection be made. 
d;=decision that the signal is the jth character 
Ca ee naire) 
fd o| ) =probability density of the length of v. 
F(v| a.) = conditional probability density for the ob- 
served data v when a; is the incoming signal. 
1,j, k=indexes. 
Ni;=a network of station No. 1,24=1, 2,:--, 6. 
Near Net wormor station NO: 2, t=1.2, °° =, 6 
No=the network of rejection channel. 
P.=probability of correct recognition. 
P,=probability of rejection (reiection rate). 
P,.=probability of misrecognition (error rate). 
bi=ea priori probability that the zth character 
occurs, (= 1,262" -, ¢). 
b= (h1, ps, gate pee 
R(p, 6) =expected risk of the system; R=Ro+ Ri. 
Ro(p) =expected risk of the system when rejection 
is made for all recognition. 
Ri(p, 6) =part of R which is dependent upon 6. 
s=number of samples. 
T =functional transformation. 
V=s-dimension observation space. 
V;=set of all v when the 7th character is present. 
V;=largest subset of V; such that ViV;=0 
for all 7 47. 
v=(v1,¥2, -- + ,Us),a vector in V. 
|| = (04.1 2,7) "?, the length of vector v. 
vt =ith component of the observed data v. 
W(a:, dj) =weight function. 
Wij= W(ai, d;). 
Ye, Wp, Wm = Weights. 
x,(v) =output of the zth channel. 


a =permissible rejection rate. 
B=constant. 
8(v) =decisionrule,6 =[6(do| v), 6(dilv) - - - 6(d.|v)]. 
6*(v) =optimum decision rule which minimizes the 
expected risk. 


6**(y) =optimum decision rule which minimizes the 
CGO tau. 
6=angle. 
+  p=a normalizing constant. 
¢;;?=statistical variance of noise. 
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To Prove Tuat 6** Has a Minriuum ERROR 
RATE FOR A GIVEN REJECTION RATE 


Without loss of generality, it is assumed that the ab- 
solute probability density of the occurrence of v, namely, 
Dat ee b:F(v| a) is nonzero over the entire observation 
space V. Otherwise, the set over which )_7, biF(v| a;) 
is zero is first deleted. 

Let m(v) be the subscript such that 


max [pF (0 | a;) | = bn (0 | Gn (33) 


and let 6!(v) be any arbitrary decision rule having the 
same rejection rate as 6**. It is to be proved that 
P01) 2 P08). . 

For every 6!(v), a decision rule 6?(v) can be con- 
structed as follows: 

For every 2, 


52(do| v) = 8'(do| v) 
52(dm |v) = 1 — 6'(do| v) = D> 8%(d;| »v) (34) 
=) 
52(d; | ?)=0  foralla 4 0 # mm, 
Evidently, 
P67 PG) S=1 0 (35) 
and 
P64) = >> é'(d; 0) pil*(v | a;)dv 
V 4=) 
Se f Ss 5'(d; | 2) Pl’ (v | Am) av 
Yat 
= ca 2) Pak (0 | Am) dv 
Vv 
= P,(8?). (36) 
It follows from (35) and (36) that 
P.(5’) < P.(8'). (37) 


That is, 6? is better than 6! (or at least as good) in the 


sense that for the same rejection rate 6? has an error 
rate smaller than, or equal to, that of 61. 


The next step is to show that P.(6**) <P.(6*). As 


shown in (22) ad (23), the decision rule 6** partitions 
the observation space V into two nonintersecting re- 
gions, Vo** and V— Vo**, so that for every ve Vo** 


bnF(v| am) < BD) piF(v 


t=] 


a;) (38a) 


5** (dy |v) =], (38b) 


25-1 
and for every ve V — V)** 
Punk (v| Be Jee Oe % Pla 1;) (39a) 
i=1 
d**(d,,| v) = 1. (39b) 


Let Vo* be the largest subspace of V such that 5°(do| v) 
is nonzero for all v belonging to V2. Vo? is not properly 
contained in Vo**. This follows readily from the condi- 
tion that P,(6**) = P,(6?). The latter may be written as: 


f ae pF (| a;)dv 
Vet —Ve"Ov ial 


0 


++ [1 — 8°(do| v)] D> pr(v| aide 
Ve ayo: cl 
= 52(do| v) >> pF (0 | a;)dv. (40) 
Jy areny?e oe 
Substitution of (38) and (39) in (40) gives: 
il Pml (v | Gn) av 
era ella 
te [1 — 5°(do| 0) | Pmi*(2 | Am) av 
OTe ll a eo 
<[ 6°(do 2) PF (2 | Am) av. (41) 
Vievior Vion ; : 


The equality sign prevails if, and only if, PnF(v| dm) 
isequal toB >-*_, pF (o| ai) throughout the region Vo**U 
Vo?. 

The probabilities of correct recognition of 6** and 82 
may be written respectively as: 
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P(5**) = if Punk’ (v| an)dv 
V—Vy"" 


I 


[ Pnl(2 | Gm) dv 
ye Wo era) 


+f Pnk (2 | a,,) dv, 
Vig Vora 


by 52(d; | v) pF (0 | a;)dv 


Vi wl 


={ Punk (0 | Gm) dv 
VW" U9") 
+f PnF (| Am) dv 
Mnpi=Vio ellos 


+f 5°(dn | 0) Pnl (0 | Gm) dv 
Veni ieas 


(42a) 


oa 5°(dn |v) Prd? (o | a,)dv. (42b) 

asian alae 
In accordance with (42), (41) is equivalent to P,(6**) 
<P.(6°). Proof that P.(6**) <6P.(6!) is thus completed. 
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Errors in Electronic 


tial \nalyzers 


{— Bandwidth Limitations’ 
PAUL C.. DOW, Jat 


Summary—When a differential analyzer or analog computer is 
set up to solve a given equation, certain errors are introduced into 
the solution because the components of the computer are not perfect. 
Stated another way, the computer produces the solution to an equa- 
tion, called the ‘‘machine equation,” which differs from the given 
equation. In this paper the computer imperfections considered are 
1) operational amplifier frequency response, 2) capacitor leakage 
resistance, and 3) stray capacitance in summing and integrating 
amplifier circuits. It is shown that when the given equation is a sys- 
tem of one or more linear differential equations with constant co- 
efficients, the machine equation can be expressed approximately as 

an equation of the same degree as the given equation and with con- 
stant coefficients which are functions of the coefficients of the given 
equation and of the computer imperfections. 


INTRODUCTION 
[i has been shown! that bandwidth limitations 


inherent in summers and integrators used in dif- 
ferential analyzers produce shifts in the character- 
istic roots of the differential equation being solved. 
The methods of analysis used to obtain these shifts de- 
pend for their simplicity on particular computer setups. 
In this paper no restriction as to the method of computer 
setup is imposed. Furthermore, the purpose here will 
be to find not the errors in the roots of the equation but 
rather the errors in the constant coefficients of the dif- 
ferential equation, whether or not the roots are known. 
The bandwidth limitations of summers and integra- 
tors have been considered by a number of authors* and 
are briefly reviewed here. The transfer functions of 
summers and integrators are developed and then used 
to find the exact equation solved by the computer. This 
in turn is reduced by certain approximations to an equa- 
tion of the same order as the given equation. 


* Manuscript received by the PGEC, June 19, 1957; revised 
manuscript received, September 17, 1957. This paper contains part 
of the results of a dissertation submitted i in partial fulfillment of the 
requirements for the Ph.D. degree at the University of Michigan, 
Ann Arbor, Mich. 
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OPERATIONAL AMPLIFIER 


Many computer amplifiers have frequency response 
characteristics as shown in Fig. 1. The open loop trans- 
fer function of such an amplifier is 


Go 
aie 


G(s) = (1) 


If wo is the frequency at which the open loop gain is 
unity, then | G(s) | =1 when s=jwo. From (1) one has, 
if Go>l, 


we: (2) 


= 7—-—-6 db/octave 
= 
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O 
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Fig. 1—Frequency response of open loop amplifier. 


It is shown in the Appendix that even those ampli- 
fiers which closely sees las the frequency response 
characteristics of Fig. 1 at high frequencies but not at 
low frequencies can still be represented by a first-order 
transfer function such as given by (1) for the purpose of 
this analysis. 
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SUMMER 


The amplifier connected as a summer is shown in Fig. 
2, in which the stray capacitance and resistance to 
ground are shown in dotted lines. 

If Z; is the impedance between each input terminal 
and the amplifier summing junction P, and Zs is the 
feedback impedance, and there is no grid current or 
amplifier unbalance, one has 


ites ewe V 


‘See. 0 7m ay 
_ es te Oe (3) 
i Zi Zs R, 
i 
and a - 
Vo = — G(s)V. (4) Fig. 2—Amplifier connected as a summer. 


Eliminating V from (3) and (4) yields 5 ae EY ’ 
In most computer applications it will generally be true 


va) ; 
t oe 2 V; that 
t Z; law) 
ee 5) alt 
te oll Ge bose S| p Ry R 
Bee LehiG.2 Sb R 
G(s) f fu Aree ty BR RC + (1 ae > eee =) 
Go frente ees 
From lig. 2 it is seen that 
ys? << 1 (11) 
R; 
“4; = Coe (6) and (9) can be written 
R : = ; R; 
= (7) (1+ sR,C) = Vv, | 
1 SRyC; i OR; 
é v= = eres} (12) 
Substituting (1), (6), and (7) in (5) yields 1 + Bs 
=< RA + SRC; 
_ Ea = a 
i i SANGO 
if —— ee (8) 


0 


WRC, Seite Re Ok ee sR,C;) 


For economy of notation, let R;C,; foreach ibe the same. The final expression for the summer is written 
(If this assumption does not apply, the following equa- é; 
\ Pp PPly § eq T;s Al { R; 


tions become more lengthy but the same method of Ves pS. V; (13) 
analysis can be used.) Then (8) becomes Las the laa Rs 
R; where 
—(1 ++ sR;:C;) ns, ae V; ; 
Tt R; IPS oS Or (14) 
Vo= ay ve (9) 
a + Bs + ys" . and 
where "fi Rees | 
r= (1+ D4 2) + Re, (15) 
1+ (1+ O24 A) Go Fidei tod 
a= = cag eee ‘ 
Go t R; Rs 
INTEGRATOR 
Bi Riek: 
Bre Rete i The amplifier connected as an integrator is shown in 
Go 2 oR: R, nae E bee ; apy 
Fig’ 3, in which the stray capacitance and leakage re- 
Le a Ry sistance are shown in dotted lines.? I°q. (9), which was 
ae ee R; v3 ! R ‘Cs fe i A _ - : z . . , 3 ; - 
1 a Crt RCs + RL ye =) obtained for the summer, will also apply for the inte- 
To ; R; 
y= —[ 4/C;+ £C,4+ RC; Ve —], (10) 7 In a second article the effect of dielectric absorption and com- 
Gy Sale ie plex permittivity of the dielectric in integrator capacitors will be 


considered. 


Dow: 


Fig. 3—Amplifier connected as an integrator. 


grator if Ry and C; are replaced by Rz and C, respec- 
tively. Then 


V; 
1 -++- SRC; = 
CEs 


Vo=— 3 (16) 
: a’ + B's + y's? 
where 
1 1 1 1 1 
(natin pi ae te ~( = =a) 
a zal =) Ge 2X ee RG 
a ah rea | { 1 
/ ee 1 a if : ae = 
B as a (ee 5 2, Ries RG 
Patt aay RC, nie 5 
oe o>) 
ra P(14S4RGLD ge) (17) 
eee SER 


In most computer applications it will be true that 
(17) can be written approximately 


i 1 1 
oe + 
oe rae po RC 
Bit 
ie 
io — (18) 
+ Go 


and (£6) can be written 
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—(1+ sRC) 2 ae 


ke CONG eee ae 3 
T eet 1 ) 
TS eh ad COW ica Radial 
& ) are pee RC 


The final expression for the integrator is written 


(19) 


Dy (7; s+ 1) a 
ne ae 
Tea: re 1)(Ts x) ae 1) 


(20) 


Differential Analyzer: 


where 


ii 
jis ae 
Go 
T, = RC. (21) 


SOLUTION OF A SINGLE DIFFERENTIAL EQUATION 


The equation to be solved, 
equation,” is 


called hereafter the “given 


(22) 


(See 


k=0 


Although a particular computer setup will be used to 
solve(22), it will readily be seen that the method will 
apply equally well to any computer setup. Fig. 4 shows 
the computer setup used in this analysis. 


Fant 
baat 


Fig. 4—Computer setup for linear differential equations. 


To determine the equation which the computer is 
solving, the voltage at y is evaluated in terms of y’, then 
y is set equal to y’ by connecting the two points to- 
gether. This yields 


agTim(Tss-+ 1m? 
beastie Tas tl)(Tisbiy(Past 1)” * 
ag iTi( Tis HI) 
+ Tas) (Post l)(Pist1) (Pest) | 
i979) 
(GT eT et) 


(23) 


where 7, and 7. are the integrator time constants de- 
fined by (21) and are assumed the same for each inte- 
grator to keep the equations reasonably compact. 7’ 
and T,. are the time constants of the first and second 
summers, respectively, defined by (15) and T;=R,C;. 
5=1 or 0 if the coefficient a, of the term in which it ap- 
pears is determined by a resistor connected to the input 
of the first or second summer, respectively, or in the last 
term if f(f) is an input to the first or second summer, re- 
spectively. 
Noting that 


y'Ti(Ts + 1) 


ee 1)"(Tes + 1)" 


(24) 


Bandwidth Limitations Dy 


(23) can be rewritten as 


= 1 k 
Gl Ss | (Py 1)P ere 
Lay pee ( =) 


( :) gaa (BT ns + 1)(Tyas + 1)(Tas + 1)" 
f(s) (Tis -{- 1)tits 


SS : i se 5 

(67.35 -+ 1) (Tos zi 1) (Tes “T: ie , 
Iq. (25) is the exact equation solved by the analog com- 
puter setup shown in Fig. 4. It is seen that the charac- 
teristic equation is at most of order 2m-+2 in s. The 
given equation has a characteristic equation of order m 
in s. Thus m-+2 extraneous roots are introduced by the 
m integrators and two summers. Macnee has shown? 
that there are m extraneous roots 


: (26) 
s=—-- 2 

T> 

and two roots 

i 
5= —- 

dS 

(27) 
Ay =—_- — -—— eo 

Lins 


Since they are large and negative they will damp out 
quickly and have negligible effect on the solution. The 
remaining m roots, if the errors introduced by the im- 
perfect components are small, are approximately equal 
to the 7 roots of the given equation. 

Ifs;(@=1, - - +, m) are the roots of the given charac- 
teristic equation and if 


1 1 is 1 
eS ? ) ) 
Ta Ts2 T> hee 


(28) 


1 
—<|s; 
T; 


then the following approximations will be valid in those 
regions of the complex plane where s&s;: 


1 
1 — Tas 
Tis + 1 


(1 1 ye 14 m 
Ps a TS. = Tis 
(1 +- Tsis)(A -L T 25) = 1 ot. CT +- T2)s, WEF (29) 


If these approximations are made, (25) becomes 


MW m—1 
E (1 +") + >> ska, i — 67 1s — T,25— (m—k) Tos 


18 k=0 


B\s 
+(m—k+1-+8) rats | me 


18 


= 1-67 .15— T.25—mT 25+ (m+14+-5) Tis} f(s). (30) 
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This can be further simplified to 


ua k+ 1 ; 7 
~ s* ap + Og ae a OL ent T's 
k=0 L 


1 
+ (m—k+1)T, —(m-k+2+4+8) ri Ja 


= f(s) — {8Ta + Tit m2 — (m +14 8)Ti}sf(s) (31) 


where it is understood that a, = 1, dui =0, a= 

Since a feedback resistor in Fig. 4 will be connected 
to the first or second summer depending on the sign of 
a, and on the order of the equation, the following gen- 
eral rule will apply: 


“ = a, if (k + m) is odd 
a; > 0 i 
= 0 if (k + m) is even. (32) 
ays = O if (k + m) is odd 
ag ; ( ) : 
= a; if (k + m) is even. (33) 


The coefficients of (31) can be computed and compared 
with those of the given equation. By this means one can 
determine the effect of computer errors on the given 
equation even if the solution is not known. It would 
also be possible to perturb the coefficient-setting resis- 
tors or potentiometers by an amount equal to the error 
in coefficients given in (31) and to observe the effect of 
these perturbations on the computer solution, thus de- 
termining the sensitivity of the solution to these errors. 


SOLUTION OF SIMULTANEOUS DIFFERENTIAL EQUATIONS 


The given set of equations in matrix notation is 


[%; ] = [axe] [xe] a [fs] 


X= AX+4+F, 


(34) 
(35) 


The computer setup to solve (34) is shown in Fig. 5, 
where it is assumed that all the coefficients aj. are 
positive. 

Using (13) and (20) it is possible, as before, to write 
the exact equation being solved by the machine. (Note 
that if any of the aj is negative, the corresponding po- 
tentiometer would be connected to the input of the 
integrator instead of to that of the summer. The symbol 
6 is used as before; 6=1 when the input is to the sum- 
mer; 6=0 when it is connected to the integrator.) Then 
the voltage x; is given by 


[-e)] a | T33054(T ys + LyFr? | lx 
‘ie (1+ oF s\f 4 Piss 
eee ee | (36) 
(1 4287 sl Test Pagal) 


where 
71;= 7; for the jth integrator, 
’.;= 1, for the jth summer, 
1’, = T> for each integrator, 
*:= R;C; for each summer and integrator. 
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Fig. 5—Computer setup to solve simultaneous equations. 


Eq. (36) is the exact equation solved by the computer 
and is seen to be at most of order 37 if all 7 summers are 
used. That is, it is of order 2” plus the number of sum- 
mers. 


Using approximations (29), (36) can be written 


ee 25 T 8 2) Past (1 8) Ts} 


peo) a 


+ [ffi 87s — Tes + (1+ 8)Ts}]=0 (37) 
where 


Te 1, oy 
0, iff 8. 
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Finally (37) becomes 


Te 
lualoT$ Te (148) Halll | 0-2 | ba 


4 [f= feTytTe—A+HT3 fil. (38) 
Let 
B= [an {67.3 + To - A+ 47 be Lin] 
Lx 
(G; —= E sae | 
Sas Vj 
G = [f,— {8% 4+ T2- A+ 5) 73} f;]. (39) 


Then (38) can be written 


a= 6 eG: (40) 


Eq. (40) is a set of first order simultaneous equations, 
To facilitate comparison of 
the two equations, both sides of (40) are premultiplied 
by the inverse of B. Then 


as was the given equation. 


MS BoC BIG (41) 
The matrix B as given by (39) can be written 
Bese (42) 


where J is the identity matrix and the elements of the 
matrix A are small compared to unity. The inverse of B 
is given by 

Bia l= Aj he (43) 


and if the terms involving powers of A higher than the 


first are neglected in (43), one has 
Bo} = (T= oy OTs = Ora 


Substituting (39) and (44) in (41) and neglecting 
ucts of small terms yields 


(44) 


prod- 


Tx 
[x] = E aa = i’ AjmOme{ OT 5; a2 T> a (1 —- 5) rs} |fo) 


ant dTy; + T, (1 Tite 
ari {87's; ples (1 =f 6) Tih fil (45) 


where the subscript ™ appearing in each term of a prod- 
uct denotes summation over m. 

Eq. (45) is identical in form to (34) and the two can 
be compared term by term to determine the changes in 
coefficients resulting from computer errors. 

Although only two computer setups have been con- 
sidered here, the method of approach can be applied to 
any setup used to solve linear differential equations to 
determine the errors in the constant coefficients. 


EXAMPLE 


As an example of the application of the method, con- 
sider the equation of simple harmonic motion 


% + wx = 0 (46) 


Fig. 6—Computer setup for simple harmonic motion. 


Fig. 6 shows the computer setup, which is like the one 
shown in Fig. 4. 
from (31), and is 


The machine equation can be found 


-f- "| “oF OE 


p 
E tie nee ole dp eo hs (47) 
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Iq. (47) has a damping ratio given by Substituting (49) in (5) gives, for a summer, 
> 


wo Ly 
= ee V; 
OY hee tS) bp , oe re =, J 
1 iP w(7', 2s vs) (48) 0 ae (Tas = 1)? (1 : > zr) 
S EES ab fe a 
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where the stray capacitance and resistance to ground 
have been ignored since they do not influence the result. 


| (50) 
t= — 
wily 


where 7;, 7;, 7, and 72 are defined by (14), (15), and 


(21). Eq. (50) is rewritten as 
R; 
—Go(Tes + 1) =v; 
Loe Re Og Cis Boos — mine: (51) 
\yf es a —, soe . 
(1 ee | ey j2ra(s aE =) i GT | + |G. ne dete =| 
I’q. (51) can be written 
R; 
= SEY: 
oe 
Ko > a ; — : (52) 


APPENDIX ae 


Since many amplifiers in use actually have more com- 
plicated gain characteristics than indicated in Bigs 1 it 
is important to provide for the analysis cf computers 
using such amplifiers. Since all ‘amplifiers must have 
high gain at low frequencies (for good static accuracy) 
and must fall off at approximately —6 db per octave in 
the vicinity of unity gain (to insure stability when used 
as summers. or integrators), Fig. 7 shows a frequency 
response characteristic which can fairly well be used to 
represent a wide variety of operational amplifiers. 


20 log G = 
e hi - 12. DB/OCTAVE 
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w= /T, 
LOG FREQUENCY 


we!yT 
Fig. 7—-Frequency response of open loop amplifier. 


The transfer function for the amplifier of Fig. 7 is 


Ca Pas Ga 
Coe 


(Tas —+ 1)? (9) 
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If 


bee E es air “ECTS = er eee || 


|Go ven; 


(52) becomes 


Ry 
ks dig 7 
* 2. Ri : 
fp = 53 
illite, oy 
where 
2T4—Tp 
ree (1 > =) (54) 
0 7 7 


From (15) it is seen that this is equivalent to replacing 
the operational amplifier of Fig. 7 by one having the 
characteristics of Fig. 1, where 


Po = 274 —Ts, Ripe! Ae (55) 
It can be shown by a similar procedure that (55) applies 
also when the operational amplifier is used as an inte- 


grator. 
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Syn thesis of Vector Networks” 
R. E. HORN¢ anv V. G. FAUQUEt 


Summary—The convenience of vector notation in formulating 
physical geometrical problems results principally because the sig- 
nificance of the problem can be isolated from the analysis used in its 
solution. When analog computer techniques are employed in prob- 
lems of this nature, the advantages of the vector methods frequently 
are lost because the computer inherently is more adaptable to solving 
problems described in the algebraic field of real numbers rather 
than in a vector space, and vector equations must be reduced to their 
scalar counterparts before a network for solving the equations may 
be synthesized. To facilitate treatment of problems of this type, a 
method is presented for synthesizing networks directly from the vec- 
tor notation. This method will simplify synthesis and analysis of the 
networks by drawing a closer analogy between the mathematics and 
the electronics. 

The application of transformations and operators represents two 
areas in which the pertinent aspects of the method can be illustrated. 
This is accomplished by presenting first the basic mathematical 
background, followed by the network aspects of the problem, and 
finally, by the illustrative examples. The more extensive problems 
usually encountered, for example, in airborne fire-control systems, 
then may be synthesized by proper application of the elementary 
“vector networks.” 


-]. TRANSFORMATIONS 
(aaa two: sets of three-dimensional Car- 


tesian coordinates. Let X41, Xo, X3,and V1, Yo, Y3 

denote mutually orthogonal axes of the two sys- 
tems. Bases for vectors in either system will be ortho- 
normal and colinear with the coordinate axes. These 
bases shall be denoted by 4, f, 73, and ji, Jo, 73 in the 
X and Y systems, respectively. Let the i and 7 bases be 
related as follows: 


ray = fis + fisse + fisJs; 
12 = fot + foofe + fost, 
13 a fait = faeJe + fasts. (1) 


In matrix notation, this relation is expressed as 
11 fi fis aes pal 
ig) = for for for) doy = (fall) je y @) 
gal i fe fale Ea 


Consider:a vector r such that 


r= gilit gole + gsls, (3) 


and 


t= Wij + hoje + ajo. (4) 


Since a vector defined by a complete and independent 
set of bases has a unique set of components, the trans- 
formation relating any two sets of bases also relate their 
respective components. Therefore, 


* Manuscript received by the PGEC, June 19, 1957; revised 
Manuscript received, October 1, 1957. 
} Westinghouse Electric Corp., Baltimore, Md. 


ey hy 


a = [ll } ©) 


Multiplying the left-hand side of (5) by the inverse 


: | * 5 
transformation [LF 1 yields 


erated 
12) = [fully ge 1- (6) 
ee | 


Since the IL Fes transformation is orthogonal, its inverse 
is equal to its transpose. Therefore 


fis fei a 
fill* = 1 faa fea foe 1s (7) 
Ie foes a 


and (6) may be written 


| ea for far 
ho} =1 fio foo fae 24. (8) 
Li. ie fos ie a 

IT. OPERATORS 


Suppose that the rate of change of the vector r is to 
be expressed in terms of its components on the f base 
system, which is rotating. That is, 


dr = my, + moto + mals, (9) 


where d is the operator d/dt! and r=r(t). With r defined 
by 


r= gil, + goto + gals, (3) 
the derivative of r is by definition 


dr = (dgi)i1 + (dgo)i2 + (dgs)tz, + @ X 4, (10) 


where 
(11) 


defines the angular velocity of the X-coordinate system 
with reference to some fixed system. Using conventional 
vector methods, 


(ON = 114 a0 Wolo ap W313 


Sle ee ee 


1 The operator d is not to be confused with a multiplier. The 
derivative operation is not distributive with respect to multiplica- 
tion. 


_——_— ae ee|)h ll CU 


which may be alternately expressed as 
dr = W28y)1y = (ws “| dgs ne 123) 12 


ay dg3)13. 


(dgi — wsg2 + 
++ (woe -| Wf (13) 


From (9) and (13) the equivalence 


my dg — wage + woes 
my | = wseit+ dgo — wg: |, (1-4) 
Ms —wee1 + wigs + dg: 


or the factored form 
3 Wy §) 


(15) 


Fils = W3 d coe £2 
m3 poe ®) d §3 


may be recognized. It should be noted that (15) also 
may be derived directly using matrix methods.? 
The matrix 


d —ax We 
ws d —w | = || Dall (16) 
— Wy | d 


shall be designated by the term “operator matrix.” The 
product of two operator matrices will provide a second- 
order operator matrix. In general, || Deol] » will be an ath 
order operator matrix, where the superscript indicates 
the number of multiplication of the operator matrix by 
itself, / 

For use in the synthesis of networks, Section IV 
shows that the matrices need not be combined before 
synthesis. 

To differentiate a vector in a fixed coordinate system, 
the operator matrix (16) may be used, but since the w’s 
are zero, the operator matrix becomes 


d 0 g 
0d -G.= |p (17) 
00d 


The integration matrix is defined as the inverse of the 
operator matrix. Similarly, the expression 


fa = dr (18) 
is defined. Now let 
dr = qty + gots + gots (19) 
Then 
eel 
2 | = [bal (20) 


*R. A. Frazer, W. J. Duncan, and A. R. Collar, “Elementary 
Matrices,” Cambridge Univ. Press, Cambridge, Eng pp! 
1938, 
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or, since the matrix || Doo] always possesses a right in- 


Verse, 
q1 £1 
|| Del] } 2) = | ge, (21) 
Ala 


which may be solved for IIal| When 
IIql| are known. 


the initial values of 


III. NEtworxs 


Synthesis of equations of the general form 


f / 
€11 €12 €18 Ry €11 12 eH hy 
@21 €22 C23 £21] = | eo)’ C22" €23" hy (22) 
€31 €32 €33 ie C31 €39" aya hs 


will be considered, since a correspondence exists be- 
tween (22) and any equations discussed previously. The 
network which approximates? the analog of (22) will be 
general in the sense that any illustrative example will 
be a special case of this configuration. The basis for the 
synthesis is the one-to-one correspondence between (22) 
and idealized admittance equation. This idealized ad- 
mittance equation may be written as 


Vile Vie V3 
ite ¥22 Yo3 


Voi, V3o Van 


oe , / , 
— Wy, Yin Viz Vig Uy) 


hal yor’ Yoo! Vox" ito |, 
ne Yar ¥30" 33 bie 
where the coefficient matrices Il || and Io’I| represent 
network admittances, the w’s the independent voltage 
variables, and the w’s the dependent voltage variables. 

Treatment of network equations of this type, as ap- 
plied to analog computers, has been covered elsewhere,*® 
and therefore, it should suffice to discuss briefly the 
schematic representation to be used. 

The basic network consisting of an amplifier with gain 
ux having bipolarity inputs, and associated admittances 
Y;, is shown in Fig. 1. Input grid voltages are repre- 
sented by w’s and output voltages by +w, and —w,. 
By representing the network in this way rather than as 
a conventional computer schematic, recognition of the 
basic matrix form is simplified. 

Fig. 2 gives the final representation, as its combina- 
tion with additional networks of the same type. The 
principal changes have been to draw the amplifier as 


numbered squares and to eliminate the representation 
of the common reference mode. The schematic for the 


(23) 


3 Networks which are exactly analogous to the equation are not 
realizable since the gain of the amplifiers would necessarily be non- 
finite. A detailed description of this problern is given by Honnell and 
Horn, footnote reference 3. 

4P, M. Honnell and R. E. Horn, 
matical machines, 
ber, 1955. 

* P. M. Honnell and R. E. Horn, “Analogue computer synthesis 
and error matrices,” ATEE Trans., Part I—Commun. and Electronics, 
vol. 75, pp. 26-32; March, 1956, 


“Matrices in analogue mathe- 
’ J. Franklin Inst., vol. 260, pp. 193-207; Septem- 
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Fig. 1—Basic computing ne*work. 


prescribed voltage sources --w, and —u, and corre- 
sponding admittances Y,,, although not drawn in Fig. 2, 
is represented by similar grid structure. By arranging 
the schematic as a rectangular grid, the direct corre- 
spondence between the matrically-formulated problem 
and the computer network is emphasized. 

The general form of the equations considered in this 
paper is given as (23). Using the notation described in 
the preceding paragraphs, the synthesis of (23) is illus- 
trated in Fig. 3. Required amplifiers are shown as num- 
bered squares, with horizontal lines (v’s) representing 
input grids. The vertical lines indicate measured vari- 
ables, those in the left section of the array representing 
bipolarity amplifier outputs (w’s), those in the right sec- 
tion representing prescribed voltage sources (2's). Ad- 
mittance elements are connected in positions corre- 
sponding to terms of the problem matrix and, although 
only positive matrix terms are illustrated in Fig. 3, 
negative signs are obtained readily by utilizing the nega- 
tive voltage variables. With the polarity convention 
indicated, principal diagonal terms on the left side of 
the network equation (23) normally must be positive 
for network stability. 

To synthesize networks corresponding to the trans- 
iormation and operator matrices already described, it 
is only necessary to relate each term in the coefficient 
Matrices to the corresponding network parameter in the 
computer. This will be illustrated by specific examples 
in Section IV. 
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Fig. 3—Network corresponding to (22) and (23). 


IV. EXAMPLES 
A. Transformations 


Transformations through Eulerian angles are the 
common types requiring synthesis. For example, rota- 
tion of the X system about the Y, axis through an angle 


of @ (in a right-handed sense) requires the |/f;;|| trans- 
formation of Section I to be 
fu fie fis eas) 0 
for fer fos{=;,0 cos@ sin ¢}, (24) 
far fso fas iP —sin ¢@ cos¢ 


so that (5) becomes 
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Fig. 6—Network corresponding to (21). 


, 


The corresponding network equation® 


is 0 0 ies il 0. 0 14 
OL Os, = 10 cosd sin 64) wp 

Me Oot | 0 —sin @ Cos ilies 
then may be synthesized directly as indicated in Fig. 4. 
Note that the formal writing of (26) is unnecessary since 
the network of Fig. 4 is obtainable directly from (25). 
The varying admittances required for the variable ¢ are 
conventionally obtained by the use of servodriven po- 
tentiometers and series conductances which yield the 
required transfer admittance. 


Frequently, transformations involving three Euler 
angles are required. These take the form 


ana 
Ll Lo 


(26) 


0 0 lees 0 —sin #6 
cos@ sin ¢ Os al ) 


—sin¢@ cos d_ ie 60 cosé 


cosy siny 0 [i 
(27) 
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5 The ois 010 i the identity matrix. It is used here so 
001 


that this example will conform to the general form described in (28 
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Tig. 7—Network corresponding to (32). 


In this case, three networks of the type shown in Fig. 4 
are connected in cascade. Of course, this is equivalent to 
multiplying the matrix by the Heal matrix, the product 
by Ge and that product by the |o|| matrix. Since the 
form of each network has already been established, it is 
unnecessary to perform the multiplication indicated in 
(27), but only to synthesize the three networks with 
terms as required. This process not only simplifies the 
initial mathematical problem, but it provides a simpler 
synthesis and analysis of the network. 

The inverse transformations, as discussed in Section ie 
have a similar form to the transformations themselves. 
Therefore, the synthesis of networks involving inverse 
orthogonal transformations reduces to the identical 
process as that explained above. This avoids the ques- 
tion of stability which necessarily arises for networks 
involving “cross coupling” between amplifiers. 


B. Differentiation 
Differentiation of a vector, described by 


1, 


[ d —w3 Poe 
Me | = @3 Gd ~arfy gol, (15) 
Es ex 4 d | ee 
and the computer equation 
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may be synthesized as illustrated in Fig. 5. The one-to- 
one correspondence between terms in (15) and network 
elements is evident again. 


C. Integration 


Integration of the r vector in the X-coordinate sys- 
tem, as discussed in Section II may be expressed as 


fic Caen W2 BY ‘| 

W3 d ue SS oh 
a z | 
representing an inverse of the operation on the |!g!! 


matrix in (15). The network corresponding to (21) is 
shown in Fig. 6. 


(21) 


— We Wy d 


D. Filtering 


The problem of filtering “noise” in a moving coordi- 
nate system is illustrative of a simple practical applica- 
tion. As previously defined, let 

tee Pity ae Q2l2 ae £313 (3) 

represent the vector input to the filter and 
P= pilit pols + pals (29) 

represent the output vector. Then the vector equation 
(1 + Tid)p = (1+ T.d)r (30) 


is representative of the vector “lead-lag” network to be 
synthesized. In matrix notation, (30) may be written 


M7] + Till Dll }llell = (lz + Tell Del} llell 1) 
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with the corresponding network illustrated in Fig. 7. 
In (32), an insight into the stability question may be 

obtained if the constant coefficient case is assumed, re- 

sulting in a determinantal equation having roots 


1 
\teer = 33 
1 a (33) 
1 
2,3 = — T. ai jVer + we? + w3*. (34) 
1 


Eq. (34) indicates the cross-coupling effects introduced 
by the rotating coordinate system. 


V. CONCLUSION 


By indicating the direct relationship between com- 
monly-used vector equations and matrix synthesis meth- 
ods, a convenience is achieved that may be lost by con- 
ventional methods. Specific examples illustrating trans- 
formation and derivative matrices are representative of 
the type problems that arise, and the corresponding 
basic configurations serve as fundamental “vector net- 
works” to be used as units in the synthesis of more com- 
plex problems. 


Switching Functions of Three Variables* 
D. W. DAVIEST 


Summary—A switching function is a function of variables which 
take only the values 0 and 1, and which takes only these values it- 
self. There are 256 different switching functions of three variables, 
but only 218 of these really depend on all three variables. 

A switching function of three variables can be expressed in terms 
of switching functions of two variables. For example 

FPA, F(B, C)], FCB, C)} 
can be shown to represent any function of A, B, and C if Fi F.,F; and 
F, are suitably chosen switching functions. 
The problem solved is: for each switching function of three vari- 
S, what is the least number of switching functions of two vari- 
S required to express it? 


ible 
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This problem leads to a discussion of the symmetry of switching 
functions. Expressions which, like the one above, can represent any 
switching function of three variables are investigated. Expressions 
which, by permutation of variables, can represent any switching 
function of three variables are also determined. 

The investigation is done by means of “‘logical diagrams,” which 
give a better intuitive understanding than the functional expressions. 


N THE logical design of digital computers, electri- 
cal signals which have only two signifcant levels 
are often represented by variables which can take 
only two values 0 and 1. Let A, B, and C be three such 
variables. Suppose now that the three signals repre- 
sented by -l, B, and C enter an unspecified unit and this 
unit produces an output which depends only on the in- 
puts and which also has only two signi®cant levels, 


represented by 0 and 1. Such an arrangement is shown 


in Fig. 1(a). The output signal is expressed as 
F(A, BC} 


which is called a switching function of the three vavi- 
yee doy een al (Oe 


A A 
ee eer 
LS B 


(a) (b) 


Nomenclature for “black boxes.” 


ables 


Fig. 1 


An unspecified unit such as we have described is often 
called a “black box.” We are reserving the term here 
for a unit with no memory; 7@.e., one in which the output 
depends only on the present values of the inputs. 

It is well known that a switching function of # vari- 
ables can be synthesized by using an arrangement of 
several black boxes each with only two inputs, like that 
of Fig. 1(b). 

The question answered here is: for each switching 
function of three variables, what is the minimum num- 
ber of black boxes (of any type) with two inputs which 
can be connected together to synthesize it? A related 
question is: given a circuit arrangement of black boxes 
with two inputs, the whole having three inputs and an 
output but nothing being known about the boxes them- 
selves, what is the class of switching functions the given 
function of three variables is known to belong to? This 
is answered for all cases of four, boxes or less. 

Since black boxes with two_inputs are important to 
our problem we shall adopt the convention that “black 
box” means one of these. 

The problem could be stated more accurately thus: 
for each switching function of three variables, what is 
the least number of switching functions of two variables 
required to express it? As an example, take the circuit 
arrangement of Fig. 17(a). This can be expressed as 


PCAPE, OC) oF PAs OFC, Bh. 


The three boxes correspond to the three functions Fy, 
Fy, and /’;. In this paper we use the logical diagrams to 
represent this type of equation. 


FUNCTIONS OF 0, 1, AND 2 VARIABLES 


The following remarks concerning functions of less 
than three variables are perhaps obvious to anyone con- 
versant with Boolean algebra or its related subjects, but 
they are given because they can be developed later for 
three variables. 

There are two functions of 0 variables, 0, and 1. 

A single variable can take two different values. If out- 
put, 0 or 1 is specified for each input we can have four 
functions of a single variable. These would be 


F(A) = 0,1, 4, ~4, 


where ~l is the negation of A. 
A pair of variables can take four different sets of 


values. If the output, 0 or 1 is specified for each input, 
we have 2'=16 functions of two variables. Table lisa 
list of the variables of F(C-1, B), giving the four values it 
takes for the cases (4, B) =(0, 0), (0, 1), (1, 0), and 
(1, 1). The functions have been numbered 0-15, a nota- 
tion which will be used later. 


TABLE I 
NOMENCLATURE FOR SWITCHING FUNCTIONS OF Two VARIABLES 


9 } 10) 11/12) 13] 14) 15 


i|Bloj1ila{afalslol7]s 

0/0/0 E folilojiloliloliljoliloliloli 
0 | 1/0. o|1 4 (0/0 ‘titfololil1lolola} 1. 
1/o/ol/ojoj/o|1li}a{i/olololol1|i |i} a0 
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It is obvious some cases listed as F(A, B) are really 
functions of just 4 or B, or of neither. Functions num- 
bered 0 and 15 are functions of no variables, and func- 
tions 3, 5, 10, and 12 are ~A, ~B, B, and A, respece 
tively. These are called “degenerate.” There are 10 
nondegenerate functions of two variables, two non- 
degenerate functions of one variable, and two of no 
variables. 

The particular black box A F(A) which gives 
F(A) =~A, will be represented -> in what follows. The 
only other nondegenerate function of one variable if 
F(A) =A, which does not need a symbol. 

If F(4, B) is a general function of two variables, then 
so are 
~F(A, B) 
~F(~A, B) 
~B(A, ~B) 


~F(~A, ~B). 


F(~A, B) 
FUAL OB) 
F(~ A ) ~B) 


Therefore, symbols such as 


risk abbeys maa 


if the function is unspecified, can all be represented as 


Up 


in which the negation has been absorbed into the black 
box. Suppose now, that a supply of negators is available 
how many black boxes with two inputs are needed to 
construct any function? 

It is easily seen that functions 1, 2, 4, 7, 8, 11, 13, and 
14 are related by negating 4 or B or F or interchanging 
A and B. : 

The remaining nondegenerate functions, 6 and 9, are 
related to one another. 

By making use of the ++ symbol in what follows, we 
need employ oaly two different functions of two vari- 
ables, one representative of each class. We choose sym 
metric functions which we symbolize as follows: 


represents function 8, 
(intuitively 4 and B), 
represents function 6, 


(intuitively 4 not equivalent to B). 


It is convenient, however, to have symbols for two 
other functions which appeal to intuition namely 


representing function 14, 
(intuitively A or B). 


representing function 9, 


(intuitively A equivalent to B). 


REPRESENTATION OF FUNCTIONS OF THREE VARIABLES 
A set of three variables can take eight different values. 

The function is specified when its value is given for each 
of the eight values of the inputs. There are, therefore, 
28=256 different functions. This includes the following 
degenerate cases: 

0,1 

EA Bet 1s, ig mC, 

all the nondegenerate functions F(A, B), F(B, C), 

F(A, C) of two variables. 

These are 38 in numbers, leaving 218 nondegenerate 
functions of three variables. 

The set of points (A, B, C) in Cartesian coordinates 
given by A, B, and C taking the values 0 and 1 is the 
corners of a cube. A representation of F(A, B, C) is 
given by attaching the value of F to each corner. For 
example, Fig. 2 isa symmetric function of A, B, and C, 
that is to say, for-this function F(A, B, C)=F(B, A, C) 
Bree Bb A)=F( CAB) =F(C2A 3B) = F(B, 6,4) 
=F(A, C, B). For this particular function, the repre- 
sentation onacube exhibitsa further kind of symmetry, 


F(~A, ~B, ~C) = ~ F(A, B, C) 


hecause reflection in the center of the cube (inversion) 
gives ~F, 


Fig. 2—Representation of F(A, B, C) ona cube. 


The representation as a labelled cube, therefore, 
shows up symmetries in the function. The operations on 
the labelled cube which map it on itself and can give 
tise to symmetry are: 


B a A CUGs 
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These form a group G of order 96. It is the direct prod- 
uct of four cyclic groups of order 2 and the group of 
permutations of A, B, and C of order 6. 

Any form of symmetry of the marked cube is given 
by a subgroup of G under all of whose operations 
ECA. BG) 1s invariant: 

Suppose now that we have considered a particular 
F(A, B, C) and know how to produce it by connecting 
together some two input black boxes. Interchange of 
inputs is a trivial operation and so is negation of inputs 
and outputs since any such negations can be absorbed 
into the black boxes to which they are connected. Acting 
in these ways from the original F(A, B, C), in principle, 
96 different functions such as, for example, 


TRG; A) 


can be formed. Actually, because of symmetry, some of 
these might be the same. If the order of the symmetry 
subgroup is 7, there will in fact be 96/7 different func- 
tions. This reduced set of functions will be called a 
“class” of functions. When one member of a class of 
functions has been realized by some arrangement of 
black boxes, all the other members of the class can be 
realized by trivial alterations to the arrangement. 

The synthesis of switching functions of three varia- 
bles has been treated,! but a restricted set of black boxes 
was employed to correspond with a certain set of elec- 
tronic circuits. The symmetry classes were enumerated, 
but there were more classes than we find here, because 
the symmetry operation 


Ve ater a 


(negation of output) was not employed. Slepian? gives 
a method for calculating the number of symmetry 
classes for switching functions of 2 variables but here 
also only input transformations were employed and the 
output transformation F>~F was not included. 

A second method of representing functions of three 
variables will be used in this paper in order to show all 
256 functions at once. If the variable C is held fixed, a 
function of A and B is obtained. There are two such 
functions 


F(A, B, 0) = $0(A, B) 
and* F(A, By) =“giCAy B); 


which together determine completely the function 
F(A, B, C). Each of the functions ¢ can have sixteen 
different values. Thus all the 256 functions can be repre- 
sented ina 16X16 array. As an example, we give in Fig. 
3 a picture of all the degenerate functions, and all those 


1 Harvard Computation Lab. Staff, “Synthesis of Electronic 
Computing and Control Circuits,” Harvard Univ. Press, Cambridge, 
IMMERSES. lel, ae aR. 

21D. Slepian, “On the number of symmetry types of Boolean 
functions of » variables,” Can. J. Math., vol. 5, pp. 185-193; 1953. 


Fig. 3—The degenerate functions and symmetric functions. De- 
generate functions are shown by a D: symmetric functions are 
shown by the values of A+B+C which make F=1. The func- 
tions (0, 0) and (15, 15) are both degenerate and symmetric. 


which are symmetric under the permutations of (4, B, 
Gy. 

The numbering of the coordinate axes'in Fig. 3, which 
represent values of $9 and ¢; corresponds to the number- 
ing of functions in Table I. 


Functions WuicH CAN BE RraizEp 
By Two BLAcK BoxEs 
Two “black boxes,” each with two inputs, can be 
connected only as in Fig. 4 to make a three-input 


“black box.” 
HH 
eh} 
3 


Fig. 4-—The only arrangement of two “black boxes.” 

To evaluate all the function realizable by this array, 
all permutations of A, B, and C among the three inputs 
must be considered. The easiest arrangement, bearing 
in mind Fig. 3, is to connect C to input 3. We shall con- 
sider this first. 

Suppose that the black box x gives an output G(A, B). 
Then setting C=0, the output is a function of G(A, B). 
This can be 0, 1, G(A, B) or ~G(A, B). 

‘Hence 


go = 0,1, Gor ~G,: 


) 
similarly 


¢; = 0,1, G or ~G. 


can be any of the 16 functions of two 


The 


variables. 


funetion Gc 
There are six cases: 


do = @,; = 0; 1, 2, OF 2a; 
go = 1 dy ='0, 1, 2) OF 14: 
¢, = 0 go = 0, 1, 2, , Or IS; 
o,; = 1 oy = 0, 1, 2)" ** Ore, 
go = $; = 0,1, 2,-++, or 15, 

go = ~ $1 = 0,1, 2,-++, or 15, 


These give 88 functions lying on both diagonals and 
around the border of the 16X16 chart. 
Next consider the case of B applied to the input 3. 


Suppose that the box y gives an output G(x, B) where 


x is the output of black box x. 
If C=0 we have 
= G(0, B), G(1, B), G(A, B), or G(~A, B). 
li C=1, also 
¢: = G(O, B), G1, B), GUA, 


The only cases of interest are: 


B), or G(~A, B),. 


go = G(0, B) ¢1 = G(A, B), 
go = G(I, B) g1 = G(A, B), 
da = GCA, Bd ep GAs BI 
which is already covered, and 
= G(A,B) ¢1= G(~A, B), 


and those obtained by interchanging ¢ and 4. 

It only remains to consider the case of A applied to 
input 3. This is the same as the last case with A and B 
interchanged. Interchange of A and B is a transforma- 
tion of the ¢’s and is easily applied. 

If all these possibilities are enumerated and plotted 
on the 16X16 chart, a total of 152 functions are found 
to be realizable by means of two black boxes. These in- 
clude all the 38 degenerate cases. 

Obviously, the 114 nondegenerate cases can be di- 
vided into symmetry classes because if any function is 
realizable in two black boxes so are its transforms under 
the symmetry operators. We list in Table II the five 
symmetry classes, the order of their symmetry group 
and the number of different functions in the class. A 
typical member of the class, one exhibiting obvious 
symmetry, is shown also. 

The use of the negator symbol is’ unnecessary, but it 
enables the black boxes to be restricted to & and ZA. 

It is known that & and # boxes, with negators are 
sufficient to produce any two-input box. Negators on 
the inputs and outputs of the boxes in Fig. 4 produce 
functions of the same class, so the five arrangements 
shown above are the only possible ones. Negators on the 
connecting link between the boxes in cases a, c, and d 

can be transferred via a # box to the inputs or outputs, 

where they are unnecessary. Table IT, therefore, gives 
another proof that only 114 nondegenerate functions 
are realizable by two boxes. 
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TABLE II 
THe SYMMETRY CLASSES REALIZABLE BY Two “Brack Boxes” 


Symbol Used Order o 


uNncitons 


Typical No. in 
wy : re Symmetry oe 
for Class Member Crone Class 
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Fig. 5—The symmetry classes a, }, c, d, and e. 


These functions are plotted in Fig. 5 on the 16X16 
chart, the class symbol of each function being shown. 


Functions WHicH CAN BE REALIZED BY 
THREE BLAacK Boxes, But Not Two 
Three black boxes can be connected in only two ways 
‘o produce new functions of three variables. These are 
shown in Figs. 6 and 7. 
To enumerate all the possibilities arising from Fig. 6, 


oO 


Fig. 6—The first arrangement of three “black boxes.” 


Fig. 7—The second arrangement of three “black boxes.” 


we start by considering the possible values of functions 
x and y. We are only interested initially in getting rep- 
resentatives from each symmetry class, so negations of 
inputs 1, 2, or 3 or of the output, and interchanges of 
inputs are irrelevant. Therefore, only the possibilities & 
and # for x and y, with a possible negation on one of 
the two branches from input 2, need be considered. The 
case of x and y both & reduces to the case of two boxes. 
The negation is irrelevant if either x or y is #. There- 
fore, only three useful cases arise: 


x a, 
en (~2) & 3 
(aes 2 DES 


ee 2 3. 


These must be considered with all nondegenerate func- 
tions Z, but regarding negation of the output as irrele- 
vant. Functions 1, 2, 4, 6, and 8 can be used. Of the 15 
possibilities, all but 8 are found to be covered by the two 
box representation. 

These eight functions and all their transforms under 
symmetry operations were plotted on the 16 X16 chart. 
It was found that only three were essentially different, 
that is there are three new classes generated by the ar- 
rangement of Fig. 6. 

The new classes, representative members of the 
classes and the order of their symmetry group are listed 
in Table II]. The classes are plotted on the 16X16 chart 
in Fig. 8. 

We still have to consider the arrangement shown in 
Fig. 7. This is relatively simple. 

Suppose C, A, and B are connected to the 1, 2, and 3 
inputs. Let the output of x be 


G(A, B). 


TABLE Tit 


Cie SYMMETRY CLASSES REALIZABLE BY THREE “BLack Boxnrs” 


Symbol Used Typical Shen lhaed No. in 
for Class Member : aca Class 
mOup 
f re 8 
| 
4 24 
g 2 
3 
! 
h 2 2 48 
Total 80 
3 


Fig. 8—The symmetry classes f, g, and h. 


Then the values of the output of y for C=0 and C=1 
are 
vo, V1, = 0,1, Gor ~G. 
lor the box Z we need consider only the possibilities 
L=A&y, Z=A#y. 


When all these are considered, the only cases not cov- 


ered by the same arrangement as Fig. 7 with box Z re 


moved are found to be the ones lke 


@o = A, oi == Aner G 


and 


do = A, $1 = AFG. 


The first of these gives only degenerate functions and 
functions realizable by two boxes. 

The second gives some of the fuuctions plotted in 
Fig. 8, those on the columns (3, G), (5, G), (10, G), and 
(12, G), and similar rows. These are only functions of 
the g and / classes. All the other functions in these rows 
and columns are degenerate or realizable by two boxes, 

We deduce, therefore, that classes g and h can be 
realized by the arrangement of Fig. 7, but not class f. 
It is arguable that class f is more difficult to realize be- 
cause it contains the symmetric functions 4+B+C 
=0, 3 and 4+B+C=1, 2. The latter is the function 
shown in Table III. 

Table IV shows typical members of classes g and h, 
produced by the alternative arrangement of boxes, 
They have been constructed and numbered to give 
functions identical to those in Table IIT. 


FuNcTIONS WHICH REQUIRE MorE THAN 
THREE BLack Boxks 


The number of functions still unrealized is 24. If the 
charts of Figs. 5 and 8 are compared with the symmetric 
functions in Fig. 3, it will be found that the symmetric 
functions 


A- B= C= 0,1 / 
and class j 


AB Ae Cie 2, 3 


have not yet been realized. These are easily seen to be- 
long to a single symmetry class. Also the symmetric 
functions 


A+B+C#=1 
A+B+C=2 

class k, 
Aa be Cte 
AD hel Cai 0e ies 


have not yet been realized, and obviously these belong 
to a single symmetry class, which is different from the 
other. These classes will be called 7 and k, respectively. 

If these classes are plotted on the 16X16 chart, as in 
Fig. 9, it is found that all the remaining 24 squares are 
filled. Therefore, the list of classes is complete. The 
question of how many black boxes are needed to realize 
the members of classes 7 and & is answered by Table V, 
which shows typical members realized by four boxes. 
These arrangements were obtained by trial and error 
because enumeration of all the functions realizable by 
four boxes would be lengthy. 


3 This is short expression for “The function which takes the value 
1 when 4+B-+-C=0 or 3.” 
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REALIZATION OF CLASSES g AND h WITH THE ARRANGEMENT 
or Fic. 7 
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Fig. 9—The symmetry classes 7 and k. 


THE SYMMETRY OF THE CLASSES 


The symmetry, exhibited by the ten symmetry classes 
described above, (and the four symmetry classes of the 
ieyenerate functions) is obviously related to the 32 point 
<toups known to-crystallographers.’ Not all the sym- 
“etry classes have different symmetry. Classes e and h 
have the same symmetry, a single plane of reflection 
“hich corresponds to symmetry in two of the variables. 


es b and & have the same symmetry, a trigonal axis 

‘h three intersecting planes which correspond to 
~Timetry in three variables. It happens that all the re- 
Mining classes have different symmetry, making 12 


” 


‘Hilton, “Mathematical Crystallography, 


“ Oxford University 
Mress, New York, N. Y.; 1903. 


N 
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TABLE V 
THE SYMMETRY CLASSES 


Symbol Used Typical : Wee" No. in 
for Class Member Be Class 
3TOup 
j 12 8 
Ate tC 23 
k 
6 16 


AtB+ C=] 


types of symmetry altogether. By no means all the sorts 
of symmetry which are conceivable actually occur. For 
example there is no completely unsymmetric case. 

Symmetry of a type not considered in the point 
groups occurs in classes a, ¢, g, and j (and two degener- 
ate ones). Because we allowed negation of output as a 
symmetry operator, relations such as 


F(~A, B, C) = ~ F(A, B, ©) 


are regarded as symmetry. This is analogous to a plane 
of symmetry such as 


F(~A, B,C) = F(A, B, C) 


and we can call it “a plane of the negation kind.” There 
are also diagonal axes of the negation kind, as well as 
tetragonal axes of the negation kind which are diagonal 
axes of the usual kind. Inversion of the negation kind 
occurs as a symmetry operator in classes a and j. This 
is simply 


F(~A, ~B, ~C) = ~ F(A, B, ©). 


Classes a@ and j also have hexagonal rotary reflection 
axes of the negation kind which are, as well, trigonal 
axes of the usual kind. All the new types of symmetry 
are shown in stereographic projection in Fig. 10, using 
the notation of Hilton’s book with operators of the 
negation kind shown dotted. 


Class <¢ 


Class 9g 


Fig. 10-——Types of symmetry with operators of the negation kind. 
The outer circle in class g and the dashed lines in class 7 do not 
represent planes of symmetry. By J’ is meant inversion of the 
negation kind. The trigonal axes in @ aud j are also hexagonal 
rotatory reflection axes of the negation kind. 
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ARRANGEMENTS OF BLAack BoxrEs Wuicu GIvE 
ANY FUNCTION OF THREE VARIABLES 


We have demonstrated that, using four black boxes 
only, any of the functions of three variables can be pro- 
duced. To achieve this different interconnections were 
employed for each function, and the least number of 
boxes was used for each function. 

The question now to be studied is whether any single 
arrangement of four black boxes is capable of producing 
(by assigning suitable functions to each box) any of the 
256 functions of three variables. We shall still allow 
permutations of the three inputs. It will be shown that 
there are several such arrangements, and they will be 
enumerated. We shall call these arrangements which 
(with permutation of inputs) can represent any func- 
tions of their inputs “comprehensive” arrangements. 


THE NONTRIVIAL ARRANGEMENTS OF 
Four Brack Boxes 


It is easy to enumerate the possible arrangements of 
four boxes. Arrangements in which the output of a box 
influences its tnput (i.e., feedback) are not allowed. 
Some arrangements have a section containing more 
than one box which, when isolated from the rest is seen 
to have only two inputs and one output. This section 
can be reduced to a single box, so in this case, the ar- 
rangement is equivalent to one with less than four boxes 
and need not be considered. 
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Fig. 11—The arrangements of four “black boxes.” 


When arrangements of these sorts have been elimi- 
nated, the nine arrangements shown in Fig. 11 are 
found. The details of the enumeration will not be given. 
(The asterisks are explained later.) 


THE FUNCTION OF CLASSES a TO h 


The functions of classes a to h can be realized by less 
than four black boxes. One way of realizing them with 
the arrangements of Fig. 11 is to allow one or more 
boxes to become trivial functions, so that the arrange- 
ments reduce to 3-box and 2-box arrangements. 

There are two arrangements of:three boxes, the ones 
shown in Figs. 6 and 7, and the arrangement of Fig. 6 
covers all the classes a to h. The arrangement of Fig. 7 
covers classes a to h except for class f. 

It can easily be seen that, by making the appropriate 
box degenerate, the arrangements 1 to 7 of Fig. 11 can 
be made into the arrangement of Fig. 6. For these, 
therefore, only classes 7 and k remain to be investigated. 

The arrangements 8 and 9 cannot degenerate to that 
of Fig. 6, but they do go into Fig. 7. or these arrange- 
ments, therefore, classes f, 7, and k remain to be in- 
vestigated. 

The classes to be investigated, f, 7, and k, all contain 
symmetrical members: 
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F(AB,c)= 6{x (A,B), 8,cf 


Fig. 12—A way of breaking down a complex arrangement. 


for, A a Boe Co = 
fory, A+"B+C = 
fork, A+ Ba C = 


OF 3 sores 2 
OF ior 25-3; 
i Ole Ore t), 5 0b 0 uleko 


us makes 1 ge ecause permutations 
This makes investigation easier becat perm tati 
of inputs need not be considered. 


MeEtTHODS FOR DETACHING ONE Box 
FROM AN ARRANGEMENT 


Suppose that an arrangement of black boxes with 
three inputs and one output can be represented by 
Big. 12. 

Suppose also that F(A, B, C) has this property: 
whatever the value of B, there is some value of C for 
which F depends on A. Then whatever the value of B, 
X depends on A. There are only the following possi- 


bilities: 
\ — A 
= (4S 


Table VI shows that functions of classes f, 7, and k have 
this property. 


B) and their negations. 


TABLE VI 
Class | Function B=) B= 
if A+B+C=0, 3 ~(AVC) LSC 
i A+B+C=2, 3 Bl eC AVC 
‘k A+B+C=1 A#C A 


The arrangements 1, 2, 3, 5, 6, 7, and 8 of Fig. 11 can 
be represented by Fig. 12, (arrangement 3 in two ways). 
[herefore, in representing functions of classes f, 7, and 
k, the starred boxes on the upper levels of these arrange- 
ments can be assumed to be The de 
Case x= 


=4 boxes. 


5 generate 
A is of no interest. 
Suppose now that an arrangement of black boxes can 


be represented in the form of Fig. 13. 
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F(ABC)=G{ax(asc)} 


Fig. 13—A further way of breaking down a complex arrangement. 


Suppose also that F(A, B, C) has this property: what- 
ever the value of A, F(A, B, C) can take both'the values 
0 and 1, depending on B and C. Then whatever the 
value of 4 ,G(A, X) depends on X. Again the only useful 
case is GOA, A)=(A#X). 

The functions shown in Table VI have this property 
and the arrangements numbers 2, 3, 8, and 9 of Fig. 11 
can be represented by Fig. 13. Therefore, in representing 
functions of classes f, 7, and k, the starred boxes in the 
lower levels of these arrangements can be assumed to 
be # boxes. 

Having found the nature of these boxes, the problem 
of whether a given (A, B, C) can be represented by the 
whole array can be reduced to the problem of represent- 
ing a different function by the smaller array of un- 
starred boxes. 

Take the case of Fig. 12. Since X =(A ¥B), it follows 
that. 4=(X SB). HencewG 4, 8, OS 4 a) 
= F(X 4B, B, C) is represented by the lower part. 

In the case of Fig. 13, since F=(A 4X), it follows 
that X =(A#F). Hence X(A, B, C)=A¥F(A, B, C). 
So F is represented by the whole array if A ¥F is repre- 
sented by the upper part. 

In this way the problem is reduced in complexity, but 
it must be borne in mind that in the modified functions 
we are not allowed to interchange the inputs. This would 
alter the array. 

Only one of the nine arrangements is unaffected by 
this method of reduction, namely arrangement 4. Fortu- 
nately, this has been found to be comprehensive by 
other methods. 


DEMONSTRATION OF THE COMPREHENSIVE 
ARRANGEMENTS 


The following results concerning arrangements 1 to 9 
were discovered by the method of reduction described 
in the previous section. It will be remembered that all 
the arrangements can, by degeneration to a three box 
arrangement, produce functions of classes a to e, g, and 
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AtBt+C=2,3 At+B+C ¥2 


Arrangement 2 
A 


A+B+C=23 


Arrangement 4 
B Gc 


AtBtC+23 AtBtC= 2 


Fig. 14—The “comprehensive” arrangements. 


h, and that all except arrangements 8 and 9 can produce 
functions of class f. 

Arrangements 1, 2, and 4 are comprehensive, and 
their representation of classes 7 and k are shown in 
Fig. 14. The other arrangements fail in the following 
classes. 


Arrangement Classes 
3 k 
5 j 
6 Gels 
7 j 
8 if Ws R 
9 fj 


As an example of the method of reduction, consider 
arrangement 9 as a possible method of realizing the func- 
tion of class 7: 


lI 
to 


Art. B= C 


This is an example of the type of Fig. 13, where -1 is 
the tnput marked 3 in Fig. 11. 

If A=0, F=B&C, and if A =1, F=B\C. Both can 
take the values 0 and 1, so the starred box can be taken 
to be a # box. Therefore, the unstarred boxes represent 


the function 
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ie 


A+B+C=23 At+BtC = 23 


Fig. 15—Alternative realizations of a function. 


X= (A # Fh, 


for which if d =0, X =B&C:andif A = 1,X=~(BVC), 
This is seen, from Table VI, to be a function of class f. 
hence not realizable by the three boxes (see Table IV). 

All the arrangements except 1, 2, and 4 give similar 
results. 

There are, therefore, three comprehensive arrange- 
ments of four black boxes, that is, arrangements which 
can generate any function of three variables by a suit- 
able choice of boxes. The choice of boxes is not unique. 
For example, the two choices in Fig. 15 give the same 
output. 


COMPREHENSIVE ARRANGEMENTS WITHOUT ALLOWING 
PERMUTATIONS OF INPUTS 


The final problem to be discussed is: which ar- 
rangements of four black boxes are capable of producing 
all 256 functions of three variables without allowing 
permutations of inputs? Such an arrangement must 
evidently be comprehensive. 

Consider the function F(A, B, C)=(A&B)#C. We 
assert that the arrangements of Fig. 16 cannot produce 
this function. For, by applying the argument of an 
earlier section, we find that the boxes marked with a 
star in each case must be either # or degenerate. The 
degenerate cases, shown in Fig. 17, are easily shown to 
be impossible by further reduction. The nondegenerate 
case of Fig. 16(a) gives a new function for the unstarred 
boxes, 


G(A, B, X) = F(A, B, X # A) 
eS (A & ~ B) A D6. 


Fig. 17(a) must represent this function with X replac- 
ing C. Further reduction shows this to be impossible. 

The nondegenerate case of Fig. 16(b), gives the ar- 
rangement of Fig. 18, where 


V(A, X, O'= AZ BACH YO 
= (Suite CAA aoe 


Thisisa function of class g and requires three black boxes. 
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A c A e B 
‘ g 
B 
x 
a b 
Fig. 16 
A c B B A c 
b 
The only remaining comprehensive arrangement, 
which is number 4, is capable of giving any function of 


three variables without permutation of inputs. For 
suppose an arbitrary F(A, B, C) is specified by giving 


oA, B) a EAS Bb, 0) 


Fig. 17 


and a’ 


$:(A, B) = F(A, B, 1) 


Davies: Switching Functions of Three Variables 


A C 
x 

\f 

Fig. 18 


Fig. 19—A completely general arrangement with four boxes. 


and let ¥(4, B)=¢.(A, B)#¥d¢i(A, B). Then the ar- 
rangements of Fig. 19 gives F(A, B, C). 
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Summary—This paper begins with Mealy’s model of a sequen- 
tial machine and introduces a ‘connection matrix” which describes 
the machine completely. The ‘‘equivalence’”’ of states of such a 
machine may be analyzed systematically by an iterative technique, 
the validity of which is rigorously established. Once equivalence is 
completely analyzed, it is a simple matter to write the connection 
matrix for the simplest equivalent machine. The process is not diffi- 
cult to execute, even in complex cases, and could be programmed 
for a computer. 


I. SEQUENTIAL MACHINES 


NE may designate as a “sequential machine” any 

device which produces prescribed sequences of 

outputs In response to given sequences of inputs. 

A number of mathematical models of such machines 

have been proposed.!~® Mealy’s model of a sequential 

machine! assumes that the machine may be described in 
terms of 


1) a finite number of states s1, 5», ba So 

2) a finite number of distinct zzputs symbolized by 
Bu, Xs; oi Mg eS 

3) a finite number of distinct oulputs sy ymbolized by 
41; V2, ** * yp. Itis assumed that the present out- 
put and the next state are uniquely determined by 
the present state and the present input. 

Each such machine MJ may now be represented by a 
weighted, directed graph or net,®’ called the state dia- 
gram of AZ. The set of vertices v1, v2, -, v, of the net 
corresponds in 1-1 fashion to the set of states si, 52, °--, 
s, of M and the s;’s are the weights associated with these 
vertices. For cach input which effects a transition of the 
machine from state s; to state s;, there is a directed 
branch from the vertex v; to the vertex v; and we assign 
this branch as its weight, the ordered pair consisting of 
the input symbol in question and the symbol for the 
corresponding output. When there is more than one 
such branch, say, k branches from a vertex v; to a vertex 
v;, we simplify matters by employing only one branch 
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and assigning it a weight (x,,, yr.) V(r Vr.) V 

NY as (where \/ means “or”) which defines the 
union of the weights of all the original branches from 
U; LO U5: 

A sequential machine satisfying the above require- 
ments will be referred to simply as a machine in the rest 
of this paper since no other types of sequential machines 
are discussed. 

It is often convenient to identify the states simply as 
peek -, m and to use unsubscripted symbols for in- 
puts and outputs. We illustrate in Fig. 1 by giving a 
state diagram of a binary pulse-divider. The states are 
1 and 2, the input symbols are 0 and 1, and the output 
symbols are 0 and 1 also. Here the input sequence, 


1, 1, 1, 


starting with the machine in state 1, produces the out- 
put sequence 


0, 1, 0, 1, 


However, no matter what the input sequence may be, if 
the machine starts in state 1, every second 1 in the input 
sequence will result in an output 1, all other outputs 
being 0’s. 


(1,0) 


(0,0) (1,1) (0,0) 


Fig. 1—State diagram of a binary pulse-divider. 


II. THE CONNECTION MATRIX 


The basic tool of analysis in what follows is the con- 
nection matrix of the net associated with the machine. 
We recall its definition from Hohn, et al.® 


[coshasen 


where ¢,; is the above defined union (\/) of the branch 
weights on all branches from vertex 7 to vertex j. If no 
branches connect vertex 7 to vertex j, then ¢i7=0.. No 
input symbol will appear more than once in any row of 
C since the next state is uniquely determined by the 
present state and the present input, 7.e., 
deterministic machine. 

The connection matrix of the binary pulse-divider 
illustrated in Fig. 1 is 


sf Oe UL,0) 
Co leaehs el 


since Af isa 
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III. Inrut-Output POLYNOMIALS 


In what follows, we shall frequently find it necessary 
to employ certain VV polynomials in the input-output 
pairs associated with a machine. We now define these 
polynomials and list their basic properties. 

In terms of the notation of Section 1, the possible 
input-output pairs (xj, y;) are mp in number. Let us 
denote them by Pi, Ps, ++ +, Py. From the P’s we then 
form input-output polynomials with the aid of finitely 
many operations of multiplication (indicating succes- 
sion) denoted in the usual way, and union (the logical 
“or”) denoted by “V”. Concerning the set S of such 
polynomials, we now make the following postulates, in 
which “=” denotes identity and f, g, h are arbitrary 
members of S: 


eines, t wp Delong to75, 
2) fVg belongs to S, 

3) FV eE=EVy, 

AV iV I =f; 

5) (Ve) Vh=fV(gVh), 

6) fg belongs to S. 


There exists in S an input-output polynomial “0” (not 
to be confused with a possible input symbol 0) such that 


7) fVO=f and 
8) f-0=0-f=0, é 
9) If f40, gX0, f¥zg, then in general, fe gf, 
10) f(gh) = (fg)h, 
11) (Vg)h=fh\V eh and 
PEEP LEN) Fico 0) 


A polynomial P;,P;,, --~-, P;, will be called a term. 
We define an expanded polynomial of S to be a term ora 
_ union of terms. Every polynomial of S may be expressed 
as an expanded polynomial with the aid of the above 
postulates. If each term of an expanded polynomial f of 
Sis a product of r factors, the polynomial will be said to 
be homogeneous of degree r. The polynomial 0 is assigned 
no degree. 

The above postulates allow us:also to define a useful 
form of matrix multiplication. Consider two matrices 
‘laxg and Bgx, whose entries are members of S. We de- 
fine 4 B in the usual way except that “\/” replaces “+.” 


AIR = [axib1; W/ aj2bo; ye Naas Ws asda; lax 


or, more compactly, 


B 
AB == | V ous | 
k=1 


Now let C be the connection matrix of a machine M. 
Vhen, from the definition of C it follows that 


B 
Ca L V cats | 
et 


5 a matrix whose 7 entry gives all imput-output se- 
quences of length 2 which take J/ from state 7 to state 7 
Via some intermediate state k (not necessarily distinct 
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from z or j). More generally we have 


: ee | 


the 17 entry of which gives all tput-output sequences of 
length r which take M from state z to state j via r—1 
intermediate states. 

The nonzero entries of C’ will be homogeneous poly- 
nomials of degree r. Moreover, since the next state is 
uniquely determined by the present state and the pres- 
ent input, no two entries in the same row of C* can have 
terms involving the same input sequence. 

It is also useful to define the union of two matrices 
Aaxg and Baxg as the matrix A\/B defined by 


AV B= [(ay V bis) axe. 
The union has the properties 
ASB. = Bue Ax 
CANEBY NEE SAD Be: 


of which we shall make use. 


n 
Cr = | V Cik Chik, °° 
ky.ko, CH Pos a 


I 


IV. Input RESTRICTIONS 


It is possible that in the state diagram there will be 
vertices such that not every input symbol is associated 
with a branch leaving that vertex. In this case, when 
the machine is in the state corresponding to this vertex, 
its output and next state are not defined for every input 
and we say that the machine has input restrictions. If 
there are no input restrictions, the machine must be 
capable of responding to any input no matter in which 
state it may be so that at each vertex, each input symbol 
will appear on some branch leaving that vertex. 

To illustrate, suppose there are three input symbols 
a, b, cand that the machine contains a state 7 such as is 
shown in Fig. 2. Here the machine would not be able 
to accept an input a if it were in state 7. 


(4, 8) 


(Gee 
(c,%) 


Fig. 2—Example of an input restriction associated with state 7. 


The presence of input restrictions is related in practice 
to the problem of designing a machine which responds 
in a prescribed way to a limited number of ¢uput se- 
quences (t.e., sequences of input symbols). Then the re- 
strictions which define the given set of imput sequences 
are reflected in the resulting state diagram. This prob- 
lem will be treated in detail in a later paper of this series. 
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At the moment, our concern is the analysis of a given 
state diagram. For this purpose the next definition is 
useful. 

Any input sequence which, starting with the machine 
in state 2, violates no input restriction of state 7 or any 
subsequent state, is called an allowable input sequence 
for state 1. 

We have noted in the last section that the zj entry of 
C" lists all sequences of r input-output pairs which take 
the machine from state 7 to state j7. Neglecting the out- 
put symbols, we have then in the 7th row simply the al- 
lowable input sequences of length r for state 7. (If only 
input sequences are of interest, we omit the output sym- 
bols entirely from C and compute powers of this simpler 
matrix. When this is done, one must distinguish at 
times between an input symbol “0” and the 0 poly- 
nomial mentioned above.) 

If no input restrictions are present, every input se- 
quence may be presented to A/ regardless of its initial 
state and hence every input sequence of length r must 
appear somewhere in each row of C’. 


V. SYNCHRONOUS AND ASYNCHRONOUS MACHINES 


The notion of input restrictions enables us to draw an 
important distinction. The basic characteristic of syn- 
chronous machines is that with the aid of clocking de- 
vices they are enabled to distinguish between succes- 
sive, identical inputs. This suggests the following defi- 
nitions. eee 

If the allowable sequences for the various states of a 
given machine JZ include at least one sequence in which 
an input symbol follows itself, J7 will be called essen- 
tially synchronous. If nosuch sequence is included, 7.e., if 
an input symbol never follows itself in an allowable in- 
put sequence, M will be called essentially asynchronous. 

When we say “essentially asynchronous,” we do not 
mean to imply that the machine may not be realized 
with synchronous circuitry. We mean only that the ma- 
chine is never required to distinguish between succes- 
sive, identical inputs. Thus the input sequence in Fig. 3 
would be recognized by an essentially asynchronous 
machine as 0 1 O whereas an essentially synchronous 
machine might recognize it as 0 1 1 0 0 0, for example. 


1 


Pi Meer, Si epee eS 


tp tH iG 


0 


Fig. 3—Input sequence. 


VI. EquivaALENT STATES AND 
EQUIVALENT MACHINES 


In problems of analysis and synthesis of sequential 
machines, we often wish to determine whether or not a 
given state diagram can be simplttied, that is, whether or 
not it can be replaced by one having fewer states with- 
out altering any of the relations between input and out- 


TRANSACTIONS ON ELECTRONIC COMPUTERS 


Deceinber 


(0,0) a =) 0 


Fig. 4—State diagram and associated connection matrix of a machine 
having two inputs and two outputs in which states 2, 3 and 4 
are equivalent. 


(1, 1) 


put sequences which characterize the machine. 

The notions of equivalent states and equivalent ma- 
chines, which ultimately go back to Moore? and Huff- 
man,* are essential to the solution of this problem. We 
therefore make the following definitions. 

Definition 1: Two states, s; of a machine M and ¢; of a 
machine N, are called equivalent states if, and only if, 
they have the same set of allowable input sequences 
and, starting with VW in s;and N in ¢;, the two machines 
will yield identical output sequences on being presented 
with the same allowable input sequence. In this defini- 
tion Jf and V may of course be the same machine. 

Definition 2: Two machines M and N are called equiva- 
lent machines if and only if for each state s; of M there 
exists at least one equivalent state ¢; of N and for each 
state t; of N there exists at least one equivalent state 
s70f Af: 

The relations of “equivalence” here defined are true 
equivalence relations. It is in fact easy to check that the 
determinative reflexive, symmetric, and transitive prop- 
erties hold in each case. 

As an example, consider the machine with inputs 0, 1 
and outputs 0, 1 whose state diagram and connection 
matrix are shown in Fig. 4. Here every input sequence is 
allowable for every state. 

Now for two states 7 and j to be equivalent, it is nec- 
essary that the same input-output pairs appear in rows 
t and j of C. (That this condition is not sufficient is 
illustrated in the example in Section I. In the above 
example, it is therefore possible that states 2, 3, and 4 
are equivalent. The input 1 takes each of states 2,.3, 4 


(A,B), (C,¥) 


(A,B) 


(C,¥) 
(2B 0 


0 (C,é) 0 Gx 


(6,4) 0 
(6,0) 0 
), 69) 
0 Ge) (8) | 


Fig. 5—Machine with three inputs a, b, ¢ and three outputs 
a, 8, y of which € is one. 


(Ca) 0 


into state 1, yielding the output 1 in each case. The 
input 0 simply permutes these states among themselves, 
now yielding the output 0 in each case. By inspection of 
the connection matrix we therefore see that any input 
sequence will result in the ‘same output sequence 
whether we start with the machine in state 2, 3, or 4. We 
now replace the class of equivalent states 2, 3, 4 by a 
single state 2’ which is taken by input 1 into a state 1’ 
and by input 0 into state 2’, the outputs being 1 and 0, 
respectively. Input 0 takes 1’ into 1’; input 1 takes 1’ 
into 2’. This yields a state diagram essentially the same 
as that previously shown in Fig. 1. The connection 
matrix of this machine, namely 


be 0) ane 
(1,1) (0,0) 
is readily obtained by replacing each submatrix of the 
partitioned matrix of Fig. 4 by the union of all its en- 
tries. This observation will presently be reduced to a 
formal rule. 

It remains to point out that the original machine and 


the new machine are equivalent, the correspondence of 
states being 


1< 1’ 


and- (2;3,4) <2’. 


In fact, for state 1(1’) the input 0 yields output 0 and 
a transition to state 1(1’) but input 1 yields output 0 
ind a transition to state 2(2’). For states 2, 3, 4(2’) in- 
put 0 yields output 0 and a transition to state 3, 4, 2(2’) 
while input 1 yields output 1 and a transition to state 
[(1’). Thus, starting with corresponding initial states, 
both machines will produce identical output sequences 
lor a given input sequence, and hence are equivalent. 
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We now give two additional examples to illustrate 
certain aspects of equivalence before attacking the prob- 
lem formally. Consider first the state diagram and as- 
sociated connection matrix shown in Fig. 5. The ma- 
chine is assumed to have 3 inputs a, 6, c and 3 outputs 
a, B, y of which € is one. There are no input restrictions. 

First we consider states 1 and 2, noting that rows 1 
and 2 contain the same input-output pairs. Hence 
equivalence of these two states is possible. Input @ takes 
both states into state 1 whereas 6 takes both into state 3. 
Input ¢ on the other hand merely permutes states 1 and 
2. From these several observations we see that, starting 
with the machine in state 1 or state 2, a given input 
sequence will result in the same output sequence in 
either case. Hence 1 and 2 are equivalent states. 

Now examine states 3 and 4 under the assumption 
that =a in which case rows 3 and 4 contain identical 
input-output pairs. The input ¢ then takes 3 and 4 into 
states already known to be equivalent whereas inputs a 
and b take 3 and 4 into state 4. It again follows that 
states 3 and 4 are equivalent. On the other hand, no 
further equivalences are possible because of different 
input-output pairs in the rows. 

We now replace each submatrix of the partitioned 
C matrix by the union of all its entries, thus obtaining 
the 2X2 matrix and associated state diagram shown 
in Fig. 6. 


ADDITIONAL EXAMPLES 


(6,X) 


(2,B), (C,9) 
(Ca) 


(g,0) | 1 
Ce), (Banya 


Fig. 6—The machine given in Fig. 5 when reduced under 
equivalence of states assuming that =a. 
As in the example of the preceding section, one can 
verify that this machine is equivalent to the originally 
given one, the correspondences between states being 


(4D )< el’ (3, 4) <> 2’, 


and 


Finally, suppose that a, say £=8. Then states 3 
and 4 cannot be equivalent, because of the different 
outputs associated with input c. In this case we par- 
tition C further: 


(a, 8) (c, ¥) | (b, a) | 0 

ae (a, B) v. (, a8 bait (bo oy) ity a 0 
Sr SAS (c, 8) pe as (a, 7) Vv (b, y) 
(c, a) Uo ahr Ors (a, os) Wi (b, Ne 


weit 


tor eae a ALZI TER 


(2,8), (C, ¥) Mey: @,%), (4,¥) 
2 
ea! 
it 2 Se Nive 
@pvicy (4% oO | 


(C,p) 0 Gaver} 2 
(C,&) 0 


Fig. 7—The machine given in Fig. 5 when reduced under 
equivalence of states assuming that §=£. 
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so that the rows of each submatrix correspond to equiva- 
lent states and similarly for columns. Again, we replace 
each submatrix by the union of all its entries and thus 
obtain the connection matrix of a demonstrably equiva- 
lent machine with fewer states, as shown in Fig. 7. 

The type of equivalence exemplified by states 1 and 2 
above, and by 3 and 4 when £=a, is readily generalized. 
Suppose 7}, - ++, % are states such that the same set of 
input-output pairs appears in each of the corresponding 
rows of C. If each input «;, of the set of inputs ae oe IC ea 
x;, takes all of the states 7),.-- - , 7, into a fixed state 1s, 
or into states equivalent to 7, and if the remaining inputs 
to the machine simply map the set 7%), - + + , im onto or 
into itself, then the states of this set are equivalent. 
More generally, there may be a number of disjunct sets 
of inputs as the set x;,, + - 

These conclusions will be justified by the theorems to 
be proved in Section VIII. A corresponding partitioning 
of the connection matrix, analogous to that used in the 
above example, may be effected by the methods of Sec- 
tion IX. 

As a second example, consider the state diagram and 
matrix shown in Fig. 8. The inputs are a, b; the outputs 
are a, B, y. The input sequences are restricted since the 
machine will not respond to either input once it is in 
state 5. (We call state S$ a terminal state of the machine.) 

State 5 is obviously equivalent to no other state here. 
However, states 1 and 2 are permuted by input 6 while 
input @ takes them into states 3 and 4, respectively. 
Since the input-output pairs in rows 1 and 2 are the 
same, we conclude from the generalization of the previ- 
ous example that 1 and 2 will be equivalent if 3 and 4 
are. Now input @ takes 3 and 4 into states 2 and 1 re- 
spectively while input b takes them into a common next 


nes Nips 


state, namely state 5. Thus we would conclude that 3 
and 4 are equivalent if 1 and 2 are! However, examina- 
tion of the possible input-output sequences reveals that 
1 and 2 are indeed equivalent, as are 3 and 4. 


Fabshs sd 2A LL 


CRAs UE ffi 
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Fig. 8—Machine having a terminal state in which there exist 
equivalent states. 


As before, we now replace each submatrix by the 
union of its entries and obtain the connection matrix of 
a demonstrably equivalent, 3-state machine, as shown 
in Fig. 9, opposite. 

The type of equivalence illustrated here may also be 
generalized. Suppose that the rows of C corresponding 
to a set of states S; all contain the same input-output 
symbols. Similarly for S, Sy, + - - , S,, Where these are 
disjoint sets of states. If each input x; to the machine 
takes all the states of each S; into the same state, or into 
equivalent states, or into states all contained in some 
S;, then each S; contains only equivalent states. Again, 
this type of equivalence is identified by partitioning the 
connection matrix in a manner analogous to that illus- 
trated in the above example. The correctness of this 
generalization will appear as a consequence of the re- 
sults in the next two sections. 

We are now ready to examine the problem of equiva- 
lence in general. 


VIII. Toe FunpamMentAL THEOREMS 


We begin with some definitions. A matrix A whose 
elements are input-output polynomials is called an 
r matrix if it has the following properties: 

1) all nonzero entries of A are homogeneous and 
of a common degree +, 

2) all nonzero terms in each row of 4 have distinet 
input sequences, and, 

3) all-nonzero terms which appear in any given 
row also appear in every other row. (The submatrices 
in the examples of the last section are 1 matrices.) 


Let the columns of a matrix 4 of order » be parti- 
tioned into g sets of Ni, No, +++, NM, columns where 
Sey eras . 
>, N;=n. If now the rows are also partitioned into 
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Fig. 9—State diagram and associated connection matrix of a 
reduced machine equivalent to the one given in Fig. 8. 
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sets of Ny, No, - - +, N, rows, such that the ordered sets 
of N; column numbers are identical with the ordered 
sets of NV; row numbers, fori=1, -- -,g, the matrix A 
is said to be symmetrically partitioned. (Examples ap- 
pear in Section VIJ1.) Ifa symmetrical partitioning is such 
that Ni=N2= = N,=1,s0 that g=n, then we call 
it the trivial partitioning. We now have: 


Lemma 1: 


The union of two r matrices-A and B of the same order 
1s again an r matrix whenever the terms of corresponding 
rows of A and B contain no identical input sequences. It 
it easy to verify that properties 1), 2), and 3) all still 
hold for the union under the stated hypotheses. 


Lemma 2: 


The product AB of anr matrix A and ans matrix B is 
an r-+s matrix. Property 1) is immediate. To prove 
property 2), it will suffice to show that aix,bi,;, and 
Qix554,;, contain distinct input sequences whenever 
kiko or jiMje or both. Indeed, if kiko, then by hy- 
pothesis the terms of a%, and ai, contain no common 
input sequences. Hence, because of the noncommuta- 
tive nature of multiplication and the homogeneity of 
the entries of A and B, it follows that no two terms of 
Qjrbe,j3, and @ix,be,j3, Can contain common input se- 
quences. If now ki =» but 7;A%j2, the desired conclusion 
follows from the fact that no two terms of 0y,;, and 
b,,;, contain the same input sequence. 

To prove property 3), consider any term ¢ appearing 
in the 7th row of the product AB. Suppose ¢ is con- 
tained in a;,r,b2,;,. Phen the factor of ¢ provided by ai,x, 
also appears in a;,x, for any given 72.747, and for some kp. 
lhe reason is that A is an 7 matrix and every entry in 
row 2 appears somewhere in row ky. Similarly, the factor 
provided by bx,;, also appears in d;,;, for some jz. Hence, 
the term ¢ also appears in the product @i,4,Di3,, 7.€., 
property 3) is established. 
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We can now prove the first fundamental result: 


Theorem 1: Let the submatrices of a symmetric par- 
titioning of a given connection matrix C be denoted by 
Ci;,7,7=1, 2, +--+, ¢. Let the submatrices of an identi- 
cal partitioning of C’, where r is any positive integer, be 
denoted by Cy, 4, 7=1,.2, 99s ~ Gal Weng cach C,. 
1s a 1 matrix, each Cy; ts an r matrix. 


The proof is by induction on r. The theorem is clearly 
true for r=1. Assume that it holds for r=k. Examining 
the matrix C**!, partitioned identically to C and C*, we 
have 


CHI = C-Ch = [ V CC® | 
St 


By Lemma 2, each product C;;,C,;” is a (k-+1) matrix. 
Since C is a connection matrix, the 1 matrices C;,, and 
Cis, have no input symbols in common. Hence no two 
products C;;,C.,;@ and C;.,C;,; can have an input se- 
quence in common. Therefore, by Lemma 1, the union 


q 
V CC. 
=1 


& 


is also a (k-+1) matrix and the theorem is proved. 
We next use this result to characterize equivalent 
states in the connection matrix: 


Theorem 2: If a connection matrix admits a sym- 
metrical partitioning in which each submatrix is a 
1 matrix, then the states contained in each set of the cor- 
responding partitioning of the states are equivalent. 


As we have seen in Section IV, each allowable but 
otherwise arbitrary tmput sequence of length r which 
may be applied with the machine beginning in state 7 
must appear in some term in the zth row of C”. If, on the 
other hand, C” is partitioned identically to C, then each 
of the resulting submatrices is an 7 matrix, by Theorem 
1. Therefore, the output sequences corresponding to al- 
lowable but otherwise arbitrary input sequences as- 
sociated with states in a given partition are identical, 
independently of which state of the set the machine is 
in initially. 


IX. AN ALGORITHM FOR DETERMINING ALL 
SETS OF EQUIVALENT STATES 


Consider the connection matrix C of any machine MV. 
The following iterative procedure will determine all scts 
of equivalent states of AZ. 

1) Partition the states of M into disjunct sets 
Si, So, > + +, S_, of maximal order such that the rows of 
C associated with the states in each S; form a 1 matrix. 
If the resulting partitioning is the trivial partitioning 
(into sets of order 1) then there are no equivalent states 
in M. 

2) If the previous partitioning is nontrivial, reorder 
ips 


if necessary, and partition the corresponding C sym- 


the states according to the partitioning S;, Ss, ++: 


a CtiCCt(itttt((#(##(#é#éw44... eee 
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metrically with respect to the S;, indicating the par- 
titioning by C=[C,;]. If the C; are all 1 matrices, then 
the process terminates at this point because the states in 
each S; are equivalent by Theorem 2. 

3) I{ the C;; are not all 1 matrices, then partition the 

states associated with cach row of submatrices Cu, Cx, 
try tre, FS, 2s > > <4-¢ Inte -disiwact sets Say Sa, 
+, Sia, of maximal order such that for each S;; each 
of the associated submatrices in Ca, Cx, +++, Cyq is a 
1 matrix. If this new partitioning is trivial for all rows 
of submatrices Cy, Cx, +--+, Cig, then there are no 
equivalent states. 

4) If the previous partitioning of the states is non- 
trivial, repartition the connection matrix with respect 
to the sets Sy, Sx, + + -, Sia,, again reordering the states 
if necessary so as to have the rows of states in the same 
set Sj; adjacent to each other in the matrix. Denote the 
new partitioning by 


C rs kon a [D:**]; OF = [D3 **], 


where the indices k and ¢ are used to denote the sub- 
matrices in the partitioning of each submatrix C,;. 

If now the D;,*' are all 1 matrices, the process termi- 
nates at this point, for by Theorem 2 the states in each 
Si; are equivalent. Otherwise we repeat steps 3) and 4). 

We have thus outlined a two-step process: 1) parti- 
tioning the set of states into maximal, disjunct sets such 
that the submatrices formed from the corresponding 
sets of rows in the submatrices of the previous partition- 
ing are 1 matrices, and 2) repartitioning the connection 
matrix with respect to this new partitioning of states. 
These two steps are repeated until a partitioning of the 
connection matrix is obtained in which all submatrices 
are 1 matrices. This is a finite process since the number 
of states is finite. Once such a partitioning is obtained, 
then by Theorem 2 we see that all states in each subset 
of states of the final partitioning are equivalent. If the 
subsets of states ultimately obtained are all of order 1, 
then there are no equivalent states. 

At each stage, the partitioning of the states is 
uniquely defined by the requirements listed above. Thus 
the end result must also be unique. Indeed, what we 
have done is simply to divide the states of the machine 
into equivalence classes. 


X. THE REDUCED MACHINE 


After completing the process outlined in the last sec- 
tion, we define a new machine by replacing, in effect, 
the states of each equivalence class by a single state 
representing that class. To account for the transitions 
of all the states of one class to states of other classes 
under the various inputs to 17, we replace each sub- 
matrix of the final partitioned form of C by the union of 
all its entries. For simplicity of notation denote the final 
partitioned form of C by 


C= [Cx] 1, J = 4 2, 2) eds 
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the corresponding sets of equivalent states being .S4, 
Se Sq Then if d,; is the union of all the entries of 


egies yO 


Ci;, the gXq matrix 
Di= [d;;] 


is the connection matrix of a machine with g states, no 
states of which are equivalent. That D is indeed a con- 
nection matrix is seen by noting that because C is a con- | 
nection matrix and the Cj;’s are 1 matrices, no two 
matrices C;;, and Cj, contain any common input sym-_ 
bols. Therefore di;, and d;;, will contain no common in- | 
put symbol, which is the only condition to be satisfied. 

The machine NV defined by D is called the reduced 
form of M*°) It remains to show that M and W are 
equivalent machines. 

To prove this, for each positive integer 7 we partition 
C” according to the sets of equivalent states of 1, and 
let Cy denote any of the resulting submatrices. The 
entries of C,; define the allowable input sequences of 
length r which take the states of set S; into those of set 
S; and also give the associated output sequences. 

Now the entry dj; of D is the union of the entries of 
C,; and hence defines the set of all inputs which take 
some member of S; into some member of Sj. On the other 
hand, d;;, the 7j entry of D’, is given by 


r 


dp = V Gig Geyky °° * Fey. 


ky ko, 8 * Kp aye 


Any term in the expansion of this product represents an 
input-output sequence taking some state of S; into some 
state of S;. Hence, this term appears in C;;, and be- 
cause Cy; is an r matrix, this term appears in each row 
of Ci. Conversely, since each state of I appears in 
some S;, every term appearing in C,;,™ has its factors 
appearing in some produce of the form dix,dijny « 
d,,,;and hence is a term of d,;™. 

From these observations we see that each of the 
states of MZ in S; has the same allowable input sequences 
and the same associated output sequences as does state 
7 of V. Thus all of these states are equivalent (by Defi- 
nition 1, Section VI) to the state 7 of N,d=1,2,---, q 
and vice versa. Hence, by Definition 2 of Section VI, the 
machines M and N are equivalent. Since the machine V 
is uniquely defined, we may now summarize our ob- 
servations thus: 


Theorem 3: Given a machine M, there exists a corre- 
sponding equivalent machine N with minimum number 
of states which is unique to within a permutation of the 
States. 


This theorem is the same as that of Mealy! which is 
an adaptation of a result of Moore.!! 


S Mealy, op. cit., p. 1054. 
* Moore, op: czt., p. 143. 
10 Mealy, op. cit., p. 1054. 
" Moore, op. cit., p. 142. 
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(6,7) V (4,8) 
(e, 8) 
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(b, v) 
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Gy 
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XJ. AN EXAMPLE 


To illustrate the ease of application of the preceding algorithm, even in the case of a complicated machine, we 
now apply it to the following connection matrix of a machine with 13 states, 6 inputs a, b, c, d, e, f and 4 outputs 


4 5 
0 0 
0) 0 
0 0 
(6,7) 0 
(e, 8B) (2, y) 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
(Ory) 0 
0 0 


5 


(c, B) 
(a, a) 
(e, 6) 


0 
0 
0 
0 


(a, a) V (¢, B) 


(0, v) 
(c, B) 
(c, 8) 


We 13 

0 Ota 

0 0 

0 0 
» 0 

0 (6, B) 
(d, a) 0 

0 0 

0 0 

0) 0 

0 0 
(¢8)) 0 
(a, a) 0 
(e,a) (a, B) 


Examining the rows of C, we first partition the states into the following maximal, disjunctive sets 


Sey oO fr) 


5 


BS UF has Hae 


Sf Selo 
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such that each submatrix formed from the corresponding rows of each S; is a 1 matrix. Next, partition C with re- 
spect to the S;. 


1 3 4 7 9 10 12 , 5 8 11 6 13 
0 (6, v) 0) 0 (a, a) V (c, 8) 0 0 | (e, 5) 0 0 0 {| O (= 
(e, 5) (6, y) 0 0) (a, x) 0 0 : 0 0 0 (c, B) : 0 0) 
(a, &) Oe AOsy)i 0 (e, 5) 0 0 0 0 0 (c, B) | 0 0 
(A)V@,s) 0 OO 0 0 (he ee Oe) (Nm ash ae ea 
(a, a) (6,7) 0 0 (c, B) 0 0 | 0 Om (ers) 0 : 0 0 
(e, 4) 0 0 0) (Gj) (a, a) 0) | 0) 0 (0, y) 0 | 0) 0 
0 0 (7 0 (c, 8) 0 (a) | (,8) 0 om 
0 Oe iS Cae 6 Gia) wo % (ae. 20. SOM aa) | Peis 
0 (d,a) (eB) (a,8) 0 0 0 0 (6,7) 0 0 0 (c, B) 
(c, 8) (8) 0 a) Ov: Oech os Menten) os UN) ee 
(c, 8) 0 0 Ga) 0 0 (2,8) (7) V G8) 0 0 0 0 0 
Oe ey et a chy ad) I) 10 5 oe Os ave Uae Ghee Ds 
0 Cow OG. | ee Deetee) 0 lo: 0 A eee 
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Examining the rows of C once again, we form the following partitioning of the states into maximal, disjunct sets 


( 
é 


Sig = {3,4}, Su = {7,10}, “Sam {28,01}, Sa = oes cones 


again such that the submatrices determined by the S;; are 1 matrices. Now, repartition the connection matrix with 
respect to the S;,;’s: 


1 9 12 3 4 i Toxo 2 Sip eet 5 6 613 
er 0 anywiee) 0 | Ga 0 7 0 Oh eee 0 6! {Ps oe Tete 
9 (a, 8) (c, B) ee eo ee ee 0 Pe aed tan Gd Mother 4 
12 0 (c, B) (ee 0 a 0 (e, 8) 0 Oa ee | 0 0 
Chee ie - 0 (ey 0 fo "OE ey Poe i Sua mee 
4 (a, @) (e, 6) 0 0 (6, y) 0 0 | 0 0 CASTER : 0 0 
Wie eas) a “ry 6. © Oy Sioa OSD) hy lOe beeen 
10 (e, 6) (<c, B) 0 0 Go: Pa | 0 ee) aie Vile Si | 0 0 
Bie = Ga Ou Gap ooh. 0 Geli cell. 0; Oe 0a) ome 
8 (c, 8) 0 uh (e, B) 0 | (dja) O | (a, B) 0 (b,y) | 0 7 0 0 
1 (c, 8) 0 GA 1 0. 0 Marl OR Oat eae SIME See eet etiam 
ie a = Ra Sah sun ee grea hire ae 
Bemese © 0 aa Gn” OF 0 Saat Marron De niet care mene 
ieee. 9.70 pial thai 0 (ee 0 0 O° [: 091 0 arem 
Since the submatrices are Hotell 1 matrices we form a new partitioning of the states: 
Siu = ae 2. 12 Sis. = eh A}, Si3i1 = a 10}, Seis uke, si, Sore = Bios Soar = io S311 = {6, 13}4 
Repartitioning the connection matrix with respect to the S;;, gives: 
1 9 12 3 4 (ee 2 8 11 5 eo ask? 
be 0 Gra) Vic, B) =~ 9 1 Gy) 0 0 ae PA A aa NE SY ee 
9 (2, a) (c, 8) Uy ren ie Nilect ep a Oan Oe 
12 0 (¢, B) (a) | ) (6, v) 0 0 | (e, 6) el Mog allan a | 0 0 
ea wie Views wate G = aia ea aca 
4 ee) (c, 8) O Eeoeth City): P Obwaiaro. | Mac iets gs 
7 | BV €,8) ‘Geet to. pe. can Cow Ge) eam 6 GN Coe ears 
10 | @a) CL ee MR Ce MRC tet rs 
i Fo ey a Se es Cie sai Go, 0, 0 gm 
8 (c, B) 0 0 (e, 8) 0 | (dja) 0 : (a, 8) 0 1,7 0 : 0 0 
MC aa mE NEM Mee PIR Re ee Re boas 
5 7 rs nes rem Cn tae me cE ro ee | 
ell@s | 2 O° G@Galyn 0 0 66 Se 0 Unig om on wn rn 
welt e Ga) — ta) | G9) ~ 0 ae, Oman eng 0 % obt om tana em 
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Since all submatrices are 1 matrices, the states corresponding to each Si;, are equivalent states. 
The reduced connection matrix is determined by replacing each submatrix by a single element which is the 
union of all elements in each submatrix. We have then the matrix 


1 3 
1 | (aa) V (¢,8) (0, v) 
3 | @e)V (e, 8) (2, v) 
: 7 | ©, 8) V ©, 8) 0 
Cre (c, 8) (e, B) 
11 | (¢,6) V (,8) 0 
5 0 (d, a) V (¢, 8) 


6 |_(@,a) V (a) (f, y) 


which defines the simplest equivalent machine. 


XII. ConcLusion 


In this paper we have developed and illustrated a 
systematic method for reducing a given state diagram of 
+ sequential machine to its simplest equivalent form. 
The method is so organized that it could readily be pro- 
crammed for a computer, which would reduce the labor 
significantly in problems of design. (The input-output 
pairs could be identified by the integers 1, 2,---, mp 
for this purpose.) 

In machines in which input restrictions exist one can 


4 


Correspondence 


7 2 i Sie 16 
0 (c, 8) Q: .0- Cane 
0 0 Cy 40 0 
(Cae) Oo cir oe 
(d, «) (2, 8) Gy) 0 
(d,a) (,7)V(@,8) 0 0 0 
(2, 6) 0 O.-> (x4), (G8) 
(b, 8) 0 0 O(a, 8s 


at times combine nonequivalent states without altering 
the input-output behavior of the machine, provided ap- 
propriate input restrictions are observed. This problem, 
and the problem of synthesis, will be treated in subse- 
quent papers. 
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The Logical Combination of 
Punched Paper Tapes* 


The early design of a digital computer 
“4¥ call for its output section always to sup- 
~’ such “extras” as line numbers, plus and 
‘us signs, tabs, and carriage return, in 
“* to bring the numerical results of com- 
petiole into a form that ultimately can be 
“ved to typewritten copy. This strict edi- 
* Tact ion may hamper the programmer in 
“ort to meet the requirements of an un- 
format, 
‘ sider, for example, the problem of 
‘ing a typewritten number field one- 


* Keceived by the PGEC, August 3, 1957. 


hundred and twenty characters wide from 
computer-edited output consisting of ten 
useful characters per word. Algebraic signs, 
word endings, and tabs must be deleted to 
adjoin the proper sequence of characters 
without overflowing the carriage of the tape- 
controlled typewriter. Fig. 1 gives an idea 
of such copy as it ordinarily appears from 
the computer. Fig. 2 shows the contour map 
that was obtained from this copy by voiding 
the signs, eleventh digits, and between-word 
tabs during the operation of transcription. 
This reduction was accomplished easily by 
having the automatic typewriter (e.g., tape- 
controlled Flexowriter) read an overlay tape 
along with the answer tape while typing the 
information. 


This procedure may be explained as fol- 
lows. The representation of each output 
symbol is by a six-hole binary code punched 
in a traverse row across the paper tape. A 
smaller feed hole is punched continuously 
near the center of the tape for the purpose 
of guiding the tape past the punched tape 
reader. In reading two paper tapes simul- 
taneously, one above the other, normally 
open contacts for each bit are closed when- 
ever a hole is sensed common to both tapes 
for that unit of the code. This is equivalent 
to finding the logical “and” of the two tapes. 
To have the automatic typewriter behave 
as though recognizing a delete, the corre- 
sponding position on the overlay tape is not 
punched. Consequently, the two tapes have 
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ooo0 
0000 00000000000 +00000000000 00000000000 +00000000000 
0001 +00000000000 +00000000000 +00000000000 +-00000000000 
0002 +00000000000 +Q00000000000 +00000000000) +00000000000 
0003 00000000000 +00000000000 +00000000000 +00000000000 
0004 +00000000000 00000000000 +00000000000) +-00000000000 
0005 +00000000000 +00000000000 +00000000000 +00000000000 
0006 +00000000000 +00000000000 +00000000000 +0000000001f 
0007 +00000000000 +00000000000 +00000000000 +0000001111F 
0008 +00000000000 +000LV0000000 +00000000000 -O0011111I1If 
0009 +00000000000 +00000000000 +00000000000 -+illllll22e¢t 
0010 +00000000000 +00000000000 +00000000000 -+#llllll22ze2r 
0011 +00000000000 +000009000000 +00000000000 +1111122222f 
0012 +00000000000 +00000000000 +00000000000 +1111222223¢ 
0013 +00000000000 +00000000000 +00000000000 +1111222233f 
0014 +00000000000 +00000000000 +00000000006 +1111222333f 
0015 +00000000000 +00000000000 +00000000000 +41112223333f 
0016 +00000000000 +00000000000 +00000000000 +¢1112222333f 
0017 400000000000 +00000000000 +00000000000) +1111222333f 
0018 +00000000000 +00000000000 +00000000000 +111112222:% 
0019 +00000000000 +00000000000 +00000000000 +1111112 
0020 +00000000000 +00000000000 +00000000000 -llllllllzef 
0021 +00000000000 +000V0000000 +00000000000 -+#llllllllilet 
0022 +00000000000 +00000000000 +00000000000 #lllilliliiie 
0023 +00000000000 +00000000000 +00000000000 ¢11lllllliif 
0024 00000000000 +Q000UV000G00 +00000000000 +1lllllilliigf 
0025 +00000000000 +00000000000 +00000000000 +lllillliiitgt 
0026 +00000000000 +00000000000 +00000000000 +lililllllif 
0027 +00000000000 +00000000000 +00000000000 +OOLI1II1IIE 
0028 +00000000000 +00000000000 +00000000000 +O000011111IIF 
0029 +00000000000 +00000000000 +00000000000 +0000011111f 
0030 +00000000000 +00000000000 +00000000000 +0000001112f 
0031 +00000000000 +00000000000 +00000000000 +0000001122¢f 
0032 +000000G0000 +00000000000 +00000000000 +0000001122f 
0033 +00000000000 +00000000000 +00000000000 +0000001122f 
0034 +00000000000 +00000000000 +00000000000 +0000001122f 
0035 +00000000000 +00000000000 +00000000000 +0000001122f 
0036 +00000U00000 +00000000000 +00000000000 +0000001122f 
0037 +00000000000 +00000000000 +00000000000 +000000i122f 
0038 +00000000000 +00000000000 +00000000000 +0000011122f 
0039 +00000000000 +00000000000 +00000000000 +0000011112f 
0040 +00000000000 +400000000000 +00000000000 +0000011111F 
0041 +00000000000 +00000000000 +00000000000 +O0000001L1If 
0042 +00000000000 +00000000000 +00000000000 +0000000011f 
0043 +00000000000 +00000000000 +00000000000 +00000000000 
0044 +00000000000 +00000000000 +00000000000 +00000000000 
0045 


0200 
0001 
0002 
0003 
0004 
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just the sprocket hole in common, which 
indicates to the machine that no action 
whatsoever should be taken. To allow the 
meaningful characters to come through, all 
six holes are punched in the corresponding 
positions on the overlay tape. 

The overlay principle can also be used in 
checking input tapes for the accidental 
omission of punched codes if some arbitrary 
code, such as carriage return, follows a pre- 
determined pattern. The carriage returns, 
say, are replaced by six-hole punches on the 
overlay tape with the proper number of in- 
tervening sprocket holes to fill out the pat- 


tern. If the pattern is periodic, the overlay 
tape may be cemented to form a loop for 
convenience. Assuming the punch to be “on” 
and a sheet of paper wrapped completely 
around the platen for recording‘ any error 
indication, the operation proceeds in the 
following manner. In the region where 
no mistake has yet occurred, the six-hole 
punches on the overlay tape will align on 
the carriage returns of the input tape with 
the resulting intersection causing the car- 
riage to advance and a carriage return code 
to be punched. The remaining symbols of 
the input tape will be placed in one-to-one 
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sprocket holes on the overlay tape causing 
both tapes to be read at a slightly faster rate 
than normal with no printing or punching. 
On the first occurrence of an omission, how- 
ever, the tapes will break out of phase caus- 
ing one of the remaining characters to print 
and be punched at each mismatch. Such a 
discrepancy then can be readily detected, 
its exact location determined, and the proper 
correction made. 

In addition to carrying blanks to require 
nonperformance of given operations and six- 
hole punches to allow reading of meaningful 
characters, the overlay tape may contain 
other six-hole symbols, in accordance with 
the binary code in use, so that variations 
appear in what normally would occur as the 
punched tape reader forms the logical inter- 
section of the two tapes. Thus for example, 
the hexadecimals 0, 1,---,9;a,b,--+,f 
might be contrived to go into the symbols 
for delete, color shift, space, carriage return, 
or any of the available type characters, thus 
giving the programmer still more control 
over the format. 

The overlay tape may be punched either 
by manual keyboard operation or by special 
perforating equipment. Somewhat similar 
effects may be obtained by placing the an- 
swer tape under the punch and the overlay 
tape on the reader. The resulting superposi- 
tion is the logical “or” of the two tapes. It 
is also possible that the operations of tape 
union and intersection described could be 
applied ina mechanical way to derive truth 
tables for Boolean expressions in several 
variables. 

Rosrrr M. MAson 
U.S. Naval Research Lab. 
Washington, D. C. 


Demonstration of Conditional Sta- 
bility on an Analog Computer* 


It is instructive to demonstrate experi- 
mentally a conditionally (Nyquist)-stable 
circuit. This is conveniently done on an 
analog computer. The setup turns out to be 
quite critical and some preliminary calcu- 
lations must be carried out. This note pre- 
sents the results. 

Fig. 1 gives a Nyquist diagram of a con- 
ditionally stable system. If the critical 
point (—1-+07) is at A or at C, the feedback 
system whose loop transfer function is H is 
stable; if it is at B, the system is unstable 


Fig. 1—Nyquist diagram of conditionally stable 
feedback system. 


Starting with a loop gain which corre- 
sponds to B as critical point, and therefore 
to an unstable system, an tucrease in gain 
can shift the critical point to C and stabilize 
the system. 

Perhaps the simplest transfer function 
exhibiting conditional stability is ~ 


He) = (1) 


The locus of I (jw) cuts the real axis at 
In {H (jo) } = 0. (2) 


The corresponding angular frequencies 
Wi, W2 aTe 
RP — 6k +34 Vk— 1IR-—9) . 
Oy ,2 (3) 
2a?k? 
‘Therefore, >9 for H to correspond to Fig. 
1. On the other hand, k must not be much 
larger than 9 if A (jor) (jor) <~10. The 
frequency range for the experiment is deter- 
mined by a, 
1.3 0.070 
oe oS (4) 
Sat Rye Viif2 
where the last expression follows for k=9. 
The following values correspond to Fig. 
2-Fig. 4: k=9.2; Vfi-fe=14 cps. From (4) 
@=35-10-3 and from (3) wi=100_ cps, 
wi=38 cps. Eq. (1) now gives / (je) 
=—| 7H (jo.)| =—1.17X10 and IT (jos) 
= —2.66X10~, corresponding to loop gains 
of A1= —1/fF7 (/jo) =85and A2= —1/H (je) 
=38. These are the gains at the transitions 
‘rom the unstable region to the stable re- 


The complete setup is given in Fig. 4. 
‘ = R;C;=5 107 sec; k= Re C;/RiC; = 0.046 
«=9.2. The maximum loop gain is 153 


* Received by the PGEC, June 7, 1957. 
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Fig. 2—Nyquist diagram of analyzed system, 
drawn to scale. 


Fig. 3—-Part A of Fig. 2, enlarged. 


Fig. 4—-Computer setup for system simulation. 


The results are best observed in the 
phase plane representation on a cathode- 
ray oscilloscope. Horizontal and vertical de- 
flection inputs are taken from the outputs of 
two of the computing amplifiers of the set- 
up. A stable system is recognized by an in- 
ward going spiral, an unstable system by a 
growing spiral. 

Amos NATHAN 

Yona MANLER 

Faculty of Elec. Eng. 

Technion, Israel Institute of Technology 
Haifa, Israel 


On the Use of Redundant Integra- 
tors in Analog Computers* 


One type of instability in analog com- 
puter setups may be traced to the use of 
integrators whose total number exceeds the 
degree of the characteristic equation being 
simulated. This phenomenon has been re- 
ported before! but appears not to be well 
understood. Ideally, the use of redundant 
integrators should not affect the solution, 
but in practice it is found that extra inte- 


* Received by the PGEC, October 11, 1957. 

1L. G. Walters, “Hidden regenerative loops in 
electronic analog computers,” IRE TRANs., 
pp. 1-4; June, 1953. 
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grators cause perturbations in the locations 
of the roots of the characteristic equation 
and may introduce extraneous roots with 
large positive real parts. These statements 
are demonstrated in the simple circuit of 
Fig. 1. Its circuit equations are 


oleh s se 


A possible setup is shown in Fig. 2. 
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Now let us suppose that the integrator 
in the upper chain has a multiplying factor 
of 1/k instead of 1. The characteristic 
equation of the computer setup is then 


(le BS) aes 
eS ea Mae 
or 
S@ Pie te 


If k=1, the computer setup yields the same 
characteristic equation as the original net- 
work, but if k#1, the transient response 
will be of the form 


4 = AeW#K+)) instead of 1 = AcW—4#/?), 


This difference arises from the fact that 
the computer setup uses two integrators to 
do a job done in the network with only one. 
The total inductor current (41—72) equals the 
integral of one quantity, the inductor volt- 
age. In the computer, however, a is pro- 
duced by one integrator and 72 by another. If 
these are not perfectly matched, (41 —72) will 
not be correct. 

Mere dislocation of the roots of the 
characteristic equation, however trouble- 
some it may be, is likely to occur even in 
the optimum computer setup due to com- 
ponent inaccuracies. These inaccuracies 
become even more troublesome, however, in 
setups which employ redundant integrators 
since there is a very real possibility that 
the degree of the characteristic equation 
may be increased with the result that not 
only are the desired roots moyed but new 
ones are introduced. For example, suppose 
that the summing resistor to which 
—dio/dt is fed at the input of the upper 
chain is of such a value that this input is 
multiplied by a factor & instead of by unity. 


The characteristic equation of the setup is 
then 


1+S —kS 
(1 -F S) kS | ie 
aS) (beers 


or 
Sub = fy) + 25+ 1 = 0. 


The roots are 


for kA~1. If & is precisely one, then no 
trouble arises since the characteristic equa- 
tion is of the first degree with its root where 


OGatiitotes. = ee 


D. D. Aufenkamp, for biography see 
page 205 of the September, 1957 issue of 
these TRANSACTIONS. 
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Joseph O. Campeau was born in Cleve- 
land, Ohio on March 10, 1930. He received 
the B.S. degree in clectrical engincering from 
Case Institute of Technology, Cleveland, 
Ohio in 1952 and the M.S. degree from the 
University of California at Los Angeles in 
1955. At the present time, he is studying for 
the Ph.D. at the University of California at 
Los Angeles. - 

From 1952 through 1956, he was engaged 
primarily in the design of digital computers 
at Hughes Aircraft Company, Culver City, 
Calif. Since 1956, he has been employed at 
Litton Industries, Beverly Hills, Calif. His 
work there is concerned with systems and 
digital computer design. 
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C. K. Chow was born in Hong Kong, on 
December 28, 1928. He attended the Utopia 
University and National Tsing Hua Uni- 
versity, China and was awarded the B.S.E.E. 
degree in 1949. He continued his education 
in the United States, and received the 
M.E.E. and Ph.D. degrees from Cornell 
University, Ithaca, N. Y., in 1950 and 
1953, respectively. 

From 1953 to 1955, he was an assistant 
professor in electrical engineering at the 
Pennsylvania State University, State Col- 
lege, Pa. Since 1955, he has been associated 
with the Burroughs Corporation Research 
Center, Paoli, Pa. 

Dr. Chow is a member of Sigma Ni, Pi 
Mu Epsilon, and Phi Kappa Phi. 


D. W. Davies was born in Treorchy, 
Wales, on June 24, 1924. He attended the 
Imperial College of Science and Technology, 
London, Eng. and was graduated in 1943, 
with the degree of B.S. in physics. In 1946, 
he received the B.S. degree in mathematics 
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it belongs at S=(—1/2), If k is slightly less 
than unity, the trouble is not severe since 
one of the roots is a large negative number, 
corresponding to a rapidly decaying ex- 
ponential term which can be neglected, and 
the other root lies very close to S=(—1/2). 
However, & is as likely to be slightly larger 
than one, and in this case the extraneous 
root is found to lie far out in the right-half 
plane. This root has a disastrous effect on 
the solution of the problem since it corre- 
sponds to a rapidly growing exponential 
term. 

Such artifices as judiciously inserting 
additional damping to counteract the effects 


from the same school, During World War 
II, he was engaged in theoretical work on 
the British atomic bomb project. In 19-6, 
he joined the National Physical Laboratory. 
He took part in the design and construction 
of the ACE pilot model which began to 
work in May, 1950. At the present time, 
Mr. Davies is concerned with the design of 
the ACE computer. 


Robert J. Domenico was born in New 
York, N. Y., on August 1, 1926. He was 
graduated from Syracuse University with a 
B.S.E.E. degree in 1951. 

He joined IBM in June, 1951, where he 
worked on the design of circuits for a gas 
tube counter, the design of a fast rise time 
pulse generator, and the design of a video 
amplifier for electrostatic memory. In 1954, 
he was assigned to the transistor circuit 
group to design logical switching circuits. 
In 1956, he was appointed to head a group 
in the Poughkeepsie Product Development 
Laboratory to study methods of mechaniz- 
ing the design of transistor switching circuits. 
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Paul C. Dow, Jr. was born in Melrose, 
Mass. on March 31, 1927. He received the 
B.S. degree from the United States Military 
Academy from which he graduated as a 
second licutenant in the Air Force in 1949. 

After completing flying training and 
serving two years as a flying instructor, he 
attended the guided missile course at the 
University of Michigan, receiving the 
M.S.E. degree in 1954. He served two years 
as a graduate assistant in acronautical 
engineering at the Air Force Institute of 
Technology while completing his doctoral 
dissertation, and received the Ph.D. degree 
from the University of Michigan in 1957. 

Captain Dow is presently assigned as a 
project officer at the Air Force Ballistic 
Missile Division of the Air Research and 
Development Command, Inglewood, Calif. 

He is a member of Phi Kappa Phi, 
Sigma Ni, Institute of the Aeronautical 
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of spurious roots have been suggested, but 
these fail to recognize that the additional 
roots are present because the computer 
setup uses in its closed loop a number of 
energy-storage elements (the integrators) 
which exceed the degree of the characteristic 
equation of the problem. Extra integrators 
not included in the closed loop cause no 
trouble and may sometimes be required in 
cases where the minimal integrator setup to 
implement the equation yields only deriva- 
tives of some of the desired quantities. 
Norman R. Scorr 
University of Michigan 
Ann Arbor, Mich, 


Sciences, American Society for Engineering 
Education, and the National Society of 
Professional Engineers. 
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Franz E, Hohn, for biography sce page 
206 of the September, 1957, issue of these 
TRANSACTIONS. 
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Robert E. Horn (S’52-A’55) was born in 
Richmond Heights, Mo., on November 30, 
1927. He received the B.S.E.E. degree from 
Washington University, St. Louis, Mo., in 
1950. Awarded a Fortesque Fellowship for 
graduate study, he received the M.S.E.E. 
degree in 1951. Previously, he served two 
years as an electronic technician with the 
U.S. Navy. 

From 1951 to 1955, he was an instructor 
in the electrical engineering department of 
Washington University. In addition to his 
teaching duties, he was responsible for 
installation and operation of the depart- 
ment’s analog computing equipment. After 
obtaining the doctorate degree in 1955, he 
was employed by the Air Arin Division, 
Westinghouse Electric Corporation, Balti- 
more, Md., as a senior engineer engaged in 
airborne fire-control system studies. In 
his present position as a fellow engineer, he 
is responsible for the activities of an analy ti- 
cal group working on missile seeker and. 
bomber defense problems. 

Dr. Horn is a member of the AIEE, As- 
sociation Internationale Pour le Calcul 
Analogique, Tau Beta Pi, and Sigma Xi, 
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Verne G. Fauque was born in Kevin, 
Mont., on June 7, 1931. He received the 
B.A. degree in mathematics in 1953 and 
the M.A. degree in 1954 from Montana 
State University at Missoula, 

In 1954, he joined the Westinghouse 
Electric Corporation’s Graduate Student 
Training Program. Upon completing this 
program in 1955, he was assigned to the Air 
Arm Division, Baltimore, Md., where he is 
now engaged in fire-control system studies. 

Mr. Fauque is a member of Pi Mu 
Epsilon. 
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MEMBERSHIP DISTRIBUTION OF 
PGEC By SECTIONS 


The membership of the PGEC increased 
1226 from January 1, 1957, to August 15, 
1957, when the total paid membership was 
6551. Noticeable increases were seen in 
practically every section with outstanding 
gains in foreign memberships. 

The distribution of paid members as of 
August 15, 1957, is shown as follows. 


In TRE Sections 


Akron 34 
Alamogordo-Holloman 13 
Albuquerque-Los Alamos 38 
Atlanta 47 
Baltimore 114 
Bay of Quinte 3 
Beaumont-Port Arthur 2, 
Binghamton 62 
Boston 506 
Buenos Aires 3 
Buffalo-Niagara 31 
Cedar Rapids 25 
Central Florida 62 
Central Pennsylvania 16 
Chicago 189 
China Lake 10 
Cincinnati 41 
Cleveland 42 
Columbus 18 
Connecticut Valley 121. 
Dallas 36 
Dayton 50 
Denver 29 
Detroit 138 
Elmira-Corning 4 
El Paso 21 
Emporium 4 
Evans-Owensboro 6 
Florida West Coast § 
Fort Huachuca 1 
Fort Wayne 17 
Fort Worth 31 
Hamilton 5 
Hawaii 2 
Houston 46 
Huntsville \ 34 
Indianapolis 23 
Israel 14 
Ithaca 14 
Kansas City 24 
Little Rock 6 
London, Ont. if 
Long Island 330 
Los Angeles 891 
Louisville 7 
Lubbock 4 
Miami 14 
Milwaukee 54 
Montreal 47 
Newfoundland 1 
New Orleans LS 
New York 748 
North Carolina-Virginia 3 

Northern Alberta 4 
Northern New Jersey 304 
Northwest Florida 4 
Oklahoma City 13 
Omaha-Lincoln § 
Ottawa 18 
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Philadelphia 512 Dallas, Texas, on October 23-25, 1958. Full 
Phoenix 46 information regarding submission of papers 
Pittsburgh 59 and other details will be published in a 
Portland 15 future issue of these TRANSACTIONS. 
Princeton BY 
Regina 0 
Rochester 39 
Rome-Utica 35 
Sacramento 6 
St. Louis Sidi 
Salt Lake City 8 SuMMARY OF PGEC ApMINISTRA- 
San Antonio 20 es ae aa 
SaniDicc 59 TIVE COMMITTEE MEETING 
San Francisco 341 The PGEC Administrative Committee 
Schenectady 31 met during WESCON in San Francisco, on 
Seattle 53 August 22, 1957. A quorum of 14 members 
Shreveport cs and proxies out of a total of 17 members 
South Bend-Mishawaka 11 were present. 
Southern Alberta 3 Secretary-Treasurer H. W. Nordyke 
Syracuse 77 reported membership and financial figures 
Tokyo 30 as of June 30, 1957: 
joes Members........ 5998 
oronto 47 

Aion 12 Studentseseerrets 544 
Tulsa Br 19 Receipts (January 1-June 30) $16,776.44 
Twin Cities 134 Expenditures (January 1- 
Vancouver 10 Juné-tO)e4asc uae qe $ 9,339.46 
Washington 301 Balances (une 0) eyecare $26,356.99 
Wichita 9 

. Winnipeg 4 Editor J. P. Nash reported that H. D. 


In Other Localities 


Alaska 1 
Argentina 3 
Australia 5 
Austria ey 
Belgium 7 
Bermuda 1 
Brazil 4 
Cuba 1 
Czechoslovakia 2 
Denmark 4 
England 40 
France 21 
Germany 10 
Greece 1 
Hong Kong 1 
India 4 
Italy 17 
Japan 30 
Lebanon 1 
Mexico 3 
Netherlands 10 
Norway 2 
Scotland 1 
Sweden 22 
Switzerland 7 
Turkey 2 
Venezuela 3 
USSR 1 

Total Paid Menibers 6551 


1958 SIMULATION CONFERENCE 


The 1958 National Simulation Confer- 
ence, sponsored jointly by the PGEC and 
the Dallas IRE Section, will be held in 


Huskey is resigning as editor of the Review 
of Current Literature section in the TRANS- 
ACTIONS, effective January, 1958. L. F. Jones, 
Chairman of the Baltimore Chapter, and 
R. W. Melville agreed to cooperate in 
organizing a Bibliography Committee which 
would recommend action on the reviews as 
well as on a possible bibliography of com- 
puter literature. 

On the motion of R. W. Melville, Stu- 
dent Relations Chairman, $100 was voted 
to his committee for working with high 
schools. 

In conformance with IRE policy, a 
citizenship restriction which was inad- 
vertently placed last year on the IRE 
Computer Fellowship will be removed in 
the future. 

The AIEE has been added to the list of 
societies whose members may become PGEC 
affiliates. 

The Administrative Committee decided 
that the PGEC Vice-Chairman should 
automatically be Chairman of the Awards 
Committee each year. This action assures 
that awards activities will continue to be 
in the hands of a member with necessary 
experience in PGEC affairs. 


Other Committee appointments are: 


R. W. Melville—Student Relations 

L. C. Nofrey—Chairman, Membership 

R. T. Blakely—Vice-Chairman, Member- 
ship 

J. C. LaPointe—Constitution and By-Laws 

N. M. Blachman—Publications 

S. B. Disson—Sectional Activities 

W. D. Winger—Tcchnical Program Repre- 
sentative, 1958 IRE National Convention 

J. P. Nash—TrRAnsActions Editor. 


W. BucunoLz 
Chairman 
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Reviews of Current Literature 


It is the intention of this section to review articles that have been published since January 1, 
1953, and to publish eventually reviews of all books pertaining to the computer field. Authors can 
be of considerable assistance in this review process by sending two reprints of their articles to 
H. D. Huskey, Department of Electrical Engineering, University of California, Berkeley, Cali- 
fornia. The editors wish to express their gratitude to the reviewers who, through their efforts, make 
this section possible—H. D. Huskey 


GENERAL 


57-159 
The Computing Revolution—E. L. 
Harder. (Elec. Eng., vol. 76, pp. 476-481, 
586-590; June, July, 1957.) This two-part 
article outlines some of the significant de- 
velopments in the computing field during 
1956. Subjects discussed include computer 
developments, engineering applications, au- 
tomatic programming, simulation, analog 
vs digital computers, computer uses in 
business, and control applications. The use 
of computers as an engineering and research 
tool is examined, and the growing uses in 
business and manufacturing are noted. 
H. T. Larson 


57-160 
The Position of the University in the 
Field of High Speed Computations and Data 
Handling—A. S. Householder. (Computers 
and Automation, vol. 5, pp. 6-10; May, 
1956.) In a rather rambling manner, the 
author describes the past and ‘present 
activities of European and American uni- 
versities, followed by a summary of his 
opinion of the proper future position of 
universities with respect to this subject. The 
main theme is that specialized courses in 
the programming of digital computers and 
numerical analysis are not the proper ap- 
proach. He prefers the complete reorganiza- 
tion of the mathematics curriculum in a 
manner which integrates numerical methods 
with other techniques starting with the 
first-year calculus courses. It is the author’s 
feeling that, if this is done, the student can 
receive a more sophisticated training in this 
technique than occurs in one semester 
courses which can only provide a collection 
of recipes. In this manner the student 
would be better prepared to obtain solutions 
to his mathematical problems in the field 
in which he is engaged. 
Gordon E. Morrison 


57-161 

The Convention on Digital-Computer 
Techniques—(J. TEE, vol. 2, pp. 475-489; 
August, 1956.) These pages contain informa- 
tion on a Convention on Digital-Computer 
Techniques, April 9-14, 1956. Included are 
a three-page suniunary of the convention 
and reviews of the following papers: “Digital 
Computers and the Load-Flow Problem,” 
by J. M. Bennett; “Whe Use of the Pilot Ace 
for Testing a New Design of Proton Syn- 
chroton,” by G. G. Alway; “The Application 
of Digital Computers to Electric Traction 
Problems,” by A. Gilmour; “The Use of 
Digital Computers in Obtaining Solutions 


to Electric-Circuit Problems Involving 
Switching Operations,” by S. J. M. Denison 
and D. G. Taylor; “Sorting of Data on an 
Electronic Computer,” by D. W. Davies; 
“Application of Digital Computers in the 
Exploration of Functional Relationships,” 
by G. E. P. Box and G. A. Coutie; “An 
Electronic Calculator for Punched-Card 
Accountancy,” by L. Knight; “The Pro- 
gramme-Controlled Computer,” by E. J. 
Guttridge and R. P. B. Yandell; “The HEC 
Computer,” by R. Bird; “The Manchester 
University Mark II Digital Computing 
Machine,” by T. Kilburn, D. B. G. Ed- 
wards, and G. E. Thomas; “A Transistor 
Digital Computer,” by E. H. Cooke- 
Yarborough, R. C. M. Barnes, J. Stephen, 
aiid G. A. Howells; “Transistor Arithmetic 
Circuits for an Interleaved-Digit Com- 
puter,” by R. C. M. Barnes, G. A. Howells, 
and E. H. Cooke-Yarborough; “The Digital 
Computer as an Aid to the Electrical Design 
Engineer,” by B. Birtwistle and B. M. Dent 
(See Rev. 57-218, this issue.) 

Dee shart 


681.142 57-162 

Application of Computers in Automatic 
Systems-—A. A. Fel’dbaum. (Awtomatika 4 
Telemekhanika, vol. 17, pp. 1046-1056; 


November, 1956.) A survey. Twenty-six 
references including several to Russian 
literature. 


Courtesy of Proc. IRE 
and Wireless Engineer 


57-163 

Analog Versus Digital Techniques for 
Engineering Design Problems—D. B. 
Breedon. (IRE'TrRANs., vol. PT-2, pp. 86- 
89; April, 1957.) Nearly every problem 
encountered in engineering at some time 
proceeds from the qualitative to the quanti- 
tative phase where the results of mathe- 
matical analysis must be applied in actual 
computation. Most often the computation 
is short enough that automatic means are 
not necessary. However, more and more 
problems are requiring powerful aids to 
calculation. This increase is due as much to 
expanded thinking encouraged by the mere 
availability of computers as to any actual 
backlog of work. Therefore it is to the 
engineer’s advantage to know what com- 
puters can do for him, even though he may 
take his problem to someone else for final 
preparation and programming. The follow- 
ing text presents some examples in which 
automatic calculation is being used. The 
logic used in choosing the computing meth- 
ods is shown based on the characteristics of 


problem and computer. As background for 

the examples the most important of these 
characteristics are presented briefly. 

Author’s Summary 

Courtesy of Proc. IRE 


57-164 
A Shaft-to-Digital Encoder—B. M. 
Gordon, M. A. Meyer, and R. N. Nicola. 
(Proc. Western Comp. Conf., February 11-12, 
1954, Los Angeles, Calif.; pp. 128-133; 
April, 1954.) A miniature commutator-type, 
analog-to-digital converter is described. 
Coding ambiguity is avoided by using two 
sets of commutators driven through an 
intermittent drive system and mechanically 
phased so that one commutator set is always 
stationary and ready for reading while the 
other commutator set is being positioned. 
Provision is made for connecting two input 
voltage sources to represent the zeros and 
ones of the parallel output. An auxiliary 
commutator is used to switch the input 
voltages to the stationary readout-commu- 
tator set. All commutators are cylindrical, 
one-quarter inch in diameter. Over-all size 
of encoder is 2} inches long and 14 inches in 
diameter. 


R. A. Gaskill 


57-165 
Digital Compensation for Control and 
Simulation—Julius Tou. (Proc. IRE, vol. 
45, pp. 1243-1248; September, 1957.) A 
very brief review of the Z-transform method 
is given. A mathematical treatment is 
given, illustrated by an example to show 
that 1) a digital-analog feedback system can 
be stabilized by digital techniques, 2) the 
system error can be made zero at the sam- 
pling instants by digital compensation, and 
3) the programs required for 2) can be made 
realizable (7.e., to require no foreknowledge) 
by a delaying program. The technique of 
synthesizing programs from Z transforms is 
not discussed nor are the physical require- 
ments of the resulting programs (program 
steps, computing and sampling rates). The 
limitations of the technique do not seem to 
be clearly described. 
A. D. Scarbrough 


57-166 

Temperature Dependence of Junction 
Transistor Parameters—Wolfgang W. Gart- 
ner. (Proc. IRE, vol. 45, pp. 662-680; May, 
1957.) This paper presents calculated typical 
values of a wide assortment of transistor 
parameters as a function of temperature for 


four representative kinds of transistors 
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—The Editor 


(Ge p-u-p alloy, Ge n-p-2 grown, Ge n-p-n 
rate grown, and Si 2-2-2 grown). The curves 
presented should prove particularly valuable 
as indicators of expected trends to those 
engineers designing transistor circuits for 
military applications. : 

A. D. Scarbrough 


681.142.002.2 57-167 
Factory for Electronic Digital Computers 
——(Lngineer, London, vol. 202, p. 593; Oc- 
tober 26, 1956.) Design features of the 
“Pegasus” and “Mercury” computers are 
discussed. 
Courtesy of Proc. IRE 
and Wireless Engineer 


ANALOG COMPONENT 


RESEARCH 
57-168 
_ An Electronic Analog Multiplier—D. C. 
Kalbfell. (IRE Trans., vol. EC-6, pp. 


100-103; June, 1957.) This multiplier uses 
the variable pulse area principle, but em- 
ploys phase sensitive circuitry to operate 
naturally in all four quadrants without bias 
voltages. The output is zero if either input 
is zero. The X and Y channels are separately 
linearized with independent feedback loops. 
The circuitry is simple and lends itself to 
either transistors or vacuum tubes. 
Author’s Summary 
Courtesy of Proc. IRE 


~ 57-169 
A Function Multiplie-—J. D. N. Van 
Wyk. (Nature, vol 178, pp. 1247-1248; 
December 1, 1956.) An electronic multiplier 
is described which has 1 per cent nominal 
accuracy. A rectangular area is displayed 
on the face of a cathode-ray tube. The 
illumination in each of four quadrants of 
the rectangle is detected by phototubes. If 
the displacement of the rectangular area 
in the X and Y directions is proportional 
to two voltages, the output of the photo- 
tubes can be summed to yield the product of 
these voltages. 


DEE. Hart 


57-170 

A New Diode Function Generator— 

T. Miura, H. Amemiya, and T. Numakura. 
(IRE Trans., vol. EC-6, pp. 95-100; June, 
1957.) With the diode function generators 
that are currently in use, generation of 
functions is made by combining straight 
lines. The main drawback of these function 
generators is that the slope of the line seg- 
ment cannot be changed independently. 
With the new function generator described in 
this paper, functions are generated by con- 
necting independent line segments. Accord- 
ingly, the slope of each segment is given 
independently and also quantitatively. It is 
possible to approximate any desired function 
without recourse to an oscilloscope for in- 
spection. These advantages are realized by 
using ganged potentiometers differentially. 
The operating principle and a_ practical 
generator experimentally built are described. 
Author’s Summary 

Courtesy of Proc. IRE 
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ANALOG SYSTEMS 


RESEARCH 

57-171 
Dynamic Accuracy as a Design Criterion 
of Linear Electronic-Analog Differential 
Analyzers—A. Nathan. (IRE TRANs., vol. 
EC-6, pp. 74-86; June, 1957.) A frequency 
error analysis of computing elements is 
presented which leads to a definition of 
their dynamic accuracy. The concept of a 
computing transfer function is introduced 
for this purpose, permitting the evaluation 
of an effective bandwidth, the latter being 
connected with the variance of the output 
for wide-band inputs. Limited bandwidth 
is considered as equivalent to finite resolu- 
tion and thus to an additional effective 
error. Single frequency errors are dealt with 
separately and are shown to be of minor 
importance. Suitable optimalization of 
dynamic accuracy yields parameters of 
design and performance such as optimum 
computing time and required base amplifier 
gain. The theory is applied to integrators 
and adders with base amplifiers of direct and 

of capacitive coupling. 
Author’s Summary 
Courtesy of Proc. IRE 


ANALOG EQUIPMENT 


57-172 

QC Computers for Machine Control, 
Parts 1 and 2—G. H. Amber and P. S. 
Amber. (Hlec. Mfg., pp. 80-88, 78-85, 298; 
July, August, 1956.) Analog computers for 
computing & and oa for either continuous or 
Ciscontinuous distribution are described. 
These computers can be used in automated 
systems to control quality automatically. 


G Courtesy of Data Processing Group 
General Motors Research Staff 
681.142 57-173 


A New Analogue Computer using Matrix 
Iteration for Determining the Roots of 
Algebraic Equations—J. Miroux. (Ann. 
Telécommun., vol. 11, pp. 226-232; No- 
vember, 1956.) The mathematical principles 
underlying an experimental computer are 
detailed. Its application and possible de- 
velopment are discussed. 

Courtesy of Proc. IRE 
and Wireless Engineer 


UTILIZATION OF ANALOG 


EQUIPMENT 
57-174 
Computers—The Key to Modern Manu- 
facturing Scheduling—J. P. J. Gravel and 
T. F. Kavanagh. (IRE Trawns., vol. PT-2, 
pp. 90-93; April, 1957.) Manufacturing 
operations with poor scheduling plans are 
headed for trouble. Load capacity analysis 
is a technique for measuring the feasibility 
and desirability of proposed scheduling 
plans. Simply stated mathematically, load 
capacity analysis consists of a series of 
multiplications and additions. However, the 
numerous computations in a typical prob- 
lem usually take more time than can be 
allowed. The paper describes a_ special 
purpose analog computer specifically de- 
signed to solve this problem in a matter of 
minutes. 
Author’s Summary 
Courtesy of Proc. IRE 


57-175 
Electronic-Analogue-Computer Study of 
Synchronous-Machine Transient Stability— 
A. S, Aldred and P. A. Doyle. (Proc. IEE, 
vol. 104, pt. A, pp. 152-160; April, 1957.) 
Authors’ Summary: The paper describes a 
method of solution of synchronous-machine 
transient-stability problems by a more exact 
electronic analog representation of the sys- 
tem equations than has hitherto been em- 
ployed. In the machine analysis emphasis is 
placed on the variation of field flux linkage 
during transient disturbances, and reasons 
are given why this is necessary. The equa- 
tions of performance are derived in the most 
convenient form for analog computing. The 
equivalent analog interconnections are 
given both for time varying and for constant 
field flux linkage. The components required 
for these analogs are considered briefly. The 
section on computer control deals with the 
control of integrators with reference to the 
introduction of “initial conditions,” “hold,” 
and “compute” states. A delay unit is 
described which enables the switching out 
of a faulted transmission line to be simulated. 
The results of various problems solved by 
the computer are presented in graphical 
form as boundaries between stable and un- 
stable states. It is shown how instability sub- 
sequent to the first rotor oscillation may 
occur, and an explanation is proposed to 
account for this phenomenon. Synchronous- 
machine and transmission-line analysis are 
excluded from the main text for clarity and 
are presented in detail in the Appendixes. 
D. E. Hart 


681.142 57-176 
Approximate Solution of Differential 
Equations with Partial Derivatives using 
Electrical Analogues—E. S. Kozlov and 
N. S. Nikolaev. (Awtontatika 1+ Telemek- 
hantka, vol. 17, pp. 890-896; October, 1956.) 
Analogs for solving Laplace, Poisson, and 
Fourier-type equations are briefly discussed. 
Courtesy of Proc. IRE 

and Wireless Engineer 


57-177 
Computing Techniques for the Sampling 
Parametric Computer—C. J. Hirsch and 
F. C. Hallden. (TRE Trans., vol. EC-6, 
pp. 108-119; June, 1957.) This paper de- 
scribes novel calculating techniques suitable 
for an electronic analog computer using 
exponential discharges to simulate the 
logarithmic scales of the slide rule. Among 
others, the device can perform the following 
Operations: YHNY, SZ == 0 / Ve Z=xy*; loge x;:a7; 
and evaluate such seriesas y= Ax*+ Bxr?+Cxe 
-+ ...where the exponents can be fractional 
and z, x, and y can be time variables. The 
problems can be solved explicitly or im- 
plicitly; thus in the above series y (or x) can 
be determined with equal ease if x (or y) are 
given. The appendix describes an actual 
computer which gave accuracies of 1 to 2 
per cent of full range from 0 to 95 per cent 
of full range and 4 per cent from 95 to 100 
per cent of full range without recalibration 
for different groups of problems. The ac- 
curacy can be increased by calibrating the 
device for a specific group of problems. 
Author's Summary 
Courtesy of Proc. IRE 


294. IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


This page has been left blank in order 


that readers may mount all reviews on cards. 
—The Editor 


= a ae 


1957 


57-178 
Trigonometric Resolution in Analog 
Computers by Means of Multiplier Elements 
—R. M. Howe and E. G. Gilbert. (IRE 
TrAns., vol. EC-6, pp. 86-92; June, 1957.) 
A method of generating sine and cosine 
functions in analog computers by means of 
multiplier elements and integrators is dis- 
cussed. Static accuracy of the method is 
analyzed and found to be essentially equal 
to the accuracy of the multipliers employed. 
The system accepts d0/dt as the input and 
gencrates output voltages of sin @ and cos @. 
Amplitude-stabilizing loops are employed 
to maintain sin? @+-cos? 6=1. Advantages 
of the method include representation of 
unlimited range in angles, dynamic capa- 
bilities far beyond that of the multipliers 
alone, and possibility of employing elec- 
tronic multipliers. The method has been 
successfully used to compute Euler angles 
in analog computer solutions of the three- 
dimensional flight equations. 
Author’s Summary 
Courtesy of Proc. IRE 


DIGITAL COMPONENT 
RESEARCH 


57-179 
Basic Logic Circuits for Computer Ap- 
plications—G. W. Booth and T. P. Both- 


well. (Electronics, vol. 30, pp. 196-200; 
March, 1957.) Basic digital computer 
~circuits designed around high-frequency 


transistors are described. The operating and 
design features of a flip-flop, gated pulse 
amplifier, dc amplifier, power pulse ampli- 
fier, and neon lamp indicator are discussed. 
The flip-flop and pulse amplifiers were de- 
signed to operate at a pulse repetition rate 
of 250 ke. The Eccles-Jordan circuit was 
chosen for the flip-flop because of its com- 
plimentary outputs. Saturated operation of 
the transistors is employed with collector 
current limited to 10 ma. Diode forward 
voltage drops and breakdown voltage of 
diodes and transistors dictated the choice of 
flip-flop output swing and gating levels of 
zero and minus 6 v. Further discussion 
deals with the type of logical structure for 
which the circuitry is best suited. Data re- 
garding operating expevience are presented 
and were gained through use of the circuits 
in connection with a general purpose com- 
puter which was specifically designed as a 
test facility. 

Norman F. Loretz 


57-180 

Magnetic Computer has High Speed— 

T. H. Bonn. (Electronics, vol. 30, pp. 156- 
160; August, 1957.) The ferroactor, a mag- 
netic amplifier used in the UNIVAC mag- 
netic computer, is capable of operating at 
frequencies as high as 2.5 mc. The basic 
‘ circuit is a series-type magnetic amplifier. 
Variations of this basic circuit produce 
complementing and noncomplementing am- 
plifiers. A carrier cycle is divided into input 
and output periods. The noncomplementing 
amplifier produces signals during output 
periods only when signals are applied during 
corresponding input periods. The comple- 
menting amplifier normally produces output 
signals when no input signals are applied. 
An input signal causes a no output during 


Reviews of Current Literature 


the following output period. The no-output 
current of the amplifier is absorbed by a 
constant current sink, thus suppressing the 
no-output signal. Details of circuit operation 
and interconnection are presented. Factors 
affecting the choice of 4-79 permalloy from 
among other materials are discussed. 
Norman F. Loretz 


57-181 
A 2.5 Megacycle Ferractor Accumulator 
—R. D. Torrey and T. H. Bonn. (Proc. 
Eastern Joint Comp. Conf., December 10-12, 
1956, New York, pp. 50-53; 1957.) Mag- 
netic amplifiers having power gains up to 
four at an information frequency of 2.5 
me are described. Two types are shown, a 
straight amplifier used primarily to provide 
a pulse time delay and an inverting type. 
Each amplifier is a series pulse type with a 
cycle of operation divided into an input 
period and an output period. There is an 
inherent pulse time delay of 0.2 usec through 
the amplifiers. Logical functions are per- 
formed by isolation diodes at the input of 
the amplifiers. By alternately defining a 
“one” as being the presence and then the 
absence of a pulse series chain of the invert- 
ing amplifiers will yield the logical configura- 
roms “Binel,” or” Sal” Sore? wane, JN 
blocking pulse of half the clock pulse 
amplitude is used at the output. This block- 
ing pulse disconnects the clock during the 
input phase and is arranged to have no 
effect during the output phase. The prob- 
iems of clock generation and distribution are 
discussed and detailed information is shown 
for the clock transformer. The demonstration 
model chosen is an accumulator consisting 
of a unit adder and a four-bit circulating 
register. Provision is made to add unity 
when a push button is depressed; the ac- 
cumulated count is displayed by means of 
incandescent lamps powered directly from 
the cores. Thirty-seven amplifiers are used 
in the model requiring a clock input power of 
approximately 15 w, or 400 mw per core. 
Lad Andrews 


57-182 

A Saturable-Transformer Digital Am- 
plifier with Diode Switching—E. W. Hogue. 
(Proc. Eastern Joint Comp. Conf., December 
10-12, 1956, New York, pp. 58-64; 1957.) 
This paper describes a digital amplifier 
employing a saturable transformer operating 
in the 100-kce to 300-ke range. The clock is 
a two-phase sine-wave voltage source of 
approximately 10-v amplitude. The over-all 
power consumed is 0.27 w per amplifying 
stage. The amplifier is suitable for use with 
two-level diode gating logic to provide “and” 
and “or” logical functions between amplify- 
ing stages. The “not” function is obtained 
by using two transformers to obtain a 
complementing action. A criterion for sta- 
bility, which insures that signals are not 
attenuated in long chains of amplifiers, is 
established. Factors affecting the gain of 
the amplifier are explored and the effects of 
winding reactances illustrated. An amplifier 
using a permalloy core transformer can 
drive 6 to 8 other stages, whereas ferrite core 
transformers are capable of driving only two 
others. A six stage circulating register with 
broad operating tolerances is described. 
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Also described is a proposed logical package 
containing an “and” and “or” input circuit, 
and both a direct and complemented output. 

R. D. Torrey 


57-183 

Multihole Ferrite Core Configurations 
and Applications—H. W. Abbott and J. J. 
Suran. (Proc. IRE, vol. 45, pp. 1081-1093; 
August, 1957.) A method is described using 
4-hole transfluxors and a special winding 
configuration to permit blocking and un- 
blocking pulses to have virtually unlimited 
amplitude. By threading the ac drive wind- 
ing through the transfluxor and a suitably 
designed single-hole core it is shown to be 
possible to use voltage source drivers. By 
using these techniques plus a winding tech- 
nique to compensate for lack of squareness 
in the ferrite hysteresis loop, the authors 
obtained operation from —50c to -+180c 
with an ac signal unblock-to-block ratio of 
at least 25 db. A number of winding con- 
figurations employing multihole cores to 
perform complex operations are described. 
For example, a 6-hole core (called a logicor) 
can be used to make a binary half-adder, a 
4-terminal odd parity checker, and other 
complex devices. The logical applications do 
not provide the compensating features 
which permitted wide temperature range 
operation of the 4-hole transfluxors. The 
experimental work used an ac signal fre- 
quency of 50 ke and pulse rates up to 20 ke. 
A. D. Scarbrough 


57-184 
Ferrite Apertured Plate for Random- 
Access Memory—J. A. Rajchman. (Proc. 
Eastern Joint Comp. Conf., December 10-12, 
New York, pp. 107-115; 1957.) For some 
time the idea of a continuous media for 
digital data storage has been considered. 
This is the first report of significant success 
in this endeavor. The development of small 
prototype ferrite plates with a 16 by 16 (256) 
aperture pattern is described and the be- 
havior pertinent to memory application is 
thoroughly discussed. Of particular impor- 
tance to memory technology, however, are 
the interesting circuit advances described 
for use in the memory array. A pair of plates 
is used for each digit plane, one storing the 
zeros and the other the ones. This provides 
noise cancellation, giving a positive signal 
for a one and a negative signal for a zero. It 
also provides a constant load on the drivers, 
since one core is always switched in each 
digit plane. A two-wire system is also 
described which requires a word sclection 
matrix in addition to the memory stack. 
One wire is plated on the ferrite and the 
other wire is threaded straight through all 
of the planes at once. This results in only 
one aperture pair in a digit plane being 
driven during the read operation which 
greatly reduces noise problems. The cost 
advantage attendant on this method is 
significant. Other interesting techniques are 
reported in the description of the memory 
media and access methods. 
Robert A. Tracy 


57-185 

A Compact Coincident-Current Memory 
—A. V. Pohn and S. M. Rubens. (Proc. 
Eastern Joint Comp. Conf., December 10-12, 
1956, New York, pp. 120-123; 1957.) This 


296 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


This page has been left blank in order 


that readers may mount all reviews on cards. 
—The Editor 


~ 


]957 


paper describes the present status of the 
work on evaporated magnetic film storage 
elements. It is a summary in nature as the 
extent of the work is much too large to be 
covered completely in a paper of this length. 
Within this restriction it is an excellent 
paper as the work is described concisely and 
accurately. The techniques described in the 
paper represent a major advance in the 
magnetic core storage art and it is fortunate 
that so competent a paper should have been 
written so early in the development of these 
techniques. The only slight criticism that 
can be made is that the authors have tended 
to convey the impression that they are de- 
scribing a working system rather than one 
under development. For those interested in 
magnetic core storage techniques this paper 
will be-a required addition to their library. 

Raymond Stuart-Williams 


57-186 
Megabit Memory—R. A. Tracy. (Proc. 
- Eastern Joint Comp. Conf., December 10-12, 
1956, New York, pp. 104-106; 1957.) An 
interesting new method of making magnetic 
cores and memory planes is described in this 
short but well-illustrated article. Thin strips 
of 4-79 molybdenum permalloy (4 inch 
wideX3 inch long by $ mil thick) are wound 
directly on the selecting and sense wires 
to form a core position. This wrapping 
technique provides a minimum inside diam- 
eter and thus keeps the driving currents 
small. The performance of sample wound 
cores is compared with an S-3 ferrite core 
in graphs of signal-to-noise ratio and tem- 
perature dependence and in photographs of 
core output waveforms. The signal sutputs 
were 30-40-mv peak amplitude. The higher 
Curie temperature metal tape core was less 
temperature sensitive and gave a greater 
signal-to-noise ratio then the ferrite. The 
evolution of a completely wired 2020 
memory plane is described and illustrated 
with photographs. All wires (Y, Y, Z, and 
SENSE) are first assembled in their proper 
relationship and then the cores are wrapped 
on at the intersections by a semiautomatic 
machine. Further cletails of this machine 
Wrapping process and how uniformity is 
maintained are not given in this article. 
Several advantages over ferrites are claimed 
and discussed that include: 1) less tempera- 
ture sensitivity, 2) negligible strain sensi- 
tivity, 3) major hysteresis loop operation, 
4) very good uniformity of tape from the 
same melt. Large differences exist between 
melts, but 10-20 million cores can be made 
from one 50-pound melt, 5) useful operation 
over a wider current range, 6) a factor of 
10-cost advantage over ferrite cores. The 
author claims that these improvements now 
make million-bit matrix memories quite 
feasible. An interesting question and answer 
discussion section follows the article. 
W. W. Lawrence, Jr. 


681.142 :621.375.3 57-187 
A Magnetic-Amplifier Switching Matrix 
—D. Katz. (Commun. and Electronics, no. 24, 
pp. 236-241; May, 1956.) The circuit de- 
scribed translates a four-digit binary code 
into an electronically displayed decimal 
output, 
Courtesy of Proc. IRE 
and Wireless Engineer 
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57-188 

A Large-Capacity Drum-File Memory 
System—H. F. Welsh and V. J. Porter. 
(Proc. Eastern Joint Comp. Conf., December 
10-12, 1956, New York, pp. 136-138; 
1957.) This paper describes the 2,500,000- 
character flying head drum memory designed 
for Univac-Larc, which provides a maximum 
access time of 2.6 seconds. The presentation 
is essentially a product description, outlining 
the basic design philosophy in terms of over- 
all system requirements and the design steps 
leading to the present configuration, which 
includes an air floating head which is posi- 
tional by a linear indexing mechanism over 
100-information bands. The authors reveal 
a clear understanding of the formidable 
problem involved in providing a high-speed 
electromechanical random access memory 
and indicate significant progress towards a 
solution, which at present seems a better 
commercial answer than the use of multiple 
magnetic disks or hydraulically lubricated 
drums. The claim that close mechanical 
precision is eliminated by the air floating 
design scems somewhat misleading in the 
case of this particular drum since the com- 
plexity of the head positioning mechanism 
more than compensates. Also, in adopting 
such a head design, one trades temperature 
differential and control problems for new 
problems of extreme cleanliness and surface 
smoothness. In view of the intimate rela- 
tionship among drum, head, and reading 
system the paper might well have discussed 


‘the associated recording system more ex- 


teusively, particularly since excellent density 
and clock rates were achieved (400 bits 
/inch). 

D. N. MacDonald 


681.142 :621.395.625.3 57-189 
High-Density [magnetic-] Tape Record- 
ing for Digital Computers—(Elec. Mfz., vol. 
56, pp. 153, 290; November, 1955.) Digital 
pulse densities of up to 700/inch can be 
recorded by the method described which was 
developed by the National Bureau of Stand- 
ards for use with the SEAC computer. 
Courtesy of Proc. IRI¢ 
and Wireless Engineer 


621.318.57 :621.387 57-190 
Dekatron Drive Circuit and Application 
—M. Graham, W. A. Higinbotham, and 
S. Rankowitz. (Rev. Sct. Instr., vol. 27, pp. 
1059-1061; December, 1956.) “A reliable, 
one-tube drive circuit for decade glow- 
transfer counter tubes is described. Applica- 
tion of the circuit is illustrated in a ten- 
channel glow-tube register with automatic 
electric-typewriter readout.” 
Courtesy of Proc. IRE 
and Wireless Engineer 


681.142 57-191 
An Accumulator Unit for a Dekatron 
Calculator—R. Townsend and K. Camm. 
(Electronic Eng., London, vol. 29, pp. 58-64; 
February, 1957.) Decimal numbers may be 
added or subtracted, the sum being returned 
to the main store multiplied or divided by 
ten or unity. The calculator has a punched- 
card input. 
Courtesy of Proc. IRE 
and Wireless Engineer 
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621.385.832.032 57-192 
The Technology of Electrostatic-Storage 
Cathode-Ray Tubes—P. Choffart. (Onde 
élect., vol. 36, pp. 815-821; October, 1956.) 
An outline of some of the manufacturing 
problems. 
Courtesy of Proc. IRE 
and Wireless Engincer 


57-193 
Polarization Reversal and Switching in 
Guanidinium Aluminum Sulfate Hexa- 


hydrate Single Crystals—H. H. Wieder. 
(Proc. IRE, vol. 45, pp. 1094-1099; August, 
1957.) The above ferroelectric material, 
GASH, is compared to the more familiar 
barium titanite. For both crystals no 
threshold field and hence no true coercivity 
was found experimentally. The upper prac- 
tical limit of driving frequency appears to 
be 25 to S0 ke for both crystals. GASH has 
lower hysteresis losses, very much lower 
dielectric constant, and lower electro- 
mechanical activity than barium titanite. 
The author feels that GASH has many 
properties of interest but that a great deal 
of additional work is required to establish 
its proper place as a reliable solid-state 
electronic component. 


A. D. Scarbrough 


‘57-104 
Fast Switching by Use of Avalanche 
Phenomena in Junction Diodes—B. Salz- 
berg and E. W. Sard. (Proc. IRE, vol. 45, 
pp. 1149-1150; August, 1957.) In this letter, 
the authors discuss the fast switching speeds 
in both “forward” and reverse directions 
which can be obtained by operating around 
the avalanche breakdown voltage. A number 
of oscillograms are presented which show 
extremely fast forward switching (transient 
undetectable with a Velktronix 545). 
A. D. Scarbrough 


681.142 :621.314.7 57-195 
The Junction Transistor as a Computing 
Element—E. Wolfendale, L. P. Morgan, and 
W. L. Stephenson. (Electronic Eng., vol. 29, 
pp. 2-7, 83-87, 136-139; January, March, 
1957.) The small-signal and transient char- 
acteristics of the transistor and their appli- 
cation to the design of basic circuits are 
described and examples of computer ele- 
ments using transistors are given. 
Courtesy of Proc. IRE 
and Wireless Engincer 


57-196 

High-Speed Transistor Computer Cir- 
cuit Design—R. A. Henle. (Proc. Eastern 
Joint Comp. Conf., December 10-12, 1956, 
New York, pp. 64-66; 1957.) This paper 
describes a novel switching-circuit philoso- 
phy particularly suited for the class of drift 
transistors with which appreciable improve- 
ment in switching speed may be obtained 
by using high currents and collector voltages, 
It is shown that high-speed operation de- 
mands small signal-voltage swings. These 
requirements are met by controlled-current 
switching, wherein a fixed current source 
supplies two or more transistor emitters in 
parallel. This entire current flows through 
that transister with the most forward base 
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bias. Because emitter current is controlled, 
the collector current is limited, so that 
collector-voltage swings can be maintained 
within any desired range by suitable circuit 
impedances. This technique avoids low 
collector voltages and hole-storage delay at 
the expense of high transistor dissipation. 
This price may be great, for it is probable 
that the ultimate transistor reliability 
needed for large-scale computers will be 
achieved only with very small dissipation. 
Furthermore, this- circuitry may require 
complementary p-n-p and n-p-n transistors 
without which some basic circuits, such as 
flip-flops, may be relatively complex. (“Mil- 
limicrosecond Transistor Current Switching 
Circuits’—Hannon S. Yourke, presented at 
the 1957 Transistor and Solid State Circuits 
Conference, Philadelphia, Pa., February, 
1957.) Nevertheless, until devices yielding 
high speed in saturating circuits become 
available, controlled-current — switching 
should fill many needs in the high-speed 
computer frontier. 


J. B. Angell 


57-197 

High-Temperature Silicon-Transistor 
Computer Circuits—J. B. Angell. (Proc. 
Eastern Joint Comp. Conf., December 10-12, 
1956, New York, pp. 54-57; 1957.) The 
switching characteristics of the Philco 
T-1159 silicon surface alloy transistor are 
described. as functions of temperature. It is 
shown that direct coupled transistor circuits 
can be made to operate satisfactorily from 
~—50°C to +140°C with these transistors. 
Switching times of 0.1 to 0.2 usec at room 
temperature and of about 0.5 ysec at 
+140°C are obtained in such circuits. The 
statements concerning the temperature 
stability of direct coupled circuits as com- 
pared to other types of circuitry are debata- 
ble, especially when the distribution of 
base drive currents due to different input 


impedances are considered. 
W. B. Cagle 


57-198 
High Temperature Printed Circuitry— 
Quarterly Report No. 2, Third Series— 
ieee Bowman, ©. Hiaek. Wallsins, eb al: 
(Quart. Rep. Computer Components Fellow- 
ship Mellon Inst., 36-+ix pp.; January 1, 
1957—March 31, 1957.) Measurements of 
resistance, capacitance, and dissipation 
factor as a function of temperature over the 
range 28° to 260°C have been taken of some 
resistor and capacitor films. The film com- 
ponents tested were prepared by vacuum 
evaporation techniques and by the screen 
processing of enamel materials. Cycling the 
films in this temperature range showed that 
the components stabilized during the first 
cycle. Details are presented for the design 
of 4 continuous kiln for the firing of enamel 
coniponents constructed by this laboratory. 
Of thirteen different solders tested for form- 
ing high temperature connections, only 
three were found to retain suitable mechani- 
cal and electrical properties after aging in an 
oven at 250°C for five weeks. Consideration 
is given to the development of pressure type 
connections for high temperature applica- 
tron, 


C. AH. tf. Wilkins 
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57-199 
High Temperature Printed Circuitry— 
Quarterly Report No. 3, Third Series— 
G. H. Young, C. H. T. Wilkins, et a/. (Quart. 
Rep. Computer Components Fellowship Mel- 
lon Inst., 21+-vi pp.; April 1, 1957-June 30, 
1957.) Enamel resistance films of bismuth 
fluxed and of glass fluxed gold have been 
found to have positive temperature coeff- 
cients between 25° and 500°C. The addition 
of TiOz to glass enamels has been demon- 
strated to increase the capacitance of printed 
capacitors by increasing the dielectric con- 
stant. A new extruded alumina sheet is 
being tested as a substrate material for use 
with both the screen printed enamels and 
the vacuum evaporated circuits. A pre- 
liminary study of the subzero properties of 
modular printed circuits has been under- 
taken. The method used and results obtained 
with vacuum evaporated gold-palladium 
resistance films in the range from —196°C 
to room temperature are presented. In the 
high temperature testing program, problems 
encountered in devising suitable soldering 
joints are discussed and present soldering 
procedure described. Difficulties encountered 
in using the large testing oven have ied to 
the development of simplified heating units 
capable of testing one sample at a time up 
to temperatures as high as 850°C. 
C. H. T. Wilkins 


DIGITAL SYSTEMS 
RESEARCH 


57-200 
Arithmetical Analysis of Digital Com- 
puting Nets—Richard C. Jeffrey. (J. Assoc. 
Comp. Mach., vol. 3, pp. 360-375; October, 
1956.) The author indicates several possible 
numerical interpretations of a sequence of 
binary digits, where the sequence may be 
temporal or spatial. They differ according 
to the assumed location of the radix point 
and the representation of negative numbers. 
He then shows how the Boolean algebraic 
description of a computing net can be 
translated into an arithmetic description of 
its operation in terms of whichever inter- 
pretation is of interest. 
A. S. Householder 
Courtesy of Mathematical Reviews 


57-201 

An Algorithm for Determining Minimal 
Representations of a Logic Function— 
B. Harris. (IRE Trawns., vol. EC-6, pp. 103- 
108; June, 1957.) For each logic function, or 
Boolean algebraic expression, there corre- 
sponds an appropriate computer circuit. 
However, the minimization of the appear- 
ances of the Boolean variables does not 
necessarily lead to the most economical 
circuit. A general approach to the problem 
therefore requires the development of tech- 
niques for the simple and rapid generation 
of a variety of near-minimal forms. This 
paper describes such a method for construct- 
ing the minimal representations of a logic 
function given as a truth table or in one of 
its canonical forms. The minimal representa- 
tions achieved are cither sums of products, 
or products of sums, such that no term 
contains superfluous variables and such 
that no term is superfluous. The utility of 
the method lies in the conciseness of nota- 
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tion, which permits the handling of a large 

number of variables and simplifies the 
process for machine computation. 

Author’s Summary 

Courtesy of Proc. IRE 


57-202 
Minimization of the Partially-Developed 
Transfer Tree—M. P. Marcus. (IRE 


TRANS., vol. EC-6, pp. 92-95; June, 1957.) 
A transfer tree is a particular type of multi- 
terminal network having a single input 
which may be connected to any one of a 
number of outputs. An -relay transfer tree 
is partially developed if it has less than the 
2” possible output terminals. Rearrange- 
ment of a partially developed tree can lead 
to a reduction in the total number of trans- 
fers required. This paper presents a method 
of obtaining a required partially developed 
transfer tree with the minimum number of 
transfers. 
Author’s Summary 
Courtesy of Proc. IRE 


57-203 
A Time-Sequential Tabular Analysis of 
Flip-Flop Logical Operation—G. W. Arant. 
(IRE Trans., vol. EC-6, pp. 72-74; June, 
1957.) In examining flip-flop response the 
principal concern of the logical designer is to 
find what input signals must be applied to 
the flip-flop in order to produce the output 
conditions that are desired. Equation meth- 
ods of analysis and a time-sequential tabular 
method of analysis are described, and some 
advantages of the tabular method are 
pointed out. 
Author’s Summary 
Courtesy of Proc. IRE 


57-204 

Synchronization of a Magnetic Com< 
puter—J. Kielsohn and G. Smoliar. (Proc. 
Eastern Joint Comp. Conf., December 10-12, 
1956, New York, pp. 90-93; 1957.) The 
digital signals which circulate in a high- 
speed synchronous computer are of short 
duration, have precisely determined shape 
and magnitude, and occur at precise times. 
Signals from computer input equipment are 
usually of much greater duration with un- 
determined shape and have random occur- 
rence times as far as the computer is 
concerned. Therefore it is always necessary 
to have a device which will inject into the 
computer a sequence of properly shaped and 
timed signals which bears a one-to-one 
informational correspondence to the se- 
quence of signals from the input equipment. 
Such a device is called a synchronizer. The 
authors describe a synchronizer employing 
diode gating and high-speed magnetic 
digital amplifiers of the series, or Ramey 
type, developed by Remington-Rand,. The 
Cambridge computer, for which this syn- 
chronizer was designed, makes use of the 
same kind of magnetic amplifier. A brief 
description of the amplifier is given. 
Logically, the synchronizer is of familiar 
design. It consists of two dynamic flip-flops 
in cascade. The first flip-flop, when triggered 
by the unsynchronized signal, sends a series 
of pulses to the and-gate input of the second 
flip-flop. Coincident with the next computer 
timing pulse introduced at the other input 
to the and-gate, one of the pulses emitted 
by the first flip-flop attempts to trigger the 
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second. A full timing pulse period (two word 
lengths) is available for the second flip-flop 
to respond, or to fail to respond, before its 
output is sampled at the time of arrival of 
the next ensuing timing pulse at the output 
and-gate. Because the response time of the 
flip-flop is considerably less than the timing 
period, a fully shaped and synchronized 
output signal is guaranteed. A description 
of the operation of the magnetic shift register 
stage of the computer is also given in some 
detail. 

E. W. Hogue 


57-205 
Designing for Reliability—Norman H. 
Taylor. (Proc. IRE, vol. 45, pp. 811-822; 
June, 1957.) A circuit design philosophy 
developed at Lincoln Laboratories is pre- 
sented which tends to produce reliable cir- 
cuit designs. The concept of marginal 
checking is discussed and illustrated by a 
great number of plots for a high-speed 
vacuum tube flip-flop. A number of useful 
rules of thumb for component derating is 
presented and the reliability of various 
components is discussed. Unfortunately, 
nothing helpful is said about transistors. 
A. D. Scarbrough 


57-206 
Topology of Switching Elements vs 
Reliability—|. P. Lipp, (IRE TRAns., 
PGRQC-10,. pp. 21-33; June, 1957.) The 
topology of switching elements is explored 
from a probability or reliability standpoint. 
Arrays are indicated which offer a greater 
degree of reliability than any one of their 
elements, hence disproving the belief that 
circuit complexity necessarily lessens reli- 
ability. A further result is production of a 
simple mathematical technique adaptable to 
the solution of complex circuits in terms of 
reliability. Such “circuits” may actually con- 
sist of redundant communication paths as a 
more general interpretation of switching 
elements. 
Author’s Summary 
Courtesy of Proc. IRE 


57-207 
Improvement of Binary Transmission by 
Null-Zone Reception—-F. J. Bloom, S. S. L. 
Chang, B. Harris, A. Hauptschein, and 
K. C. Morgan. (Proc. IRE, vol. 45, pp. 
963-975; July, 1957.) This paper discusses 
the improvement in transmission of binary 
information in the presence of noise that can 
be effected by decoding into 3 or 4 states 
instead of the usual two. In the systems 
discussed the receiver would indicate state 
1, state 2, or undecided (for single null re- 
ception); or state 1, state 2, undecided but 
probably state 1, or undecided but probably 
state 2 (for double null-zone reception). The 
tse of null-zone reception permits a single 
parity check which normally gives single 
error detection to have substantially single 
error correcting properties since if the single 
error is registered as a null, the parity check 
tn be used to correct the error. Tt is shown 
that single null reception greatly invproves 
the effectiveness of the single parity check. 
Double null reception provides only a small 
“(ditional improvement in per-unit equivo- 
cation but permits much less critical adjust- 
ment of the null level. 
A. D. Scarbrough 
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57-208 
Die Mathematischen Grundlagen fiir 
die Organisation der Elektronischen 
Rechenmaschine der Eidgendssischen Tech- 
nischen Hochschule—John Robert Stock. 
(Mitt. Inst. Angew. Math. Ziirich, no. 6, 73 
pp.; 1956.) This gives an account of the 
logical design of the ERMETH. This is a 
one-address serial floating decimal machine 
with a magnetic drum memory of 10,000 
words of 16 digits. A floating decimal number 
is of the form @:-10', where —20056S5199 
and |a|<10; alternatively a word can 
represent a fixed decimal number with 14 
decimal digits. There is another decimal 
digit in the word, which serves two purposes 
as a parity check (mod 3) and as an aid 
(Q sign) for a programming system devised 
by Rutishauser [same J/itt. no. 3 (1952); 
MR 15, 64]. An instruction consists of 7 
d.d.: four for the address, two for the opera- 
tion, and a third indicating one of nine index 
registers. The machine arithmetically is 
roughly a 10-m one and has an average 
access time of about 5 m. The normal input 
is by means of punched cards. The first 
introductory chapter includes a list of the 
instructions. The second describes the 
actual form of the operations. The third 
chapter discusses the principles of realiza- 
tion of the machine. The fourth chapter is 
concerned with the logical design of the vari- 
ous ,components. There are appendixes 
which give block diagrams of the arithmetic 
unit and the control. There are a few simple 
examples of programming for the machine 
which illustrate the use of index register 
and the method of reading in a program; 
there is also a subroutine for square root and 
a program for the solution of a system of 
equations by an elimination process, which 
uses the index register and the Q sign. 
J. Todd 


Courtesy of Mathematical Reviews 


57-209 
The Transac S-1000 Computer—J. L. 
Maddox, J. B. O’Toole, and S. Y. Wong. 
(Proc. Eastern Joint Comp. Conf., Decem- 
ber 10-12, 1956, New York, pp. 13-16, 
1957.) This paper presents a description 
of a scientific computer designed at the 
Philco Corporation. The machine uses the 
well-known direct-coupled  surface-barrier 
transistor computer circuits reported on at 
the IRE Convention in March, 1955, by 
R: H. Beter, W. E. Bradley, R. B. Brown, 
and M. Rubinoff. Because of the packaging 
with printed circuits and the small volume 
associated with the transistor circuits the 
machine occupies only the usable space in a 
standard office desk. The memory, obtained 
from the ERA Division of Remington Rand 
UNIVAC, has 4096 words of 36 bits each 
contained in magnetic cores. The machine 
requires 1.2 kw of power. Much of the paper 
is concerned with the logical system of the 
machine and a list of instructions which 
has some similarity to the instruction list 
fora UNIVAC 1103A. The word length of 
36 bits when used as an instruction contains 
six bits for the command, three bits for each 
of two address modifiers, and two 12-bit 
addresses. 
Arithmetic is performed in the 1’s com- 
plement parallel binary system by a 72- 


stage subtractor working with a 72-bit ac- 
cumulator. A repeat instruction is avail- 
able. The construction of the machine was 
80 per cent finished when the paper was 
written. The machine may be considered to 
be a relatively straightforward application 
of the surface-barrier transistor circuits. 
R. E. Meagher 


57-210 

TX-O, A Transistor Computer with a 
256 by 256 Memory—J. L. Mitchell and 
K. H. Olsen. (Proc. Eastern Joint Conf., 
December 10-12, 1956, New York, pp. 
93-101; 1957.) This two-part paper describes 
the TX-O magnetic core memory and a 5- 
megapulse transistor logic circuitry. The 
TX-O memory provides storage for 65,536 
19-bit words. It is a parallel, conventional 
coincident-current magnetic-core unit hav- 
ing a 7.0-wsec read-write cycle. Random 
access to the memory is provided by two 
256-position tape-wound, magnetic core 
switches. A block diagram of the memory 
system and a timing diagram are given. The 
memory and switch cores, the memory ar- 
ray windings and organization, and the de- 
sign and operation of the access switches are 
fully described. The circuits of the memory 
system, which use 425 dual triodes and 625 
transistors, are given with schematics. The 
system margins are shown. The TX-O 
utilizes transistor logic in two basic circuit 
configurations; the saturated inverter and 
the saturated emitter follower. Reverse 
base bias and speed-up condensers are used. 
Logic is performed by series and parallel 
combinations of these configurations. A 
quite complicated flip-flop circuit is used in 
an effort to have one standard flip-flop. Cir- 
cuit schematics, waveforms, margins, and 
packaging are shown. While the description 
of the memory is quite complete, the descrip- 
tion of the rest of the computer is restricted 
to the circuitry. Since neither the logical 
design nor the system characteristics are 
discussed the title of this paper is misleading. 


J. A. Githens 


681.142 57-211 
The ‘‘Bizmac” Digital Data Processing 
System—J. C. Hammerton. (Electronic 
Eng., vol. 29, pp. 174-180; April, 1957.) 
As used at the central agency for American 
military supply depots. 
Courtesy of Proc IRE 
and Wireless Engineer 


57-212 

Design Objectives for the IBM Stretch 
Computer—S. W. Dunwell. (Proc. Eastern 
Joint Comp. Conf., December 10-12, 1956, 
New York, pp. 20-22; 1957.) This paper is 
a progress report on a very impressive com- 
puter system in its research phase. 
STRETCH is a multiple arithmetic unit 
computer system in the very large, ten- 
megapulse class. It is intended for both 
scientific and business applications. It con- 
sists of three separate stored program com-~ 
puter sections especially tailored to the re- 
quirements of handling 1) input-output, 
2) data format editing, and 3) computa- 
tion. All communicate with one multi- 
section memory and all operate simulta- 
neously to execute about one million instruc- 
tions per second. Computer 1 will service 
individual input-output devices of a very 
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wide class. Computer 2 is to be serial and 
handles variable word length entries in 
either binary or alpha-numeric notation. 
Basic speeds of this section are to be about 
3-3 usec for addition and 5-15 psec for mul- 
tiplication. Computer 3 is to be parallel and 
will perform floating point addition in 0,6 
usec and multiplication in 1.2 psec. House- 
keeping operations for this computer section 
proceed concurrently with computation. The 
instruction structure will allow for one mil- 
lion words of random access memory. Ex- 
ternal memory can be expanded to one 
hundred million words. The ten-megapulse 
solid-state logical elements and ferrite mem- 
ories are described in other papers in the 
referenced proceedings. The fast memory 
will respond in 0.2 usec. The project is re- 
ported to be substantially on schedule. The 
first system will go to Los Alamos. 

William F. Gunning 


57-213 
A Transistorized Transcribing Card 
wBunch——C. i iCole, jr, Ko i Chien, and 
C. H. Propster, Jr. (Proc. Eastern Joint 
Comp. Conf., December 10-12, 1956, New 
York, pp. 80-83; 1957.) This paper describes 
a card punch used for the conversion of large 
volumes of data stored on magnetic tape in 
the Bizmac code into punched cards. The 
unit will transcribe information at a rate of 
150 cards per minute and provides accuracy- 
checking features. Three cards’ worth of 
information are in process at any one time — 
onesct of information being read into storage, 
the second set being punched on the card, 
and the third set being read back from a 
punched card and compared with the stor- 
age. The punching and card handling tech- 
niques appear to be conventional in nature. 
The electronics of the unit are completely 
transistorized, with two-input resistor logic 
being used as the standard logical eiement. 
Combinations of these elements, which are 
basically a two input “or” gate, are utilized 
to generate an “and” gate (by inverting the 
inputs to the standard “or” gate), a “not” 
gate, an inverter, a power amplifier, and a 
flip-flop. 
M. J. Mendelson 


57-214 

Automatic Input for Business Data- 
Processing Systems—K. R. Eldredge, F. J. 
Kamphoefner, and P.'H. Wendt. (Proc. 
Eastern Joint Comp. Conf., December 10-12, 
1956, New York, pp. 69-72; 1957.) The 
paper presents a semitechnical description 
of a machine which reads printed matter 
(presently, only numerical digits plus a few 
symbols) printed with magnetic ink, de- 
veloped by the Stanford Research Institute. 
Described are three interrelated problems 
connected with the design of the equipment: 
1) electrical and magnetic techniques for 
reliable reading of characters; 2) develop- 
ment of suitable magnetic ink for printing; 
3) development of an electromechanical ma- 
chine for reliably handling documents. The 
reading circuits employ a magnetic head 
with a gap width approximating the width 
of characters, a tapped electrical delay line 
allowing simultaneous examination of the 
electrical signals produced by the head over 
the period of time during which the charac- 
ter passes the reading head, and a correlation 
analysis network for recognizing the par- 
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ticular distribution of signals in the delay 
line. The printing of the information on the 
document is merely a good quality printing 
process employing a specific type of ink 
made up of a magnetic pigment. The docu- 
ment handling techniques employ a pneu- 
matic system and are capable of handling 
documents which have been mutilated to 
quite an extent. The equipment has been 
employed principally inreading identification 
numbers from bank checks. Figures on re- 
jection and error rate when reading these 
documents are presented. Following the in- 
teresting but brief presentation, a list of 
questions and answers serve to give further 
insight into this outstanding development. 


L. S. Michels 


57-215 
A Tank Farm Data Reduction System— 
D. J. Gimpel. (RE Trans., vol. PT-2, pp. 
94-100; April, 1957.) A tank farm data re- 
duction system has been developed for the 
new Tidewater Oil Company installation in 
Delaware by the Armour Research Founda- 
tion and Panellit, Inc. The function of the 
unit is to secure the temperature corrected 
volume of fluid in each of the approximately 
100 tanks in the field. The inputs to the sys- 
tem are the fluid height and average tank 
temperature. Fluid volume is tabulated digi- 
tally as a function of the height on magnetic 
tape. The system automatically searches the 
tape for the indicated volume and multiplies 
the number by the temperature correction 
factor. This paper describes the operation of 
the multiplier, the tape search elements, and 
the sensing instruments employed in the 
field. The factors governing the selection of 
the specific elements in the storage and com- 
puting system are also discussed. 
Author’s Summary 
Courtesy of Proc. IRE 


681.142 57-216 
A High-Speed Data Processing System 
—M.L. Klein, R. B. Rush, and H. C. Mor- 
gan. (Electronic Eng., vol. 29, pp. 158-163; 
April, 1957.) Data from 20-100 independent 
channels, at a full-scale level of 100 my, are 
taken at a rate of 100 kc and recorded 
digitally on a magnetic tape together with 
the source information and time. 
Courtesy of Proc. IRE 
and Wireless Engineer 


UTILIZATION OF DIGITAL 
EQUIPMENT 


57-217 
Chemical Computing—(Chem. Eng. 
News, vol. 34, pp. 3542-3545; July 23, 1956.) 
This article gives a brief survey of current 
computer applications in engineering, data 
reduction, research, and personnel problems. 
It is of interest primarily to those unfamiliar 
with the computer field. 
Courtesy of Data Processing Group 
General Motors Research Staff 


57-218 

Digital Computer as an Aid to the Elec- 
trical Design Engineer—B. Birtwistle and 
Beryl M. Dent. (Engineer, vol. 201, pp. 440- 
442; May 4, 1956.) Two applications are 
given of the digital computer: 1) impulse 
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voltage distribution in transformer winding, 
and 2) starting torque calculations for sali- 

ent pole synchronous motors. 
Courtesy of Data Processing Group 
General Motors Research Staff 


681.142 :621.316.11 57-219 

Automatic Network Analysis with a Digi- 

tal Computation System—S. Y. Wong and 

M. Kochen. (Comin. and Electronics, no. 

23, pp. 172-175; May, 1956. Discussion, p. 

176.) A general outline of the procedure and 
a method of solution are given. 

Courtesy of Proc. IRE 

and Wireless Engineer 


57-220 
Fitting the Digital Computer into Process 
Control—Montgomery Phister, Jr., and 
Eugene M. Grabbe. (Control Eng., vol. 4, 
pp. 129-136; June, 1957.) This article de- 
scribes some of the techniques and problems 
involved in applying the digital computer 
to process control. The authors point out 
that in addition to knowledge of computer 
techniques a rather broad knowledge of the 
process itself is required. The authors de- 
scribe a five step procedure of the systems 
development: 1) theoretical analysis of 
process, 2) correlation of variables, 3) in- 
vention, 4) specification of system com- 
ponents, 5) system operation. They give 
a design case history of a chemical process 
including a reactor, a heat exchanger, a 
catalyst separator, and a  fractionating 
tower, and describe in some detail how the 
operational system can be synthesized. Also 
included is the development of the appropri- 
ate variable relationships from empirical 
data, together with diagrams of optimal 
operating loci. They describe how the com- 
puter program might be designed to ac- 
complish the necessary control and how the 
digital control would be connected with the 
system. 
Neal J. Dean 


57-221 
Tape Identification and Rerun Pro- 
cedures for Tape Data Processing Systems 
—Leonard Eallson. (Computers and Auto- 
mation, vol. 5, pp. 12-13; April, 1956.) This 
is a tutorial article which briefly describes 
widely used techniques for detection of cer- 
tain tape handling errors and the recovery 
of the required information where reruns are 
necessary. The errors are detected by use of 
an information block on the tape containing 
identifying information, block counts, and 
other control information. The information 
for reruns is provided by a memory dump 
at the end of each tape output. This article 
is of no value to anyone familiar with rou- 
tine data-processing procedures. 
Gordon E. Morrison 


57-222 

The Digital Approximation of Contours— 
Robert M. Mason. (J. Assoc. Comp. Mach., 
vol. 3, pp. 355-359; October, 1956.) The 


problem is to program a computer to trace 
a contour f(x, y)=& on a grid. Suppose, for 
definiteness, f<k inside the contour line. 


The sign of f—k is examined at a sequence 
of points differing by Ax or by Ay. The se- 
quence is rectilinear until f —& changes sign, 
after which the direction rotates by 7/2 if 
f —k has become positive, —7/2 if f—k has 
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become negative. For getting onto the con- 

tour initially the machine secks out a path 
of steepest ascent or descent. 

A. S. Houscholder 

Courtesy of Mathematical Reviews 


57-223 
Some Experiences on the Manchester 
Computer with the Simplex Method—-D. G. 
Prinz. (Conference on Linear Programming, 
Ferranti Ltd., London, pp. 73-89; May, 
1954. Discussion, pp. 90-91.) This paper 
presents an account of the coding of the 
simplex method on the Manchester com- 
puter by the author. It is free from fussy de- 
tails pertinent to a particular machine, but 
it does not shirk explanation of the numeri- 
cal difficulties encountered by round of er- 
rors and how they were overcome. Mathe- 
matically, the most interesting point is the 
attempt to show some justification for the 
fact that the simplex method seems to take 
a number of iterations much lower than the 
conceivable (though never properly esti- 
mated) upper bound. The author arrives at 
the average result log, Cy, where a is the 
number of rows, 7 the number of columns 
of the tableau. The reasoning is heuristic and 
not particularly convincing, but it is, so far 
as the reviewer knows, the first brave effort 
to estimate the number of iterations. 
A. J. Hoffman 
Courtesy of Mathematical Reviews 


681.142:51 57-224 
Introduction to a Theory of Ensembles 
Based on the Prime Numbers Assimilable 
by Electronic Computers—S. Sabliet. (C, R. 
Acad. Sci., Paris, vol. 244, pp. 35-38; Janu- 
ary 2, 1957.) Computers of “ordinator” type 
can be made to provide algebraic solutions of 
equations by coding the functions to be op- 
erated on as prime numbers. 
Courtesy of Proc. IRE 
and Wireless Engineer 


BOOK REVIEWS 


57-225 

Elektronische Rechenmaschinen und In- 
formationsverarbeitung—Edited by A. Wal- 
ther. (Verlag Friedr. Vieweg & Sohn, 
Braunschweig, Germany, 229 pp.-fviil, 
illus.; 1956.) Like many other proceedings of 
computer conferences, this volume covers a 
wide range of subjects, from the design to 
the application of digital computers. This 
one, however, has a decidedly cosmopolitan 
flavor. The authors come from eleven coun- 
tries and their papers are printed in three 
different languages. The title Electronic 
Digital Computing and Inforniation Process- 
ing, was that of a conference held in Darm- 
stadt, Germany, in October, 1955. The or- 
ganizer of the conference, Prof, A. Walther 
of Darmstadt, is also the editor of this book 
containing the sixty-four papers delivered at 
the conference. He has ably performed a 
most difficult job. The authors represent 
most of the major computing centers in 
Western Europe, as well as the United 
States and three eastern countries (the 
USSR, Czechoslovakia, and Eastern Ger- 
many). Prof. Walther has collected their 
papers into a readable and interesting vol- 
ume, Forty-five of the papers are in German, 
seventeen in English, and two in French. 
English-speaking readers will find that the 
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linguistic hurdles have been lowered by the 
inclusion of English abstracts for papers not 
in English. A curious sidelight is the fact 
that the two Soviet papers appear in German 
while the three Czech papers are printed in 
English. The book records an important 
conference and is thus of historical valuc. 
Many of the papers describe machines and 
techniques employed at various European 
computing centers. Others present research 
work of more general interest. The volume 
contains some unusual items, such as a 
demonstration that all arithmetic and logi- 
cal computer operations can be built up 
from just one essential instruction contain 
ing only an address. The technical material 
will be largely familiar to those who have 
followed developments in the U. S. What is, 
perhaps, most interesting to such readers is 
the account of the many computer projects 
which are carried on in Europe, and the re- 
sourcefulness displayed in the face of limited 
budgets. 
Werner Buchholz 
Courtesy of Proc. IRE 


- 57-226 

Elektronnye Vychislitel’nye Mashiny i 
Obrabotka Informatsii (Electronic Comput- 
ing Machines and Information Processes) — 
S. A. Lebedev. (Akad. Nauk S.S.S.R., 
Moscow, 47 pp.; 1956.) This pamphlet is one 
of the series of educational booklets issued 
by the Academy for popular consumption. 
Tts author is the Academician in charge of 
BESM (this name is made up of the initials 
of the four Russian words standing for High- 
speed Electronic Computing Machine) and 
principal compiler of the Academy’s Index 
of Mathematical Tables, Moscow, 1956 
(Rev. 49, pp. 104-106, this issue). The 
booklet is couched in delightfully lucid lan- 
guage and is profusely illustrated with dia- 
grams, photographs, and sample programs. 
The modest price, equivalent to about 15 
cents, makes it available to anying having 
an interest in the subject, and its high quali- 
ty guarantees complete satisfaction to each 
purchaser. The introduction points out the 
benefits of high-speed computation and em- 
phasizes, in particular, the achievements 
made possible through the use of BESM not 
only in numerical computation but also in 
the initial—and quite spectacular—at- 
tempts to translate scientific English into 
Russian. The first chapter is devoted to a 
careful explanation of how electronic ma- 
chines function. In the second, we find an 
excellent discussion of the binary vs the 
decimal systems with a complete explana- 
tion of the engineering principles which are 
used to represent the former system. This is 
followed by some very timely observations 
on the pros and cons of floating radix point 
representation of numbers. The last chapter 
deals with a detailed description of BESM. 
Since the features of this machine are al- 
ready familiar to the readers of this maga- 
vine, there is no need to repeat the facts 


here. 


The High-Speed Electronic Calculating 
Machine of the Academy of Sciences of the 
U.S.S.R.—S. A. Lebedev, J. Assoc. Comp. 
Mach., vol. 3, pp. 129-133; July, 1956. See 
Rev. 57-24; March, 1957.) The pamphlet is 
enriched by an excellent appendix, which 
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can easily serve as a textbook for tyrosin 
coding. It outlines the basic principles of 
program preparation, and illustrates them 
by actual examples. It touches upon the use 
of subroutines and the choice of suitable 
checks to insure the accuracy of the com- 
putations. An altogether charming booklet 
for 65 kopacks! 
Ida Rhodes 
Courtesy of Mathematical Tables 
and. Other Aids to Computation 


57-227 
Digital Computer Programming—D. D. 
McCracken. (John Wiley & Sons, Inc., New 
York, N. Y., 218 pp.+30 appendices-++5 
index+-vii; 1957.) Primarily an introductory 
text on digital computer programming, this 
book is also a survey of programming tech- 
niques. It is the first general treatment of 
programming without excessive concern for 
a particular aspects of the subject or for a 
particular computer. A successful exposition 
of computer programming must necessarily 
refer to a specific machine. The author has 
chosen to use a hypothetical computer, 
TYDAC. Lile the alternative solution of 
using an existing computer, the decision to 
use TYDAC is not ideal but it is, however, 
skillfully implemented. The book remains a 
book on programming and not on program- 
ming TYDAC. The reader with no computer 
at his disposal will find TYDAC easy to 
understand and to program, although an 
element of reality is absent because YDAC, 
unlike most real computers, has no peculi- 
arities or irregularities. The reader witha 
computer will find that the illustrative pro- 
grams do not depend on the detailed charac- 
teristics of TYDAC, that the exercises are 
readily adaptable to other machines, and 
that chapters not applicable to his machine 
may be omitted without loss of continuity. 
The author, fortunately, considers program- 
ming as the entire problem preparation 
effort and not merely as the detailed writing 
of machine instructions. The various aspects 
of programming are, with one exception, 
well covered. Such topics as subroutines and 
indexing are well developed. A brief intro- 
duction to automatic coding is included. The 
illustrations and exercises are very effective. 
The organization of the book suffers in the 
separation of decimal point location meth- 
ods from floating decimal methods, but is 
otherwise successful. The text is written 
simply and is quite readable, despite occa- 
sionally faulty grammar or imprecise lan- 
guage. Although a good chapter on program 
checkout is included, the important com- 
panion subject of checks and reruns is not 
considered. Fuller discussion of the advan- 
tages and disadvantages of some of the 
various programming techniques presented 
would also represent an improvement, The 
book is well conceived, however, and gen- 
erally well executed. It isa good introduction 
to an increasingly popular subject. 
Peter Calingaert 


57-228 

Transistor Engineering Reference Hand- 
book—H. E. Marrows. (John F. Rider, Inc., 
New York, N. Y., 288 pp.-+-2 indexes-+-viil, 
illus.; 1956.) The author states in the intro- 
duction that, while many publications deal- 
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ing with transistors are available, “the need 
to assemble and coordinate information on 
all the commercial aspects of the industry 
under one cover had led to the writing of this 
book.” The volume contains five sections. 
Section I is a “General Survey of Tran- 
sistors.” Following a chronology of impor- 
tant transistor developments, a short de- 
scription of transistor materials, structures, 
and fabrication techniques, and a list of 
pertinent references is given. Characteris- 
tics and circuit properties of junction tran- 
sistors are discussed next. A classified bibli- 
ography on transistor applications concludes 
the section. The material is complemented 
by forty figures and tables related mainly to 
transistor circuit design. Section II occupies 
most of the volume and presents detailed 
data sheets (including characteristics and 
performance ratings) on some two hundred 
types of commercial transistors which were 
available at the time when the book was 
written. Section III contains reference data 
on circuit components used with transistors, 
such as audio transformers of various cate- 
gories, IF and pulse transformers, capaci- 
tors, batteries, thermistors. A list of com- 
mercially available transistor test sets is also 
given. Section IV lists the specifications and, 
in many cases, shows the circuit diagrams 
of some one hundred products utilizing 
transistors. Products described are various 
amplifiers, pulse generators and other oscil- 
lators, power converters, meters, flip-flop, 
communication equipment, and radio re- 
ceivers. The volume is completed by a Man- 
ufacturers’ Directory (Section V). This re- 
viewer is unable to display a great deal of 
enthusiasm for Section I. The material is too 
condensed to serve either as an understand- 
able introduction to transistor physics and 
applications for the beginner or as a usable 
reference for engineers acquainted with the 
art. As the section stands, not even the little 
space occupied by it is used efficiently. For 
example, the author describes such expert- 
mental devices as filamentary transistors, 
field-effect transistors, and fieldistors but 
but makes no mention of drift transistors 
and unijunction transistors. The model cir- 
cuits presented are not selected most ju- 
diciously to be representative of the con- 
- temporary art. There is some duplication: 
for example, instead of presenting a larger 
variety of flip-flops with a description of 
their respective merits, three essentially 
identical configurations are shown on pages 
1-25; similarly, the configurations of Figs. 
1-19 and 1-20 and those of Figs. 1-20A and 
1-20B are identical. The author did a very 
conscientious and useful job in compiling the 
material of Sections IT, III, and IV. These 
sections will be of considerable use to tran- 
sistor circuit designers. It certainly is prefer- 
able to have information on transistors and 
associated components available in a single 
volume rather than to have to search for the 
material in countless specification sheets as- 
sembled in a multiplicity of folders. It 
should, of course, be remembered that, due 
to the steady release of new and improved 
transistors and components, the value of the 
book is necessarily ephemeral. However, for 
the present, it will be found useful by many 
engineers. 
A. P. Stern 
Courtesy of Proc. IRE 
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An Introduction to Cybernetics—W. R. 
Ashby. (John Wiley & Sons, Inc., New York, 
N. Y., 287 pp.-+7 index--ix, illus.; 1956.) 
Dr. Ashby, who is known to many engineers 
as the author of Design for a Brain and the 
father of the “homeostat,” has written what 
he considers an easy tntroduction to cyber- 
netics and not “merely a chat about cy- 
bernetics.” The preface informs us that the 
book addresses itself to the “many workers in 
the biological sciences” who “are interested 
in cybernetics and would like to apply its 
methods and techniques to their own special- 
ty.” It is Ashby’s contention that the basic 
ideas of cybernetics can be treated without 
reference to electronics and that they are 
fundamentally simple. Ashby expresses, fur- 
thermore, the conviction that although ad- 
vanced techniques may be needed for 
advanced applications, a great deal can be 
done in the biological sciences by the use 
of quite simple techniques, provided they 
are used with a clear and deep understanding 
of the principles involved. He expresses the 
belief that by refining common-place and 
well-understood concepts step by step, the 
biological worker who has no knowledge of 
mathematics beyond elementary algebra 
can be introduced to such cybernetics 
topics as feedback, stability, regulation, ul- 
trastability, information, coding, noise, and 
so forth. Ashby divides his book into three 
major parts: Part I deals with the principles 
of mechanisms (including those that apply to 
determinate machines and very large sys- 
tems). Part II, entitled, “Variety,” deals 
with what is meant by information; it is de- 
signed to enable the reader “to proceed 
without difficulty to the study of Shannon’s 
own work.” Part III deals with regulation 
and control and provides “an explanation of 
the outstanding powers of regulation pos- 
sessed by the brain, as well as the principles 
by which designers may build machines of 
like power.” The book is well written and 
equally well produced (with the exception of 
a rather curious index which includes entries 
such as borrowed knowledge, chair chame- 
leon, coffee, fly-paper, flying saucer ghosts, 
Hitler, etc.). And yet, this reviewer cannot 
rid himself of a certain uneasiness in recom- 
mending it for the purpose for which Ashby 
wrote it, to wit: for those who want to 
achieve by self-study an actual mastery of 
the subject. Many admittedly easy exercises 
are scarcely designed to relieve this reluc- 
tance. When, in the chapter entitled, “The 
Black Box,” exercise 6/16-2 asks the ques- 
tion, “To what degree is the Rock of 
Gibraltar a model of the brain?,” it is with 
a sense of keen disappointment that one 
reads the answer, “It persists, so does the 
brain; they are isomorphic at the lowest 
level.” Such exercises hardly do justice to 
the often insightful exposition that Ashby 
provides elsewhere. The book’s strongest 
point is, perhaps, that it brings the mathe- 
matically unsophisticated reader into some 
elementary contact with a broad spectrum 
of mathematical raw materials for concep- 
tual models. In this respect Ashby’s book fills 
a void in the cybernetic literature in spite 
of its exceptionally scant bibliography. 
Ashby’s views of and aspirations for cyber- 
netics as a mathematical discipline are re- 
vealed when he writes, “the truths of cy- 
bernetics are not conditional on their being 
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derived from some other branch of science.” 
In Ashby’s view, cybernetics, whose subject 
matter is the domain of all possible ma- 
chines, can, like mathematical physics be 
indifferent to the criticism that it gives 
prominence to the study of nonexistent sys- 
tems. Here Ashby becomes a victim of his 
own semantics. The position of esteem that 
mathematical physics occupies today must 
in large part be attributed to the fact that 
throughout the past three centuries scien- 
tists have learned to abstract from their ex- 
perience aspects that are both mathemat- 
ically manipulable and yet relevant to pre- 
dicting and controlling an ever-widening 
range of physical phenomena. Workers in 
the “soft” sciences (it might be more ap- 
propriate to call them the not-so-vertebrate 
sciences: their problems are tough rather 
than soft, but the fields lack a strong theo- 
retical backbone) were deeply stirred by the 
cybernetic wave of the future. Here was a 
suggestion that they use new mathematical 
models and techniques that promised to be 
more applicable to problems of information, 
organization, perception, and perhaps even 
learning than the calculus has been. Here 
was a bridge to an era of computer technol- 
ogy with its potentialities of handling large 
quantities of data and of trying out complex 
conceptual models in reasonable time spans. 
Here was, finally, new hope for a rapproche- 
ment between scientists in these so-called 
soft areas with their more secure colleagues 
under the aegis of a possible unity of the 
communications sciences. Many of these 
hopes have been dashed, and much good 
will has been dissipated by those who 
honestly proclaimed that the problems of 
the biological and social sciences were on 
their way toward “solution” each time they 
had to their own satisfaction been able to 
transform “organized complexity” by the 
use of a more or less common-sense verbal 
device. This reviewer is convinced that we 
must not relax our efforts but that we must 
on the contrary dig in for a much longer 
pull. There are some tender shoots in areas 
in which models of modest scope and meas- 
urable phenomena are actually on speaking 
terms. We must intensify our explorations 
of mathematical models (such as in the the- 
ory of automata) and continue to search for 
significant experiments that will take ad- 
vantage of our newly acquired armamentar- 
ium. There are passages in Ashby’s book 
that are not in marked disagreement with 
the above views. And yet, the tone of the 
whole book does not seem to convey the 
proper perspective of the road ahead. There 
are scientists whose current work is more or 
less directly inspired by the writings of 
Wiener, Shannon, etc., Should their work 
be labeled as cybernetics and be introduced 
accordingly? This is certainly a matter of 
choice. Ashby has a perfectly good right to 
ignore these rather specific experiments and 
to substitute his own vision of the spectrum 
of problems that cybernetics will ultimately 
be capable of dealing with. But by not giving 
us at least one example of a truly successful 
application of a cybernetic model to a really 
tough biological problem, Ashby’s book 
claims too much and proves too little. 


W. A. Ro senblith 
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